K 7 3 23 376 

A 

NASA CR-1122'.O 


STC - SAB PROGRAM USERS MANUAL FOR 


THE TURBULENT BOUNDARY LAYER AND TURBULENT 
SEPARATION PREDICTION METHODS EMPLOYED IN 
THE NASA LANGLEY STREAMTUBE CURVATURE COM- 
PUTER PROGRAM 


CA COpy L£ 


BY: D.R. Ferguson 


Sponsored by NASA Langley, Contract NAS 1-10804, and Prepared 
by General Electric Company, Aircraft Engine Group, 
Cincinnati, Ohio 

for 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



1. Report No. 2. Government Acceuion No. 

NASA CR- 11 2240 


4. Title end Subtitle 

STC-SAB Program Users Manual for - The Turbulent Boundary Layer 
and Turbulent Separation Prediction Methods Employed in the NASA 
Langley Streamtube Curvature Computer Program. 


7. Author(s) 

D.R. Ferguson 


9. Performing Organisation Name and Address 

General Electric Company 
Aircraft Engine Group 
Cincinnati, Ohio 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 


3. Recipient's Catalog No. 


S. Report Date 

D 


6. Performing Organisation Code 


8. Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 
NAS 1-10804 



14. Sponsoring Agency Code 



16. Abstract 

The Streamtube Curvature Program (STC) has been developed to predict the inviscid flow 
field and the pressure distribution about nacelles at transonic speeds. The effects of 
boundary layer are to displace the inviscid flow and effectively change the body shape. 
Thus, the body shape must be corrected by the displacement thickness in order to calculate 
the correct pressure distribution. This report describes the coupling of the Stratford 
and Beavers boundary layer solution with the inviscid STC analysis so that all nacelle 
pressure forces, friction drag, and incipient separation may be predicted. The usage of 
the coupled STC-SAB computer program is outlined and the program input and output are de- 
fined. Included in this manual are descriptions of the principal boundary layer tables 
and other revisions to the STC program. The use of the viscous option is controllod by 
the engineer during program input definition. 


17. Key Words (Suggested by Authorfs)) 

Streamtube Curvature, Computer Program, 
Users Manual, Viscous Flow, Separation, 
Inlets, Nacelles, Transonic Flow 


19. Security Classif. (of this report) 



Unclassified 


20. Security Classif. (of this page) 

21. No. of Peges 

Unclassified 

155 



For sale by the National Technical Information Service, Springfield, Virginia 22151 






















TABLE OF CONTENTS 


PAGE 

1.0 INTRODUCTION 1 

2.0 PROGRAM DESCRIPTION 2 

2.1 Coupled STC-SAB Solution Procedure 2 

2.2 SAB Subroutines - STC Interface 3 

2.3 Program Nomenclature ... 3 

3.0 STC-SAB PROGRAM INPUT 9 

3.1 Overall Input Data 9 

3.2 Boundary Coordinates 10 

4.0 STC-SAB PROGRAM OUTPUT - SAMPLE CASES 11 

4.1 Standard SAB Boundary Layer Output 11 

4.2 Related Boundary Layer Output 13 

4.3 Sample Cases 13 

5.0 GENERAL OPERATING PROCEDURES 14 

6.0 REFERENCES 16 

7.0 FIGURES. . 17-78 

8.0 APPENDICES 

8.1 STC-SAB Input Sheets 79-89 

8.2 Program Listings . 90-155 



1.0 INTRODUCTION 


The Streamtube Curvature Program (STC) has been developed to pre- 
dict the inviscid flow field and the pressure distribution about nacelles 
at transonic speeds. While this program is a basic tool for the calcu- 
lation of the overall forces on nacelles, it is well understood that the 
effects of boundary layer friction drag and displacement of the inviscid 
flow must be included to obtain accurate performance predictions. The 
displacement of the inviscid flow effectively changes the body shape, 
thus altering the inviscid pressure distribution. The summation of pres- 
sure-area forces taken over the body can be seriously in error when this 
displacement effect is not included. 

The coupling of a boundary layer solution with the inviscid STC 
analysis allows inclusion of displacement effects as well as a complete 
evaluation of all nacelle forces, including friction, and the prediction 
of boundary layer separation. The turbulent boundary layer method selected 
for use with the STC program is the integral method of Stratford and Beavers 
described in detail in References 1 and 2. This report is concerned with 
usage of the coupled STC-SAB (Stratford and Beavers) program. Included in 
this manual are descriptions of the principal boundary layer tables, and 
the revised program input and output. 
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.0 I’ROCJtAM DESCRIPTION 


The Stroamtubc Curvature (STC) and the Stratford and Beavers 
Boundary Layer (SAB) routines have been combined as a composite program 
for evaluation of boundary layer displacement effects and friction dra;, 
as well as solution of the inviscid problem. The following sections 
deal with the coupled STC-SAB solution procedure, including a description 
of the principal SAB routines and the storage regions for input and cal- 
culated boundary layer data. 

2.1 Coupled STC-SAB Solution Procedure 

The STC program and its use are described in References 2 and 3 - 
The SAB section represents an extension of the current capabilities of 
the existing version of STC to include boundary layer effects, namely 
inviscid flow field displacement, friction losses, and turbulent boundary 
layer separation. The general calculation procedure for the coupled 
STC-SAB solution is as follows: 

1. ) Initially, the STC program is executed to obtain the pressure 

distribution on the solid boundaries at a given refinement or 
convergence level. 

2. ) Upon completion of the inviscid solution, the SAB routines are 

entered to calculate the boundary layer growth on specified solid 
surfaces and store the "least squares" smoothed displacement 
thickness information in the boundary layer tables. Boundary 
layer parameters such as momentum thickness, displacement thick- 
ness, and skin friction coefficient are printed immediately 
after the normal STC boundary output. If the boundary layer 
has separated, the location of the separation point is indicated. 

3. ) The STC solution may now be restarted using boundary information 

corrected for displacement effects. The adjusted inviscid pres- 
sure distribution is then used for another pass through the SAB 
boundary layer calculation routines. 

4 . ) Steps 2.) and 3 *) may be repeated until a desired level of con- 

vergence is attained. Note that no program control, such as 
weighting of consecutive & * values, is imposed on the coupled 
solution iteration. Continuation of the solution is entirely 
controlled by the user. 

Normally, three ( 3 ) to four ( 4 ) consecutive passes through the coupled 
STC-SAB program should provide adequate convergence of the adjusted in- 
viscid pressure distribution. 



2.2 SAB Subroutines - STC Interface 


lhe SAB control and calculation routines are located in the output 
overlay (2,2) and are called by subroutine WRIBDY after the inviscid 
boundary data are printed. The sequence of subroutine calls is as follows: 


1. ) Initially, WRIBDY calls LBDYBL (located in overlay (0,0)) to 

determine if a boundary layer is to be calculated for the given 
solid boundary. If yes, subroutine SAB is called with the bound- 
ary name given. 

2. ) Subroutine SAB determines the number of active points on the 

boundary using the information (surface distance SW) calculated 
in WRIBDY. Three types of situations may exist: 

TYPE = 1 Boundary layer initiated at a stagnation point SW(l). 

TYPE = 2 Axisymmetric spinner on the axis. Boundary layer is 
initiated at the last zero radius point, SW(NI). 

TYPE = 3 Boundary layer is started at a stagnation point down- 
stream of the start of the boundary. An initial 
"equivalent" flat plate distance must be supplied. 

The boundary layer history upstream of the stagnation 
point is ignored. 

3. ) Subroutine SABBL is called by SAB to calculate and print boundary 

layer information. The output quantities are described in detail 
in section 4. During the course of calculating the normal bound- 
ary layer parameters, the Stratford separation parameter F (de- 
scribed in Reference 2) is evaluated for all adverse pressure 
gradient regions. If the preceding parameter indicates a bound- 
ary layer separation, the calculation is terminated and the re- 
maining locations in the i* table are filled with the displace- 
ment thickness at the point preceding the separation point. The 
output in this case includes only the data prior to and including 
separation. 

4. ) Subroutine SAB finally calls BLTBBL to store the smoothed dis- 

placement thicknesses and their derivatives in the boundary layer 
data table. In this process, the flow adjustment table and the 
station table are moved up in memory to accommodate the boundary 
layer data. If enough room in the TABLES storage area is not 
available, a comment to this effect is printed and the boundary 
layer input table is altered to reflect no boundary layer presence. 
Upon restart, no displacement thickness correction will be applied 
to this boundary. Ihe above condition may be alleviated by pro- 
viding additional storage in the TABLES area, as described in 
Reference 4. 
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5.) The interface with STC is provided on restart through subroutine 
BDYPTM, which determines the location of orthogonals intersecting 
the solid boundaries. Interpolated corrections are applied to the 
boundaries using the following relations: 



Zc, Rc, tp 

c 


di 

*c ds 


R = R, - S cos ( <t> ) 

CD c 

Z = Z,+ i* sin ( <t> ) 

c b c 



R = R + 5 * cos ( <t> ) 

c b c 

Z = Z - S* sin ( <(> ) 

c b c 

If the boundary layer on the solid surface is separated, a comment 
to this effect is printed each time the boundary is accessed to 
calculate an orthogonal intersection. When this situation is en- 
countered, the user is advised to discontinue the calculation, 
since the displacement thickness information downstream of the sep- 
aration point is in error. 

Listings of the SAB subroutines and the STC routines altered for bound- 
ary layer capability are given in Appendix 8.1. 
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2. 3 I'rogram Nomenclature 


Communication between the STC routines and the SAB routines is ac- 
complished by the use of labeled common. The principal data storage >roas 
utilized by both STC and SAB are the boundary layer input table and e 
boundary layer data tables. These regions are described in detail in this 
section. The remaining labeled commons are given alphabetically following 
the description of the main boundary layer tables. Within each block, 
variables are listed according to the position occupied in the block. 
Erasable temporary storage is used for communication between WRIBDY and 
the boundary layer subroutines. Variables listed in this block are the 
names used in SABBL. In all cases, the pertinent dimension and type in- 
formation are included with the variable name (R = Real, I = Integer, 

L a Logical). Variables normally containing BCD data are typed as 
H s Hollerith, even though they may have real or integer names. 

Boundary Layer Input Table 

The boundary layer input table is stored in labeled common BLBDY. 
Input boundary layer information as supplied on page STC/Sheet-2 of the 
input sheets are read by routine RBD. Tie resulting table consists of 
the following three items stored serially for each boundary: 


Variable Name Type 


BLB(I) 


H 


BLB( I + 1) 


I 


BLB( I + 2) 


Description 
Boundary name 

Indicator designating whether a 
boundary layer calculation is to 
be performed. 

0 - No 

1 - Yes 

Initial equivalent flat plate 
distance to first point on boundary. 


The subscript I ranges from 1-58 and is incremented by 3 for each boundary. 
Information for a maximum of twenty (20) boundary layers may be stored; 
viz, 

Common /BDBDY/BLB ( 60 ) 


Boundary Layer Data Tables 

The framework of the STC program is designed to allow flexibility 
as to the configuration to be analyzed. For example, very weak limits 
are placed on the number of flow boundaries and the number of channels 
into which the flow is split. To allow flexibility, the bulk of the cal- 
culation data is saved in singly dimensioned arrays. Within each array, 
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tho data arc packed together to maximize storage efficiency. Information 
such as the boundary coordinates and flow properties is compactly stored 
in a single array, TABLES, so that only the total amount of information 
saved is limited by the array size. No limit is placed on the amount of 
information to be placed in any one Table of which there are currently 
six: 

o Channel input table 
o Boundary table 
o Table of convected properties 
o Table of wake displacement thickness 
o Flow adjustment table 
o Station table 

Boundary layer data in the STC-SAB program are stored in the TABLES 
region immediately before the flow adjustment table. A subtable for each 
boundary layer is constructed and the information is stored in the following 


order : 



Variable Name 

Type 

Description 

BNAME 

I 

Boundary name 

LBLNXT 

I 

Pointer to the next boundary 
layer table 

NSEP 

I 

Index pointing to separation 
location (normally 0) 

DUMMY 



SWREF 

R 

Reference distance for alteration 
of coordinates in the boundary layer 
table; viz, boundary origin 

SIGN 

R 

Boundary type 
-1 Upper boundary 
+1 Lower boundary 

SW(1) 

R 

Distance along surface 

DSTAR(l) 

R 

Smoothed displacement thickness ( 4 

DDS TAR ( 1 ) 

R 

Slope of smoothed displacement 
thickness (d4 /ds) 


SW(2) 
DSTAR(2) 
DDS TAR ( 2 ) 
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The boundary layer data table is located between LDO and LDE in 
the tables region. The index limits LDO and IDO are stored in /IXORIG/ 
after LSE. 


STC-SAB Labeled Commons 


Common 

Variable Name 

Tjr£e 

Description 

BLDTA 

EDNA ME 

I 

Boundary name 


LOWER 

L 

T = Lower boundary 
F = Upper boundary 


IBTYPE 

I 

Initial condition indicator 

1 = Boundary layer initiated 

at stagnation point SWBL(l) 

2 = Axisymmetric spinner 

3 = Boundary layer initiated 

at a mid-boundary stagnation 
point 


N1 

I 

Index of first point in SWBL 
table 


NI 

I 

Index of last point in SWBL 
table 


CAPX1 R Equivalent flat plate distance 

from boundary layer origin to 
first calculated boundary layer 
point. 


BLSEP 

NSLOC 

I 

Index of separation point in 
SWBL array (normally O) 

CGRAV 

CG 

R 

Gravitational conversion con- 
stant in consistent units 


ERASE 2 

DSTAR(IOO) 

R 

"Unsmoothed" displacement 
thickness & * 


SWBL (100) 

R 

Distance along surface 


ZW(100) 

R 

Axial coordinate (Z, X) 
on solid surface 


RW(100) 

R 

Radial (R) or normal (Y) 
coordinate on solid surface 
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Common 


Variable Name 


DSTR(IOO) 

DDSTR(lOO) 

VE(IOO) 

MACH ( 100 ) 

MACHSQ(lOO) 

CP(IOO) 


PQPT(IOO) 

PW(IOO) 

REXP(lOO) 

PR(IOO) 


CAPX ( 100 ) 


R 

R 

R 

R 

R 

R 


R 

R 

R 

R 


R 


"Smoothed” displacement 
thickness i* 

Derivative of "Smoothed" 
displacement thickness 
(d i* /ds) 

Velocity at edge of boundary 
layer 

Mach number at edge of boundary 
layer 

Square of MACH 
Pressure Coefficient 

c = 2(p-p )/yp M 2 M > . 1 

p s ' s o o 


C = O 
P 


M < .1 
o 


Static to total pressure ratio 


Static pressure on surface 
1-2 


R ** 


1.25 


used in boundary 


layer calculation 
Boundary layer parameter 


MACH 

1+ . 2 * MACH 


] 


4 

* REXP 


Equivalent flat plate distance 
along surface 


REBL RESTBL L 

VISCOS TREF R 

MU REF R 

SCON R 


Boundary layer restart indicator 
F - Normal option 
T - Restart - restore tables 

Reference temperature for dynamic 
viscosity calculation 

Viscosity at reference temperature 

Sutherland constant for fluid; 
used in viscosity calculation 


* 
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3.0 STC-SAB PROGRAM INPUT 


Program input to the SIC program is unchanged for the coupled STC-SAB 
version and is described in detail in Reference 3. The four (4) distinct 
card input sets read by the program are: 


1.) 

Input 

sheet 

0 

2.) 

Input 

sheet 

1 

3.) 

Input 

sheet 

2 

4.) 

Input 

sheet 

3 


Identification information 
Overall input data 
Boundary coordinates 
Channel flow properties 


The above input sheets, revised for the STC-SAB program, are given in 
Appendix 8.2. Since the revised sheets are essentially identical to the 
original SIC input sheets, only the changes applicable to the boundary layer 
calculation will be discussed. 


3.1 Overall Input Data 

The normal option exercised in the STC program is to specify pressures 
and temperatures in dimensionless form normalized by the free stream ambient 
pressure and temperature. When the SAB boundary layer option is chosen, 
however, pressures and temperatures must be given in dimensional form. 

Also, several additional parameters must be specified to define the fluid 
viscosity and thermodynamic quantities. 



Parameter 

Description 

Preset Value 

— 

TREF 

Reference temperature for 
viscosity calculation 

518.688 °R 


MUREF 

Reference viscosity at 
IREF 

-7 

10 lbm/in sec. 


SCON 

Sutherland Constant for air 

198.6 °R 

— 

CG 

Gravitational conversion 
constant 

32.174 ft-lbm/lb 


Within the program, viscosities are calculated using the following 
Sutherland relation: 


U - 


ref 




( 


T _ + S \ 
ref J 

T + S / 
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Since length units of inches are normally used, no attempt lias been made 
to preset necessary quantities to their metric MKS values. Tin; units used 
must l>e consistent, and the gas constant RG and necessary pressures (PS0 or 
PTO) and temperatures (TSO or TTO) must be given in the proper units. With 
the viscosity and CG specified as above, the proper values and units arc: 




Units * 

R, Z 

inches 

PTO, PSO 

lb f /in 2 

TTO, TSO 

0 

'r 

RG 

1716.2 ft lb f /lb m 

VMG1 


10^ ft/sec . 

VMG2 


looj ft/sec 


3 . 2 Boundary Coordinates 

Input to specify a boundary layer calculation is supplied in the 
$A namelist with the boundary coordinates. The preset option is to not 
calculate a boundary layer on the given surface. The necessary input to 
specify a boundary layer calculation is as follows: 


Parameter 


Description 


Preset Value 


BL 


F - No boundary layer 
T - Boundary layer 


F 


CAPX1 Equivalent flat plate distance 0 

from boundary layer origin to 
the first calculated boundary 
layer point 


The latter parameter CAPX1 may be calculated as follows: 


CAPX1 = — 

_ a 
Pr 



S . 
orig 


Pr a ds 


where a = 0 Planar flow 

a = 1 Axisymmetric flow 



1+.2M 2 


♦ 

See 


Table I for other sets of consistent units 
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.O STC-SAi> PROGRAM OUTPUT - SA^LE CASES 


The standard output from the STC program is divided into six sections: 

1. ) Card input and preliminary printout 

2. ) Input and calculated boundary coordinates and angles 

3. ) Solution history 

4. ) General input and output data 

5. ) Flow field data along orthogonal lines 

6. ) Calculated flow data along field boundaries and final channel 

momentum balances 

Output from the SAB portion of the STC-SAB program appears after the in- 
viscid boundary data for each solid surface with a boundary layer. The 
designation of BL = T or BL = F appears with the above printout of the 
"input and calculated boundary coordinates and angles". Additional related 
boundary layer output will be printed by BLTBBL if the TABLES area will not 
accommodate additional boundary layer data. Also, on restart, printout will 
be produced by BDYPTM if a separated boundary layer is present on the perti- 
nent boundary being accessed. 

4.1 Standard SAB Boundary Layer Output 

The standard output from the SAB portion of STC consists of the 
boundary layer parameters at each orthogonal intersection of the boundary. 
The initial output consists of a bold heading specifying BOUNDARY LAYER. 

This is followed by the boundary layer parameters. 


Variable 


Description 


XW 

THETA 


DSTAR 


DELTA 

REX 


Axial coordinate 
Momentum thickness 



Displacement thickness 



Boundary layer thickness i 
Local Reynolds Number R^ = 


PV S 
e e 
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Variable 


CAPX 


CF 


SW 

DSTR 

DDSTR 

SEP 


F 


Equivalent flat plate distance along 
surface 


S 



where P 


_ a 
Pr ds 



a = 0 Planar flow 
a = 1 Axi symmetric flow 


Skin friction coefficient 
T w 

C _ 

f 1/2 PV 2 
e e 


x = Shear stress at solid surface 
' w 


Distance along surface 

"Smoothed" displacement thickness b* 

* 

Derivative of DSHt (dA /ds) used for 
correcting the local flow angle 

Separation flag. Normally blank, 
appears as SEP if the boundary layer has 
separated 

Stratford separation parameter 


The Stratford separation parameter is defined as follows: 


dC 

1/2 

-6 „ 

8 E 


10 Re 

x 

ds 




0.1 



M 2 


M. 


2 


This parameter is only calculated in an adverse pressure gradient 

(dP /ds > 0) and appears as 0.0 in regions of favorable (dP /ds<0) or 
W 

zero pressure gradient. The distance along the surface S is taken as the 
distance from the beginning of the adverse pressure gradient region. 

Also, is the Mach number at the minimum pressure point. For practical 
purposes, separation is assumed to occur if F attains a value of 0.5 or 
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greater. Wien separation occurs, the boundary layer calculation is termi- 
nated and printout only appears up to the separation point. The integrated 
friction dray for the surface is printed below the tabular data and is not 
calculated f or a separated boundary layer. 

4.2 Related Boundary Layer Output 

The boundary layer data tables are built by subroutine BLTBBL at the 
end of the problem. If a boundary layer on a given surface is specified in 
the boundary layer input table and there is insufficient storage in the 
TABLKS region to accommodate all the data, the following output will occur: 

TABLE SPACE EXHAUSTED— BOUNDARY LAYER DATA 

FOR I BOUNDARY "boundary name" NOT SAVED 
I LOWER I 

The above comment serves as a warning to the user to increase the available 
size of the TABLES area. Following the printing of this comment, the cal- 
culation switch in the input table is turned off. 

If subroutine BDYPTM detects a separated boundary layer in the course 
of interpolation for the displacement correction, the following comment 
with boundary name and separation location appears each time the boundary 
is accessed: 

** WARNING** SEPARATED BOUNDARY LAYER, BOUNDARY = "boundary name" 
SW = XXXXX 


Since displacement thicknesses downstream of the separation point are in 
error, the user is advised to terminate the problem at this stage. 

Intermediate diagnostic printout of the boundary layer data table and 
the interpolated displacement data may be obtained by setting PDUM(l5) * 0, 
however, this printout is seldom useful to the normal user. 

4.3 Sample Cases 

Two sample cases have been selected to illustrate the input/output 
for the STC-SAB program. The first case consists of the unseparated 
boundary layers on the inlet and cowl surfaces of a NASA Inlet No. 8 oper- 
ating at M = .8 with a mass flow ratio of 0.8093- This example case is 
identical ?o the one given in the STC users manual (Reference 3 ) with 

exception of the presence of the boundary layers. The input and output for 
this problem is given in Figure? 1. Three (3) passes through each program 
are shown with MAXIT = 8. On the third pass, the STC tolerance was reduced 
from 1. to .001 to obtain more accuracy in the inviscid solution. 
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Tito second c.iso illustrates the boundary layer on a 2 -l> circular 
arc in a wind tunnel at M = . 663 . 11*0 input/outpul is given in 
Figure 2. As indicated in the output, boundary layer separation occurred 
at SW = 19 . 87 b. Upon restart, the previously described warning continent is 
printed each time the boundary is accessed. Figure 3 shows the comparison 
of the STC-SAB results with test data, and illustrates good agreement v Ln 
the separation location. Discrepancies in the pressure distribution may 
be accounted for by the presence of upstream boundary layer fences and 
bleed parts which caused the flow in this region to be three dimensional. 


5.0 GENERAL OPERATING PROCEDURES 


The STC-SAB program described herein may be run on any Control Data 
6400/6600 machine operating under SCOPE 3*0 or a higher level operating 
system. In general, operating procedures will differ from site to site. 

The following comments on deck set-up and operating instructions are re- 
stricted to the program as installed at the NASA Langley Research Center. 
Minimal changes should be necessary for successful installation at other 
CDC sites. 

A large version of STC-SAB resides in source and absolute binary 
form on permanent (data cell) disc storage. This version is capable of 
handling 768 field points and has a TABLES storage region 2200 locations 
in length. Should this space be inadequate, the length of this latter 
region may be changed using the procedures given in Reference 3- The 
installed version will execute from the absolute binary file using a field 
lengch of 106K. 

As indicated in Reference 4, a partially completed STC problem may 
be restarted using the output data from a previous STC execution. When 
used for a boundary layer run, the output tape also contains the boundary 
layer data from the previous run. If TAPOT and TAPIN are specified as T 
and 'CAPE 2 and/or TAPE1 are not assigned to a tape file, the system will 
assign them to disc. Hence, consecutive boundary layer restart cases may 
bo run by simply setting TAPIN and TAPOT to T on all cases after the first. 
The general procedure for carrying out a boundary layer iteration would be 
to run a given STC problem to a specified refinement/or convergence level 
and then run successive restart cases at the same level to converge the 
combined inviscid - boundary layer problem. For example, 
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BL 

itorat j on 
1 


BL 

iteration 

2 


w V V 

1 STC F T 

SA 

2 

MAXIT = 8, 

$ 

original STC input 


I I STC T T 

SA 

MAXIT = 8, 

It 


1 1 STC T 

$A 

MAXIT = 8, 

$ 
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Figure 


Comparison of Measured and Predicted Separation Points 


• Assumed DL thickness = 0.3 in. @ liX/C = -J . r > 

• M - U.Ch3 

i i 

• Hoference (AiAA 1’APfcH 7l“5&5) 


Measured 



Measured Pressure 

Distribution 
Predicted Pressure 
Distribution 



STREAMTUBE CURVATURE PROGRAM 
WITH BOUNDARY LAYER 

Oct. 1972 Page of 

Overall Input Data STC/Sheet-1 



input tape? output tape? 
T or F T or F 


1 STC 

Mach number, ambient pressure and temperature, fluid properties 

$A (1.) (1.) (1.) (1.4) 

MACH0= , TS0= , PS0= , RG= , GAM* 

Highlight radius, maximum body radius, body closure tolerance 

RHL= , RM= , TTE* , 

axisymmetric or planar? 

(T) or F 


AXI= 


spacial grid refinement criteria, see notes 

GR ( 1 ) = , , , , » » > 

SGR(1)= , , » * » j , 


NGR= 


GZ ( 1 ) = 


» 


SGZ (1 )= 


NGZ= , 

maximum Mach number increment between grid points 

streamwise normal 

direction direction 

( 0 . 1 ) ( 0 . 1 ) 


VMG1* , VMG2= 


maximum number of refinements 


MAXIT* 


fluid reference temperature, reference viscosity, Sutherland constant, 
units conversion constant 

(518. 7°R) (10~ 6 lbm/in.sec) (198. 6°R) (32.174 ft-lbm/lb f sec 2 ) 


TREF= , MUREF= , SC0N= , CG= 


$ 
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Oct. 1972 


STREAKTUBE CURVATURE PROGRAM 


Boundary Coordinates 


Page of 

STC/sheet-2 


boundary channel 
name name 

W 'W V 


2 BDY 


upper boundary? angle input? 
T or F T-no, F-yes 


$A 




PPER= 


ZR0NLY= 


boundary layer? equv. flat plate distance 
T-yes, (F-no) to boundary layer origin 

( 0 .) 

/ < ■. 

BL= , CAPXI® , 


Z R ANGD 


B ( 1 ) = 


$ 
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Feb. 1972 


STREAMTUBE CURVATURE PROGRAM 


Page of 

Channel Flow Properties STC/Sheet-3 


channel 

name 

W V 

3 CHN 

$A 

ratio of 
specific 
heats 
(1.4) 


gas 

constant 


( 1 . 0 ) 


flow rate 
may be 
adjusted? 
(T) or F 


GAM= , RG= , VARY= , 

stagnation properties, see notes 3 and 4 
total temp total pressure 

TT0= , PT0= , 

Mach no. static temp static pressure 


MACH0= , TS0= , PS0= 

flow area normalized by A^. . 


A 0= 


$ 


General Instruction and Notes for Sheet-1 of the STC Input Forms 


1) The STC Program computes the subsonic and transonic field of inviscid flow 
past (and within) arbitrarily shaped planar and axisymmetric bodies. Inlet 
and exhaust nozzle flows wherein there may exist jet streams with differing 
energies are typical applications. An optional boundary layer analysis 
(SAB) is included to evaluate friction losses and displacement of the 
inviscid flow. 

2) The total flow is composed of one or more streams, the properties of 
which are to be listed on Sheet-3 (except as noted below). Each stream 
occupies a "channel" which is identified by a one to six character 
alphanumeric word. Each channel must be bounded, at least in part, by 
an "upper boundary” and a "lower boundary". Each boundary is also given 
an identifying one to six character name and the coordinates are listed 
on Sheet-2. The following sketch illustrates the naming of channels 
(CHN) and boundaries (BDY) . 


L t I—U L ! C — L i 1 1 < J t 1 C 1 


T 


-A i tit t 4. « It 


CHN=EXT 


BDY=FF 



iz i 


CHN=W9 


An external flow channel must be named EXT, the recommended name for 
the inlet capture flow channel is INT , and the far-field interface 
boundary must be named FF. Otherwise the selection of the channel and 
boundary names is arbitrary. The special channel names EXT and INT cause 
extra streamlines to be placed in the first refined grid. The boundary 
name FF indicates that the boundary condition on FF is to be obtained 
from an analytic far-field solution. 

There is no specific limit to the complexity of the flow field in regard 
to the number of channels or the number of boundaries. Limits are set 
on the total amount of data which may be input. 
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The solution method consists of constructing a grid of streamlines and 
orthogonal lines. Starting with two streamlines per channel (one for 
each boundary) and an orthogonal passing through the first and last 
point of each boundary, the grid is automatically refined by dividing 
the grid intervals in half and in half again as required. The numerical 
resolution, the solution accuracy and the computer execution cost are 
all directly related to the extent of grid refinement. The input variable 
MAXIT determines the maximum number of refinements. Providing this limit 
is not exceeded, the grid will be refined, locally as required, until 
the spacing of orthogonals and streamlines is less than the value deter- 
mined from the SGR and SGZ tables and the Mach number difference between 
any two points on a streamline or an orthogonal line is less than VMG1 
and VMG2, respectively. Grid size values versus radius (or y-ordinate) 
are to be tabulated after SGR and GR, respectively. NGR is the number 
of entries in each list. Grid size versus the axial coordinate is to 
be tabulated after SGZ and GZ, respectively. NGZ is the number of GZ 
values. If dimensional values of RG, TSO and PSO are input (see Note 6), 
then VMG1 and VMG2 must have units of velocity rather than Mach number. 

See supplemental notes for additional details. 

A partially refined grid may be saved on tape by specifying a T in 
column 24 of the first card, or read from a previously created tape 
by specifying at T in column 14. 

If TAPE 1 and/or TAPE 2 are not assigned via a REQUEST card, they 
are assigned to disc. This allows the user to obtain output for a 
given refinement level and provides the option of changing input 
parameters on the restart. For the restart case, specify a T in 
column 14 of the first data card and include in the $A list only 
those input quantities (viz; MAXIT) which differ from those originally 
input . 

In the initial calculation grid, an orthogonal line will pass through 
each leading and trailing edge point and through each sharp corner point 
(i.e. a point on the boundary with an angle discontinuity). It is not 
possible to analyze a configuration in which two or more of these points 
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are approximately opposite to each other. For example, if a configura- 
tion contains more than one leading edge, the edges must be staggered 
relative to the streamwise direction. 

5. A free stream Mach number is specified by supplying a value of MACH0. 

6. Perfect gas assumptions are employed and the levels of ambient pressure 
and temperature may be dimensionless (TS0=PS0=RG=1) or dimensional. 
Dimensional values in a consistent set of units as described in the 
STC-SAB User Manual must be supplied if boundary layer calculations are 
requested. 

7. A reference (or highlight) area is calculated from the input value of 
Rut as follows: 

= 1 tR hl 2 


axisymmetric: 




pianar: A^ - Ay RL - R^ 

This reference area (or Ay in the planar case) is used for defining the 
mass flow for each channel. See STC/Sheet-3 note 5. 

8. Computed pressure drag forces are normalized by the (maximum) body area 
where 

axisymetric: Am ■ tt R ^ 

planar: A * Ay » R 

m m m 

9. Finite trailing edge thickness is permitted; the maximum thickness, 
or body closure tolerance, is to be supplied after TTE. 

10. On this and the following input sheets, the values in parenthesis are 
used if other values are not input. 

11. For boundary layer cases, supply a reference temperature, reference 
viscosity and a Sutherland constant if different from air values. The 
following sets of units may be used: 
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Units 


Parame ter 

Dimensionless (STC) 

English (in. ) 

English (ft . ) 

MKS 

L 

any 

in. 

ft. 

M 

PSO.PTO 

*atm 

psla 

psfa 

N/M 2 

TSO.TTO 

*atm 

°R 

°R 

°K 

TREF 

- 

°R 

°R 

°K 

MUREF 

- 

lbm/in.sec 

Ibm/ft.sec 

kg/m. sec 

SCON 

- 

°R 

°R 

°K 

RG 

1 

ft 2 /sec 2 °R 

ft 2 /sec 2o R 

J/kg °K 

CG 

- 

ft-l.b/lb^sec 2 

ft-lbm/lb^sec 

2 (unity) 

VMG1 ,VMG2 

** 

ft/sec 

ft/sec 

m/ sec 

*atm 

- Normalized by ambient conditions 




** - Dimensionless (values are approximately equal to a 

Mach number difference) 

Notes for Sheet-2 of the STC Input Forms 

1. Use one of these sheets for each boundary. Supply a one to six character 
name to identify the boundary in column 14 of the first card. Also 
indicate the name of the channel to which the boundary is adjacent in 
column 24. On the second card indicate whether the boundary is above 
(UPPER = T) or below (UPPER = F) the channel. 

2. The upper or lower "contour" which bounds a given stream may be composed 
of several "boundaries". In this case, an input sheet must be completed 
for each boundary; the last point of the first boundary must have the 
same coordinates as the first point of the second boundary, and so forth. 

This option is useful when considering variable geometry configurations 
such as flaps or movable nozzle parts. The movable part may be translated 
and rotated, as indicated by Note 8, while the fixed part is held stationary. 

3. List values of Z (or X) , R (or Y) and the surface angle in degrees at 
discrete points along this boundary contour after the symbol "B(l)=". 

Points at sharp corners must be listed twice, one time for each angle 
which exists at that point. In each interval, the STC Program fits 

a locally rotated cubic polynomial. The input points must be smooth 
and consistent with the specified angles. 
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All points are to be listed In the streetwise direction. For an inlet 
lip, the points are listed by starting at the highlight point and then 
proceeding around the nose to the trailing edge or downstream boundary. 

The internal and external surfaces are listed separately under different 
boundary names. However, the coordinates of the first point must be the 
same with ANGD equal to +90° for the external surface and -90° for 
the internal surface. 

It is recommended that the boundary coordinates and angles be obtained 
from an analytic definition of the contour, and that around the nose, 
angle variations between points be 20° or less. 

4. Pressure and Mach number distribution data will be printed at each 
orthogonal intersection with the boundary, and not at each input boundary 
point. Orthogonal stations, however, will be placed at any repeated 
point in the boundary table. List the same points twice if it is desired 
to have an orthogonal placed in that position. (This option is modified 
slightly when ZRONLY = T.) Orthogonal stations are always placed at the 
beginning and end of a boundary and at a juncture point between boundaries 
along the same contour. 

5. If the coordinates but not the angles are known, the third column in the 
B-input array may be omitted. In this case specify ZRONLY - T and list the 
coordinates twice at any point where a curvature jump or an angle jump 
exists. The double points will later be deleted if the angle discontinuity 
is less than 0.01 degrees. These double points are removed because extra 
calculation stations (see Note 5) are usually not desired at such points. 
However, the double point angle tolerance, DBLPTS, preset as 0.01, may be 
input as zero if such double points are to be retained. 

6. With either input option, care must be taken to specify the coordinates 

with precision. The round off or reading error of the coordinate data 

2 

should be less than AS / (10*L) , where AS is the local distance between 

points and L is some characteristic length, say the length of the cowl. 

1/2 

Conversely, the spacing between points should not be less than (10 6L) 
where 6 is the relative accuracy of the coordinate data. The tabulated 
output curvatures may be consulted to verify the smoothness of the input 
data. 
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7 . 


NACA Series 1 Cowl coordinates are stored internally. With the ZRONLY - T 
option they may be selected by listing: 

B(l) * 991, 1, 

X^, Y^, Series 1 Segment 

Y 2* Y 2’ 

Xj t ^2 * 

, , Cowl Aft Segment 

where X^, Y are the highlight coordinates and X 2 , Y^ is the position 
of the maximum diameter at the end of the Series 1 contour segment. 

8. The input coordinates of a boundary may be adjusted by supplying the 
following input quantities not shown on the front of this sheet: 


ROTATE 

ZPIVOT 

RPIVOT 

SCALE 

ZTRANS 

RTRANS 


angular rotation in degrees 
pivot coordinates 

multiplicative constant to be applied to the coordinate data 

translation increment in the axial direction 

translation on increment in the radial/vertical direction 


The order of transformation is rotation, scaling and translation. Hence, 
the pivot coordinates are in the same coordinate frame as the input data 
and the translation increments are in the rotated coordinate frame after 
scaling. It is only necessary to input data for the transformation 
operations to be executed. 


9. The normal option assumes no boundary layer (BL » F) . If a boundary 
layer calculation is desired, input BL ■ T. Also, supply an initial 
"equivalent" flat plate distance for the boundary layer origin to the 
first calculation sation if different from 0. (stagnation point) 

Si 


CAPXI 


+1 


“1% 



RP a dS 


orig 




1.25 

1.2 


ReXj 

Rex, 


10 

10' 
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Notes for Sheet-3 of the STC Input Forms 


Use one of these sheets for each channel to supply entrance flow properties. 
(See exception under Note 5). 

Of the input items shown on the face of the input sheet, use only those 
which are required for the selected options. 

The total pressure and total temperature may be input by either of the 
following two procedures: 

a) Specify TTO and PTO if the stagnation properties are known. 

These values may be normalized by the free stream ambient 
temperature and pressure. 

b) Specify MACHO, TSO and PSO if the static properties and Mach 
number are known. Again TSO and PSO may be normalized by the 
free stream ambient values. If only MACHO is supplied (TSO & 

PSO are omitted) the TSO and PSO values from STC/Sheet-1 will 
be used. 

If neither of the above is input, free stream values as supplied on 
Sheet-1 are used for MACHO, PTO and TTO. 

If the gas constant, RG, is different from the value supplied on 
STC/Sheet-1, supply the value which applies to this channel. RG, TS0 
TTO, PSO and PTO may all be given as dimensionless (normalized by free 
stream ambient) or dimensional using a conslstant set of units. 

Input a value AQ for the determination of the channel flow rate. AO is 
an area fraction mormalized by A^ 88 defined under Note 7 of Sheet-1; 
the (dimensional) channel flow area is then the product AO*A^ . The 
flow rate for the channel is computed by using one-dimensional relations 
from the total properties (as determined under Notes 3 & 4) , the 
supplied Mach number, MACH0, and the flow area. For internal inlet 
channels, specify R^ L as the highlight radius and AO as the mass flow 
ratio. 




5. If for any channel the Input data on this sheet is not supplied, the 
reference properties on STC/Sheet-1 will be employed and the frontal 
area calculated at the entrance station will be used. This option is 
suggested for an external stream. 

7. Although approximate flow rates must always be supplied according to 
Note 5, the program will adjust channel flow rates as required 
to meet the zero pressure loading conditions at a trailing edge or to 
meet a maximum (choked) flow rate. The number of channels which require 
flow rate adjustment is equal to the number of trailing edges. If the 
flow rate is not to be varied for this channel, specify VARY = F. 



PROGRAM LISTINGS 



♦DECK BL6L0K 

BLOCK DATA RLBLOK 

COMMON /1XOR1G/ I DUM 1 ( 14 » , LOO , LDE » IDUM2 ( 17 ) 
COMMON /BIBDY / IBLB(60) 

COMMON /VISCOS/ TREE , MUREF , SCON 
KLAL MUREF 

COMMON /Kt BL / RtSTBL 
LOGICAL RE STBL 

DATA IBLB/60*0/ 

DATA TREF, MUREF, SCON /5 18 . 688 , 10 .E-7, 201 .6/ 
DATA LDO, LDE/ 1,0/ 

DATA RESTBL/F/ 

ENID 


♦ DECK LBDYBL 


FUNCTION LBDYRL1BNAME, LOWER) 

CLBDYBL LOCATE INDEX IN bL INPUT TABLE 

LOGICAL LOWER 


INTEGER 

COMMON /CHDATA/ 

* 

* 

LOGICAL 
I NIT EGcR 
DIMENSION 

equivalence 

COMMON /CBITS 


BNAME 

HOT ( 1 ) » L B NE X T (1 ) » L B Z 1 ( 1 ) * 

CHNAME ( 1 ) tUPI 1) .LEDEX(l) , 

Z B T ( 1 ) ,RBT( 1 ) ,ANGB T ( A2 ) 

UP 

BDTtCHNAME,BDNAME 

BDNAME ( 1 ) » LBA (1) ♦LB8I1) 

( 8DNAME » ZB T ) » ( LB A, RBT > . < LBB , ANGBT ) 
/ I BI TS » I BLANK 


C LBDY BL =0 IF NO BOUNDARY LAYER 

C L BDY BL = I NDEX OF BOUNDARY IN BL INPUT TABLE 


COMMON /BL BDY / BLBI60) 

DIMENSION IBLBI60) 

EQUIVALENCE ( I BLB , BLB ) 

C CHECK FOR LOWE R= T — ( POSSIBLE COLLATED BOUNDARY) 

LBDYBL = 0 

I F ( LUWER ) GO TO 100 
1 I BL = -2 

I BL = IBL + 3 

I F ( IBLB( IBL ) .EQ. BNAME ) GO TO 10 

IF( IBLB( IBL) .EQ.IBITS .OR. I8L.GE.58 ) CALL ERRORK ( 6HLBDYBL ) 

GO TO 2 

10 I F ( I BLR( IBL + 1 ) .EQ.O ) GO TO 12 
LBDYBL = IBL 
12 RETURN 

C LOCATE POSITION IN BOUNDARY TABLE 
1JJ LB = LBF(BNAME) 

I F ( LEDEX(LB) .EQ.O ) GO TO 1 
LBZ = LBZi(LB)+3+LB 

IB = BDT(LBZ) 

1 1 J IBL = -2 

120 IBL = IBL + 3 

I F ( IBLBI IBL ) .EQ.IB ) GO TO 10 

I F ( I BLBI IBL ) .EQ. IBITS .OR. IBL.GE.58 ) GO TO IAO 

GO TO 120 

1 AO CALL ERRORK ( 6HLBDYBL ) 


END 
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♦DECK SAB 

SUBROUTINE SAB( ENTRY ) 

CSAB MAIN SUBROUTINE FOR BOUNDARY LAYER CALCULATION 

INTfcGER ENTRY 


C ON ENTRY=F I RST , SAVE B,S2 ON TAPE4 
C ON ENT RY =L AST, RESTORE B,S2 


* 

* 

♦ 

♦ 


COMMON /BLDTA / 

INTEGER 

LOGICAL 

COMMON /CB / 

COMMON /CS2 / 

COMMON /IXORIG/ 


COMMON /ERASE2/ 

COMMON /ALLCOM/ 
LOGICAL 


BDNAME, LOWER* IBTYPE»NI»NI ,CAPX1 
BDNAME 

LOWER 

B ( 300 ) 

S2I300) 

LHO, LHE, LBTO ,LBDE , LTOtLTE, LWO,LWE, LFO 
LOtLESTA, LSO , LSE , LDO,LOE , LOUM ( 4 ) * 
MO»NM»NJ , NFCOLSfMAXNJ,MAXOL , MAXNM,MAXLE, 
LEO, LEE, LRO, LRE , LRD 

XIK 100) ,SW( 100 ) ,ZW( 100) ,R W ( 100) ,DUM(200) 
Vfc ( 100 ) , 0UM1 ( 800) 

DUM2I5), AXIA,DUM3( 14) 

AXIA 


GO TO (1,2,45) , ENTRY 

1 REWIND 4 

WRITE (4) ( B C I ) , 1=1, NM) , <S2( I) , I*1,NM) 

C SCAN TABLES TO SET N1 

2 I BTYPE= 1 

IF( KWCD.EQ.O .AND. AXIA ) IBTYPE-2 
GO TO (5,8) , IBTYPE 

5 DO 6 I = 1 , N I 

I F ( V E ( I ) .EQ.O. ) GO TO 20 

6 CONTINUE 

I BTYP E= 3 

N 1 = 1 

GO TO 30 
« DO 1 j I = 1 , N I 

IF ( Rw( I ) .NE.O. ) GO TO 12 
10 CONTINUE 
RETURN 

12 N 1 = I-i 

GO TO 30 
20 N1 = I 

C CALCULATE BL FOR BOUNDARY — (BDNAME) 

30 CALL SABBL 

C INSERT SMOOTHED DATA INTO /BL TAB/ 

40 CALL BLTBBL 
GO TO 50 
45 REWIND 4 

READ (4) ( B ( I ) , I = 1 , NM ),(S2(I),I s l »NM ) 


50 RETURN 



*OLCK SABBL 

SIJRKUUT INK S ARAL 

* S ABBL 


COMMON /CB1TS / BITS, BLANK 
EQUIVALENCE (HITS, IBITS) * ( BLANK , I BLANK ) 


COMMON /ALLCOM/ 
REAL 
LOGICAL 

COMMON /BLDTA / 
LOGICAL 

COMMON /ERASE2/ 

1 

2 

DIMENSION 

EQUIVALENCE 

REAL 

COMMON /VISCOS/ 
REAL 


MACHA,PSA,TSA,PTA,TTA,AXI , RG , GAM , MACHC ( 12 I 
MACH A , MACHO, MAC HOS 
A X I 

BDNAME .LOWER, IB TYPE, Ni ,NI , CAP XI 
LOWER 

DSTAK( 100 ) ,SW(100) ,ZW(100) ,RW<100) ,DSTR( 100) , 
DDSTR ( 100 ) , VE ( 100) ,MACH( 100) .MACHSQI 100 ) ,CP( 100 ) , 
PQPTI 100 ) ,PW( 100) , REXP(IOO) , PR ( 1 00 ) ,CAPX< 100) 

XW(1 ) ,YW( l) 

( ZW,XW) , ( RW , Y W ) 

MACH , MAC HSQ 
TREE .MUREF, SCON 
MUREF 


COMMON /CGRAV/ CG 


COMMON /BL SEP / NSLOC 


DIMENSION 


3 

EQUIVALENCE 


REX( 100) , THE TA( 100 ) , DELTA ( 100 ) ,P ( 100) ,F1< 100) , 
F2(10j) ,F3( 1 JO) tCF(lOO) , I SEP ( 100 ) , DCPQDX ( 100 ) , 
F ( 100) ,AVG( 100) 

,CPK( 100) ,DC PK ( 100) 

( DCPQDX, DCPK) 


DATA PI/3. 1^159/ 
DATA KSEP/3HSEP/ 


*A* 

NSLOC = 0 
N2 = Nl+1 

NT = NI-N1+1 

PSO = PSA 

TSO = TSA 

MACHO = MACHA 

I F < MACHO.EQ.BITS > MACHQ=MACH(N1 ) 

I F < MACHO. EQ.O. ) MACHO=MACH I N2 ) 

MACHO S = MACHO*MACHO 
1B0 GAM1 = GAM/ I GAM- 1 . ) 

CAPX2 = 0. 

IF< CAPX1.NE.0. ) CAPX2 = CA PX1 
CVP = RG/ ( GAM- 1 . ) 

I F ( TTA.EQ.BITS .OR. TTA.EQ.l. ) TTO=TSO*( 1 • + . 5* ( GAM-1 • )♦ MACH 0**2) 
1 F ( PTA.EQ.BITS .OR. PTA.EQ.l. ) PTO=PSO* ( TTO/TSO ) **GAM1 
VMAX = SORT ( 2 .*GAM1*RG*TT0 ) 

PTOQPU= ( 1.+(GAM-1.)*.5* MACHO*MACHO ) **GAM 1 

CPT = . 5*GAM*MACH0*MAC HO 

DO 190 I =N 1 , N I 

MACHS Q( I) =MACH( I )*MACH< I ) 

PQPT( I ) = (CP( I ) *CP T ♦ 1 . ) /P TOQPO 
PWII) = PTO*I l.-(VE(I»/VMAX)**2)**GAMl 
190 CONTINUE 
* CALCULATE EXP 

RHOT = PTO/ ( RG*TTO) *CG 
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GAMM = 1 • ♦ ( GAM- 1 . )* . 5*M ACH0*MACH0 
RHOS = RHOT*GAM**(-(l./(GAM-i.Hl 
TSO = TTO*GAMM 
V = MACHO*SQRT (GAM*RG*TSO) 

AMU = MUREF*( TSO/TREF) **l.5*( TREF+SCONI/ C TSO+SCON) 
AL = ( SW(NI >-SW(Nl) )/2. 

RE * RHOS*V*AL/AMU 

EXP = 1.25 

IF( RE.GT.2.E7) EXP=1.2 

IF(fcXP.EQ.1.25) GO TO 205 

CONI = .23 

C0N2 = .022 

CON 3 = .028 

CON A = -( 1./6. ) 

GO TO 210 
2u5 CONI = .37 
CON 2 = .036 

C0N3 = .046 
C0N4 = -.2 

210 I F ( .NOT .AX I ) EXP=0. 

DO 215 I =N1 »N I 
REXP( I ) =0 . 

I F C .NOT. AX I ) REXP(I) = 1. 

IF(RW( IKGT.O. ) REXP(I)sRWm**EXP 

PR ( 1 ) = (MACH(I)/( l.+MACHSQ( I )* . 2 )) **4*RE XP < I> 

215 CONTINUE 

*B* CALCULATE SW,CAPX,REX 
GAM 12 = (GAM-1.)*. 5 

AMU = MUREF* ( TTO/TREF ) **1.5* ( TREF+ SCON) / (TTO+SCON) 
Z2 = SORT ( GAM/ ( (GAM-1. )*CVP<‘TTO) I 

GAMP = ( GAM-2. )/ (GAM-1. ) 

Z4A = SCON/TTO 

Z4D = 1 . / ( 1 . + Z4A ) 

Z1M = PTO*CG/ AMU 

CAPX( N1 )=CAPX2 

CALL SETM ( 1 » I BLANK » I SEP* 100) 

DO 220 N=N2 » N I 

I = N-l 

SWD = SW(N)-SW( I ) 

AINT = (PR(N)+PR( I) )*.5*SW0 

C APX ( N ) = AINT/PR(N)+CAPX( I)*PR( I I /PR < NJ+CAPX2 
TTOT = l.+GAM12*MACHSQ(N) 

Z1 = M ACH ( N ) *Z 1M 

Z 3 = TTOT**GAMP 

Z4 = ( 1./TT0T+Z4A )*Z4D 

REX(N)=Z2*Z1*Z3*Z4 
220 CONTINUE 

CALL LSPFIT(SW(Nl)»PW(Nl),NTt SW( N1 ) t F 3 ( N1 ) ,NT , 1 ) 

TTOT = l.+GAM12*MACHSQ(Ni) 

Z1 = MACH(Nl) 

Z 3 = TTOT**GAMP 

Z4 = ( l./TT0T+Z4A)*Z4D 

REX(N1)=Z2*Z1*Z3*Z4 
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*C* CALCULATE THETA, DSTAR, DELTA 
K/. - 0 

CALL SETM< 1 ,0. , E , 100 ) 

T H L T A ( M 1 ) = 0 . 

DSTAK ( N1 ) =0 . 

Of. LTA(N1)=U. 

FIND = 0. 

FMAX = — 10 • **6 

DO 2 JO I =N 1 » N 1 

I F ( I.NE.Nl ) GO TO 225 

IF ( CAPX2.EQ.0. ) GO TO 230 

225 CAPXX = CAPX( I ) * ( REX ( I ) *C APX ( I) ) **C0N4 

THETA ( I ) =C0N2* ( ( l.+MACHSQU>*.n**(-.7))*CAPXX 
DS T AK ( I ) = CON 3* ( l.+MACHSQ( I) * • 8 ) ** ( .44) *C AP XX 
DELTA ( I )=C0N1*CAPXX 
IF ( I . cQ.Nl ) GO TO 230 

* CHCCK FOR SEPARATION 

I F ( P*i{ 1 + 1 ) .LE.PWl I ) .OR. I.LE.K2 ) GO TO 2290 
K * I 

1225 K * K + l 

I F ( K.GT.NI ) GO 10 1226 
I F ( PW(K) .GT.PW(K-l) ) GO TO 1225 

1226 K 1 = 1 

I F ( Kl.EQ.(K-l) ) GO TO 2290 
K 2 = K-l 

KIM = Kl-1 

I F ( MACH(KIM) .EO.O. ) K1M=K1 
MACHOS = MACHSQl K1 ) 

CPK ( K 1M ) = i. 

DO 226 K=K1M,K2 
I F { MACH( K ) .EQ.O. ) GO TO 226 
C PK ( K ) = 1.-MACHS0(K ) /MACHOS 

226 CONTINUE 

DO 2^7 K=K 1 »K 2 

227 DCPKl K ) = ( CPK( K )-CPK( K-l ) i/<SW(K) — SWI K— 1 ) > 

K 2 M = K 2-1 

DO 226 K=K1,K2M 

2 26 DC PK ( K ) = ( DCPK ( K ) +DCPK( K+ 1 ))*. 5 
DO 229 K=K1,K2 
SWK = SW(K I -S W ( I J + CAPX ( I > 

F(K) = CPK(K)*(SQRT(ABS(SWK*DCPK(K> ))*<( 1 E-6 ) *RE X I K ) *SWK )**(-. 1 > ) 

229 CONTINUE 

2290 FMAX = AM AX 1 ( F ( I ) , FMA X ) 

I F ( FMAX.GE. .5 ) GO TO 232 

230 CONTINUE 

N3 = I 

GO TO 2 34 

* SEPARATION 

232 I SEP( I )=KSEP 
N3 = I 

NSLOC = N 3 
I = NI-N3 

CALL SETMl 1,DSTAR (N3) , DS T AR ( N3+ 1 ) * I ) 

* D* CALCULATE P FOR TOD 
234 P ( N 1 ) = 0. 

DO 240 I =N2 * N 3 
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K = 1-1 

A1 = ( RW ( K ) +RW ( I ) )* *5 

A 2 a (DSTAR(K)+DSTAR| im.5 

I F ( AX I ) GO TO 235 

P(I) * A2«(PW( I )-PW(K) )*PIK) 

GO TO 240 

235 P ( I I * 2.*PI*A1*A2*(PW( I )-PW(K) )*P(K) 

240 CONTINUE 

* CALCULATE TOD, TOTAL SKIN FRICTION DRAG 

IF( AX I) GO TO 250 

DRM = GAM*( (PW(NI )*MACHSQ(NI )*THETA(Nin-(PW(Nl)*MACHSQ(NI) 

I *THET A( N1 ) ) ) 

GO TO 255 

250 DRM = GAM*( ( P W ( N I ) *MAC HSQ( N I) *THETA ( Nil *2. *PI*RW < NI II- 
; (PW(N1)*MACHSQ(N1)*THETA(N1)*2.*PI*RW(N1>)) 

255 TOD = DRM-P ( N I ) 

♦E* CALCULATE CF 
300 DO 310 I=N1,N3 
RX= 1 . 

IF( AX I ) RX = RW ( I ) 

F 1 ( I ) = RX*PW( I )*MACHSQ(I) 

F2 ( I ) = F 1 ( I )*THETA( I I 
310 CONTINUE 

NN = N3-N1+I 

CALL L SPF I T ( SW ( Nl ) , F 2 ( Nl ) , NN , SW ( Nl ) ,CF < Nl ) , NN , 1 ) 

N 1 1 =N i 

I F ( MACH ( 1 ) .NE. 0.) GO TO 319 
Nl 1=N2 
CF ( 1 ) =0 . 

319 DO 320 I =N 1 1 » N3 

CF C I ) = 2.*CF(n/Fl(I)-2.*DSTARm*F3(n/IGAM*PWm^MACHSQU)) 

320 CONTINUE 

CALL LESTSQ(SW,DSTAR,Nl,NIt3»5tDSTR) 

NN = NI 

DO 327 I =N1 , N I 

327 DSTR(I)= OSTR ( I ) +DSTAR ( I I 

CALL LSPFIT(SW(N1)*DSTR(N1) *NN« SW( Nl I »DDSTR(N1 ) * NN, 1 ) 

* WRITE OUTPUT 

WRITE (6,1002) 

1002 FORMAT ( //37X, 30H BOUNOARY LAYE R//) 

WRITE (6,1004) ( I ,XW( I ) » THETA (I ) ,DSTAR( I ) , DELTA ( I ) ,REX( I ) , 

1 CAPX( I ) ,CF(I ),SW( I),DSTR(I )« DDSTR (I)»ISEP(I)*F(I),I*N1, N3 I 

1004 FORMAT (4X,1 HI , 5X, 2HXW , 4X . 5HTHET A, 5X, 5H0STAR,4X»5 HOELT A* 5X, 3HREX , 

* 7X,4HCAPX,6X,2HCF,8X,2HSW,6X,4HDSTRt4X,5H00STR # 5X,3HSEPt8X* 

* 4HFSEP/ 

* (2X, I3,F9.4,3F9.5,F9.0,F9.4,F9.5fF10.4*2F9.5 t 2X,A6,F13.6). ) 

WRITE (6,1003) TOD 

10 j 3 FORMAT (/6X,20HT0TAL FRICTION DRAG*»Fi4«5) 

900 RETURN 
END 
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♦DECK BLTBBL 

SUBROUTINE BLTBBL 
CBLIBBL BUILD BOUNDARY 


LAYER TABLES 


COMMON /BLDTA / 

INTEGER 

LOGICAL 

COMMON /IXQRIG/ 

* 

♦ 

* 

STATION TABLE 
COMMON /CHDATA/ 

* 

COMMON /ERAS E2/ 

♦ 

DIMENSION 
EQU IV ALENCE 
BOUNDARY LAYER 


BDNAME, LOWER, IBTYPE,N1,NI ,CAPXl 
BDNAME 

LOWER 

L HO, LHE , LBTOtLBDEt LTO»LTE « LWO»LWE» LFO.LFE, 
LO,LESTA, LSO«LSE , LDO »LDE , LDUMU), 

MO,NM,NJ »NFCOLS,MAXNJ,MAXOL *MAXNM »MAXLE» 
LEO,LEE»LRO*LRE»LRD 


X1(1),LNEXT(1),MLB(1) »MUB ( 1 ) ,PRI M( l ) , 

TYPELB ( l ) tNAMELBI 1)*ILB(1)»FLB(1) 

XI 1( 100) tSWBL(lOO) fZW(lOO) »RW<100) fDSTR(lOO) f 
DDSTRI 100), VE( 100), MACH (100) ,DUM(700) 

LEDEXI 1) *LBZ1 (1 ) 

I LEDEX,TYPELB)» (LBZltMLB) 

TABLE 


INDEX- LD=LDO* LDE INITIALLY 1,0 


DIMENSION 

INTEGER 

EQUIVALENCE 


* 

* 


COMMON /BL SEP / 

COMMON /REBL / 

LOGICAL 

COMMON /CPRINT/ 

COMMON /CTABPR/ 

COMMON /BLBDY / 

INTEGER 
LOGICAL ENTRY1 
DATA ENTRY 1 /T/ 

DATA UPPER tLOWR/SHUP PER* 5HL0WER/ 
DATA SWSAVE/O./ 

I F ( RESTBL ) GO TO 1111 

GO TO 1 


BNAME(l) ,LBLNXT(l) ,NSEP(2) ,SWREF( 1) ,SIGN(1), 
SW ( 1 ) • DS TAR ( l)fDDSTAR(l) 

BNAME 

( BNAME , X 1 ) * ( LBLNXT ,LNEXT ) , ( NSEP* MLB) , 

( SWREF,PRIM) , (SIGN, TYPELB), (SW,NAMELB), 
(DSTAR,ILB), (DDSTAR,FLB) 

NSLOC 

RESTBL 

RESTBL 

PDUMM( 6) , PDUM(20) 

I ITAB 
I BLB ( 60 ) 

UPPER 


- C 


RESTORE TABLES 


1111 NUM = LDE-LDO+1 

NMOVL = LESTA-LFO+1 

CALL MOVE( 1,X1(LF0),X1 (LDO) , NMOVE , 1 ) 


LFO 

LESTA = 

LO 

LFE 

LDO 

LDE 

RESTBL= 


LDO 

LESTA-NUM 

LO-NUM 

LO-1 

1 

0 

.FALSE. 


RELOCATE FLOW ADJUSTMENT AND STATION TABLES 


1 NUM 


= 3*( NI-N1+1 )+6 
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MAXI = LESIA+NUM 

I F ( ( MAX T -L HU ) . GT .MAXLE ) GO TO 1000 
L F (JNk W = LFO+NUM 
NMOVf = LESTA-LFU+1 

CALL MOVE ( WXKLFO ) , X l ( L FONE W ) , -NMOVE , 1 ) 

LD = LDE + 1 
I F ( LDE ) 2,2,5 
2 LOO = LO 
LD = LOO 
LDE = LDO+NUM-1 
GO TO 6 

5 LDE = LDE+NUM 

6 LFO = LFONEW 

LESTA = M AX T 
LO = LO+NUM 

LFE = LO- 1 
LBLNX T ( LD ) = LD+NUM 

C DETERMINE SWREF FOR BOUNDARY BONAME 

LB = LBF(BDNAME) 

I F ( LB.NE.O ) GO TO 15 
I F ( .NOT. LOWER ) CALL ERR0R1 
SWREF (LD)= SWSAVE 
SWS AVE= 0. 

GO TO 20 

15 SWREF ( LD ) = 0. 

I F ( LfcDEX(LB) . EQ. 0 ) GO TO 20 
IV1 = i 

IV2 = (LED£X(LB)-LBZl(LB))/3-*-l 

SWREF ( L D ) = BARCS< BDNAME, IVl,IV2) 

SWS AVE = SWREF ( LD ) 

20 BNAME l LD ) = BDNAME 
S I GN ( LD ) = -1. 

I F ( LOWER ) SIGN(LD)=i. 

NS E P l L D ) = 0 

I F ( NSLOC.NE.O ) NSEP ( LD ) =LD+3* I NSLOC-Nl+l ) -3 

C MOVE BL PARAMETERS TO TABLE 

30 DO AO LD 1 =N 1 , N I 
SW ( LD ) = SWBL(LDl) 

DSTAR(LD)= DSTR(LDl) 

DDSTAR( LD)=DDSTR( LD1 ) 

LD = LD + 3 
40 CONTINUE 
GO TO 2000 
1000 LUP = UPPER 

I F ( LOWER ) LUP=LOWR 
WRITE (6,1001) LUP , BDNAME 

IOjI FORMAT ( //2X»48H TABLE SPACE EXHAUSTED— BOUNDARY LAYER DATA FOR , 
* A6 , 2 X , 8HB0UNDAR Y ,2X,A6,2X, 9HN0T SAVED//) 

DO 999 LL =1,58,3 
I F ( IBLB(LL) .EQ.O ) GO TO 2000 
I F ( I BLB( LL ) .EQ. BDNAME ) IBLB(LL + 1)=0 
999 CONTINUE 
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2000 IM PDUM( lb ) . EU.O. ) GO TO 2001 
I1TAB = L HO 

CALL TARPRT ( 6HSTABLT ,X1,LESTA,6) 
CALL l A R P R T ( 3HBLB»BLB*60t 10) 

2001 ENTRY i = .FALSE. 

RE TURN 

END 


no ooooo non 


♦ DECK LFSTSQ 

SUBROUTINE L L ST SU ( X , Y , I A , 18 , NOC ,NS * DY) 

♦ LEST SO 1ST/2ND ORDER CURV FIT BY LEAST SQUARE DEV -LESTSQ- 

♦ VERSION 2 

♦ NO ROTATION OF AXIS 

♦ AUTOMATIC REDUCTION OF NS AND NO NEAR THE END PTS 

OIMtNSION X( 10) ,Y< 10).DYI 10) 


INPUT- 

X(I)»Y(I)»I=IA,IB ARE ENTRY COORDINATES 
NOC = ORDER OF CURVE FIT + 1, =2 OR 3 

NS = NUMBER OF POINTS INCLUDED IS EACH LEAST SQUARE FIT 
MINIMUM NS IS =N0+N0-1. ALSO* NS MUST BE ODD. 


OUTPUT- 

OY ( I > = DEVIATION OBTAINED FROM THE CURVE FIT 
COMMON /ERASE / B13),A13,3) 


MIS 

= (NS-D/2 

I AA 

= IA+1 

IBB 

= IB-1 


DY ( I A ) = 0 . 

I F ( IAA.GT.IBB) GO TO 160 
DO IS J 1= I AA, IBB 

C INITIALIZE TO ZERO 

DO 1 L 0 J= 1 » 12 
110 B ( J ) = 0 . 

C SET UP MATRIX < A)( X)=(B) 

A( It 1 ) =NS 

MI = MINO( I-IA,MINO(MIS,IB-I )) 

NO = M I NO ( NOC * M I + 1 ) 

JA = I-MI 

JB = I+MI 

DO 12 U J=JA, JB 
XP=X( J ) -X ( I ) 

YP = Y ( J ) -Y ( I ) 

XP2=XP**2 
At 1,2 )=A( l » 2 ) +XP 
B ( 1 ) = B ( 1 )+YP 
AI2,2)=A(2,2)+XP2 
B(2)=B(2)+YP*XP 
I F ( NO-2 ) 115,120,115 
115 A(2,3)=A(2,3)+XP2*XP 
At 3, 3 > = A ( 3, 3)+XP2**2 
Bt 3 ) = Bl 3)+YP*XP2 
120 CONTINUE 

A ( 2 , 1 )=At 1,2) 

I F < NO-2 ) 125,130, 125 
125 At 1,3 )=A( 2,2) 

At 3,1 )=A( l, 3) 

129 At 3,2 ) = A ( 2, 3) 

130 CALL S I MEQ ( NO , A , B , 3 ) 

DY ( I ) =8 ( 1 ) 
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l'ju continue 

16o DV ( Iti ) = 0 


RETURN 

END 
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♦DECK SIMEQ 

SUBROUTINE S I MEQ ( NN , A , B , MP ) 

CSIMEQ PRO NO F3494A 

THE EQUATIONS WHICH ARE SOLVED ARE AX=B. THE MATRIX 
SIMEQ SIMULTANEOUS EQUATIONS 

A AND THE VECTOR B ARE DESTROYED. FOR PRINTOUT OF 
THE MATRIX TO BE SOLVED SET MP NOT EQUAL TO ZERO 
NN IS THE NUMBER OF EQUATIONS 
DIMENSION A ( 3 » 3 ) * B ( 3 ) 

2b DO l'.u K = l , NN 
30 P= A ( K , K ) 

35 ASSIGN 85 TO MT 
40 DO 55 I =K , NN 

45 IF ( ABS( P )-ABS ( A( I ,K) ) ) 50,55,55 
50 P = A ( I ,K ) 

52 ASSIGN 65 TO MT 

53 l = I 

55 CONTINUE 

60 GO TO MT , < 65, 85 ) 

65 DO 80 J =K , NN 
70 P=A(K,J) 

75 A ( K , J ) = A ( L , J ) 

80 A( L, J )=P 

81 P= B ( K ) 

82 B ( K ) = B ( L I 

83 B ( L ) = P 

85 B(K) = B(K)/A(K,K ) 

IF(K-NN) 90,145,90 
90 L=K+i 

DO 100 J=l,NN 
L JO A ( K , J )=A(K,J)/A(K,K) 

DO 140 I =L , NN 
I F ( A ( I , K ) ) 120,140,120 
120 DO 125 J = L , NN 
125 A( I , J ) = A( I , J ) - A ( I , K ) * A { K , J ) 

140 B( I ) = B( I ) — A ( I ,K )*B (K ) 

145 L=NN- 1 

DO 170 KK= 1 , L 

K=NN-KK 

P=0.0 

DO 165 J = K» L 
165 P = P + A IK, J + l )*Bl J + i ) 

170 B ( K ) = B ( K ) -P 
1999 RETURN 
END 
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oooooonooooooooon oooonnooon 


♦ DECK ERRO-. m 

SJBKUUl INt ERROR! 

CEDUMPN STC EDUMP - INPUT LINK -EDUMPN* 


C 

C 


COMMON /ALLCOM/ M ACHA , PS A , T S A , P TA , TT A , 
i MACHC,PSC,TSC,PTC,TTC, 

DAXI T,SCALEA,TTE,CHOTST 
REAL M ACHA ( 1 ) , MAC HC 

LOGICAL axia,axic 

LOo I CAL CHOTST 

CHANNtL INPUT DATA TABLE 


AX I A , RGA f GAMA 
AXICfRGCt GAMC 


♦ 

t 


I 

<L 

A 


INDEX- LH = L HO , L HE 

COMMON /CHDATA/ CHNAM (1) , LHNE XT < 1 ) , WTF LO W ( 1 ) , TTO (1) , P TO 1 1 ) 

TSO( 1) tPSO(l) , MACHO (1) , AO ( 1 ) , VAR Y ( 1 ) , 

RG ( 1 ) * GA M ( 1 ) . N R ( 1 ) i NC ( I ) ,TAB(6) , 

BB ( 75 ) 

LOGICAL VARY 

INTEGER CHNAM 


» 


DIMENSION VOID 

REAL MACHO 

EQUIVALENCE (VO, MACHO) 

BOUNDARY TABLE 
INDEX- LB=L8DQ,LBDE 

L BNEX T = INCREMENT TO NEXT BOUNDARY 

LBZI = INCREMENT TO THE FIRST BOUNDARY POINT ( = 0 BEFORE COALL AT 1 0 
CHNAME= CHANNEL WITH WHICH THE BOUNDARY DATA IS ASSOCIATED 
UP = T OR F FOR UPPER OR LOWER BOUNDARY 

LEDEA = RELATIVE INDEX OF L.E. POINT WHEN LOWER AND UPPER SURFACE 
CONTOURS ARE CONNECTED 

BDNAME, LBA,LBB=NAME AND INDEX LIMITS OF SPECIFIC BOUNDARY 

OATA WHEN BOUNDARIES ARE COALLATED 
DIMENSION BDT( I) »LBNEXT(i),LBZi(l)» 

. CHNAME( I ) ,UPl I) tLEDEX(I) , 

^ ZBT( I ) ,RBT(I) ,ANGBTU2) 


LUGICAL UP 

INTEGER BDT , CHNAME, BDNAME 

DIMENSION BDNAME ( I ) ,LBA(1 ),LBB(l) 

EQUIVALENCE ( BDNAME , ZBT ) , ( LB A ,RBT ) , (LBB.ANGBT) 

FLOW ADJUSTMENT TABLE 
INDEX- LF=LFO,LFE 
NFCOLS= B 

X 1 F = ORTHOGONAL COORDINATE 

X 2 F = STREAMLINE COORDINATE OF SL EMINATING FROM T.E. 

X1BF = X I-COORD I NATE OF CHOKE STATION OF FLOW BELOW T.E. 
XIAF = X 1-COORDINATE OF CHOKE STATION OF FLOW ABOVE T.E. 

SIF = SI-COORDINATE OF T.E. (UPPER SURFACE). THIS ITEM 
IS USED WHEN INTERPOLATING FOR WAKE DELTA-STAR. 
LFB,LFA=INDICES OF STATIONS BELOW AND ABOVE T.E. 
NCHB,NCHA=NUMBER OF CHANNELS BELOW AND ABOVE T.E. 

LRF = INDEX OF DUMMY ORTCHN LIST FOR THE T.E. 

LRXF = INDEX OF LAST CHANNEL BELOW THE T.E. 

JORDER= 0 IF TOTAL FLOW AT XiF IS GIVEN 
= 2 IF FLOW ABOVE T.E. IS GIVEN 
= I IF FLOW BELOW T.E. IS GIVEN 
J OR DE R = -I IF FLOW AT XIF IS CHOKED AND SINGLE CHANNEL 
DIMENSION X1F( I ) ,X2F ( 1 ) ,XIBF ( I) , XIAF (I ) , 

i S1F(1) , NCHB ( i ) f NCHAl l) , JORDER(I) ,VNR< 12) 

EQUIVALENCE ( L FB , X IB F ) , ( LFA , X I AF ) , ( L RF , NCHB) , (LRXF.NCHA) 
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DIME NS ION L E B ( l),LFA(l)»LRF(l) » LRXF { 1 ) 

c table nr convicted properties 

C INDEX- LT=LTO,LTE 

C CH CHANNELNAME 

C L TNt X T - INDEX INCREMENT TO THE NEXT CHANNEL 

C LPSI = RELATIVE LOCATION OF PSI LIST 

C NpT = NO. OF PSI, TT, PT ANO RCU VALUES 

C LTT - RELATIVE LOCATION OF TT LIST 

C LPT = RELATIVE LOCATION OF PT LIST 

C LRCU = RELATIVE LOCATION OF RCU LIST 

DIMENSION CH( 1 ) ,LTNEXTU I ,NPT(l) .LPSIU I ,LTT(1 ) ,LPT( 1 I , 

i LRCU ( 1 ) » 

t CRGI II ,CPGJ( i),C2CP( 1) ,QGAM(1) ,FGT( 1) ,FGPI 1) , 

3 FGR(l) ,AREAT8(485) 

INTEGER CH 
DIMENSION XCH( L I 
EQUIVALENCE (CH,XCH) 

C STREAMLINE TABLE 

COMMON /SLTAB / W ( 1 2 8 ) , X 2 1 1 2 8 I , SLC HN ( 12 8 ) 

INTEGER SLCHN 
STATION TABLE 
INDEX- L=LO,LESTA 

S C HUK t = STATION CHOKE INDICATOR ( ADJWF ,BRHS , WRIOUT I 
MCL = SHARP CORNER INDICATOR IBLDTBSI 

MCL = FIELD INDEX OF CONTROL STREAMLINE ( PT MOV E , FLOBAL ) 

DIMENSION X I ( 1 I , LNEXT ( 1 ) ? MLB ( 1 ) ,MUB( 1 I ,PRIM<1 I , 

1 TYPE LB 1 1 I ,NAMELB ( 1 ) . ILB (II , FLB (1 ) , S 1LB (II , 

i TYPEUBI 1 I .NAMEUBI 1 ) , IUB(l) ,FUB(1 ) ,SIUB{ 1) , 

3 VMBl II »DWDV( 1 ) » X2CLU) , VCL ( 1 ) ,MCL ( 48 1 > 

LOGICAL PRIM 

INTEGER TYPELB, TYPEUB 
DIMENSION SCHOKE(l) 

EQUIVALENCE ( SCHOKE , DWDV I 

TABLE OF WAKE DISPLACEMENT THICKNESS 
INDEX- LW=LWO,LWE 

DIMENSION X2W(1> , L WNE XT ( i I ,SIW(47) 

di mens ion dsti i i 

EQUIVALENCE (DSTtSlW) 

SUBTABLE ARRANGEMENT IS- 

X2W,LWNEXT(=2 + 2N I , S1WI1),S1W(2I...S1W(N), DSTI 1 ) ,DST( 2 ) . . .DST ( N I 
X2W = STREAMLINE COORDINATE 

S1W = DISTANCE ALONG STREAMLINE FROM T.E. 

DST = WAKE DISPLACEMENT THICKNESS AS A FUNCTION OF S1W 

C FIELD TABLES 

C INDEX- M=MO,NM 

COMMON /CZ / ZI300I 
COMMON /CR / R < 300 I 
COMMON /CS2 / S2 ( 300 I 
COMMON /CS1 / S l ( 300 I 
COMMON /CPHI1 / P H 1 1 ( 300 I 
COMMON /CM / JMSI 300) 

COMMON /CCURV / CURV(300) 

COMMON /CB / B ( 300 I 

COMMON /CIDEX / M , J , MU ,MD , I S TAG 

C TABLE UF INDEX LIMITS 

COMMON /IXQRIG/ LHO, LHE, L8D0,LBDE , LTD ,L TE , LWO,LWE, LFO»LFE , 
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* 

DIMENS IUN 
TABLt OF LEADING 


L 0 ,LE STA , LSO » LSE *LDO,LDfc , LDUM (4) , 

NM , NJ.NFCQLS, MAXN J , MAXOL , MAXNM , MAXLE t 
LEO, LEE. LRO , LRE , LRD 
L I M I T$<24 ) 

EDGE AND TRAILING EDGE POINTS 


INDEX- L E =L EO , L EE , 10 

NLE»NTE=NO. OF L.E. AND T.E. COINCIDENT PTS, RESPECTIVELY 
CHL,CHU=NAME OF CHANNEL ABOVE ANO BELOW PT, RESPECTIVELY 
BDL,BDU=BOUNDARY NAMES ASSOCIATED WITH THE POINTS 

NUSED = COUNT OF TIMES THAT POINT USED IN CONSTRUCTION OF /ORTCHN/ 
COMMON /LETEPT/ XE ( I ) , YE (1) , ANGE ( 1 ) , NLE ( l ) ♦ NT E { I ) , 

I C HL ( 1 ) ,CHU( 1 1 ,BDL ( I ) , BDU( 1 ) » NUSED (49 1 ) 

INTEGER CHL,CHU,BDL,BDU 

TABLE OF CHANNELS EMBRACED BY EACH ORTHOGONAL 
INDEX- LR=LRO*LR£ ,LRD 

LRD = NUMBER OF CHANNELS PLUS ONE, LR INDEX INCREMENT 
LEDGE = INDEX OF THE ORTHOGONAL POINT IN THE LETEPT-TABLE 
LRPRL V= POINTER OF LINE OF UPSTREAM CHANNELS IN ORTCHN-TABLE 
CHN A = CHANNEL NAMES 

COMMON /ORTCHN/ LE OGE < l ) , LRPREV (1) , CHNA ( 479 ) 

INTEGER CHN A 


DIMENSION JCHNA(l) 

EQUIVALENCE ( JCHNA ,CHNA ) 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

equivalence 


( CHNAM , 8DT ,CH»X2W»X1F ,Xl) 

<LHNEXT,LBNEXT,LTNEXT,LWNEXT,X2F,LNEXT ) 

( WTFL0W,LBZ1,NPT,S1W,X1BF,MLB> 

( TTO,CHNAME,LPSI ,XIAF,MUB) , I PTO, UP , LT T , SI F , PR I M > 
(TSO,LEDEX,LPT,NCHB,TYPELB) 

( PSO, ZBT, LRCUfNCHA, NAMELB) 

(MACHO, RBT,CRG,JORDER,ILB) , ( AO, ANGBT ,CPGJ , VNR, F LB ) 
(VARY,C2CP,SILB) , I RG ,QGAM, TYPEUB ) 

( GAM, FGT, NAME UB ) , ( NR ,FGP , I UB ) , ( NC »FGR*FUB) 

( T A B ( 1)*AREATB,S1UB), ( TAB{ 2 ) , VMB ) , ( T AB ( 3 ) » DWDV ) 

( TAB(4) ,X2CL) , ( TAB ( 5 ) , VCL ) , (TAB(6),MCL) 


COMMON /CBITS / 
COMMON /CREDIN/ 
EQUIVALENCE 
I 

COMMON /CTABPR/ 


BITS, BLANK 

Z TRANS, R TRANS, ROT ATE ,Z PIVOT, RPIV0T,SCAL E,NB,TBB(9 ) 
( XTRANS,Z TRANS) , I YTR ANS , RTR ANS ) , ( XP l VOT , Z P I VOT ) , 

( YPIV0T,RPIV0T) 

I1TAB 


i 


COMMON /CSEGME/ 


COMMON 

COMMON 

COMMON 


/CSMOOB/ 
/CSMOOC/ 
/BLBDY / 


I A ( 1 0 ) » I B ( 
N ,NSEG, N I 
XA(IOO), YA 
DUM1 ( 200 ) , 
IBLB (60) 


DATA TXA/2HXA/, TZA/2HZA/ 


10) ,IMA( 10) ,IMB( 10 
, N I M 

(100) , DEVI (100) 
ANG(IOO) ,DUM2(400) 


) , JT YPE (10) , 


,DEV( 100) ,CURVB( 1 00) 


I GGO = 1 
GO TO 1777 
ENTRY EDUMP 
I GGO = 2 

1777 CONTINUE 

1100 FORMAT (/// 1X36HCHANNEL INPUT DATA TABLE, /CHDATA/ ~) 
WRITE (6,1100) 

I1TAB = LHO 
NC X = NC 
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IHNCX.LT. 3) NCX = 5 

CALL TAHPRT( BLANK tCHNAM, LHEtNCX ) 

112U FORMAT ( ///IX 5 4HB0 UNO AR Y COORO INAT6 S AND ANGLES IN RADIANS* /BDYTAB 
*/ - ) 

WRITt (6,1120) 

UTAH = LBDU 

CALL I ABPRT (BLANK, BOT ,LBDE,3) 

1 1 1 j FORMA f ( / / / 1 X 4 1HTABLE OF C QNVECTED PROPERTIES, /CONVTB/ -) 

WRITt- (6,1110) 

UTAH = L TO 

CALL T ABPRT ( BLANK , CH ,LTE ♦ 7) 

IF(LEE.LT.LEO) GO TO 140 

1130 FORMAT ( / / / 1 X 1 25H0R DERED LIST OF UPSTREAM BOUNDARY PNTSt L.E. PNTS, 

♦ T.t. PNTS, ANO DOWNSTREAM PNTS WITH REFERENCES TO CHANNELS AND BO 

♦UNDAR IES./1X10H/LETEPT/ - //4 X2HLE 6 X , 2H XE 1 OX , 1 5HY E ANGE12X, 

♦3HNLE9X, 12HNTE CHL9 X , 3HCHU9X , 3HBDL9X , 3HBDU1 OX , 5HNUSED ) 

WRITE (6,1130) 

I1TAB = LEO 

CALL T ABPRT ( BLANK , XE , LEE , 10) 

140 IF ( LRfc .LT .LRO ) GO TO 150 

1140 FORMAT ( /// 1X98HTA8UL AT ION OF CHANNELS EMBRACED BY THE ORTHOGONALS 
♦WHICH PASS THROUGH THE ABOVE POINTS, /ORTCHN/ -//4X26HLR 
♦LE LR-PREV) 

WRITE (6,1140) 

UTAB = LRO 

CALL I ABPRT (BLANK,LEDGE»LR£» LRD ) 

1150 FORMAT! / / / 1 X . 7H S TRE AMLINE TABLE-/17X32HJ X2 SLCHN 

♦ w/< U8,F12.6,6X,A6,F12.6, ), ) 

150 WRITE (6,1150) ( J , X2 ( J ), SLCHN (J ), W ( J ), J=1 , NJ ) 

1190 FORMAT! ///1X37HWAKE DISPLACEMENT THICKNESS, / WAKETB// / 1 IX 19HX2W/ S 1 
*W DST ) 

WRITE (6,1190) 

IiTAB = L WO 

CALL TABPRT(BLANK,X2W,LWE,2) 

UBO F0RMAT(///1X43HTABLE OF FLOW ADJUSTMENT STATIONS, /CADJWF/ //15X3HX 
*1F9X, 3HX2F8X,4HX1BF8X,4HX1AF9X,3HS1F8X,4HNCHB8X, 16HNCHA JORDE 

♦ R ) 

WRITE (6,1180) 

IITAB = LFO 

CALL T ABPRT ( BL ANK , XI F , LF E ,NFCOL S ) 

1160 FORMAT ( ///1X25HSTATI ON TABLE, /STATAB/ -) 

WRITE (6,1160) 

UTAH = LO 

CALL T ABPRT ( BLANK,Xi,LESTA,5 ) 

CALL JMSPRT 

1170 FORMAT ( ///1X19H FIELD COORDINATES -) 

WRITE (6,1170) 

CALL TABPRT ( 1HZ ,Z ,NM, 10) 
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CALL TABPRT ( 1HR ,R , NM , LO) 


PRINT UVLRALL DATA 
CALL f ABPRT ( 6HALL COM , MAC HA, 2 0, 8 ) 

in I B L B ( 1J.NL.0 ) CALL TABPR T < 5HBLB0Y, I BLB , 60 , 3 ) 
I F ( LDE.EQ.O > GU TO 1321 
UTAH = LDO 

CALL T ABPRT ( 5HBL TAB, CHNAMtLDE » 3 ) 

1321 CONTINUE 

1 F ( IGG0.EQ.2 ) RETURN 
LSTOP = 5 

GO TO ( 900 . 1777 ) , LSTOP 
90u RETURN 
END 
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-REDINP- 


* DE CK RED IMP 

SUBROUTINE REDINP 
♦REDINP STC READ INPUT 

COMMON /BCOMMN/ PROGM ( 8 ) , PROGSV , F I L I N , F I LOT 
LOGICAL FILIN, FILOT 

COMMON /ADAMOL/ N AME I 6 ) , ADDRb S( 6 ) , TI TLE C 6 ) , I DENT ( 6 ) 

COMMON / ADAMO 2/ ENDJOB ,NUMPLT ,PLOTED ,ENDCRD 
LOGICAL ENDJOB, PLOTED ,ENDCRD 

COMMON /ALLCOM/ MACHA, PSA ,TSA ,PTA, TTA, AXI A,RGA, GAMA, 

A MACHC,PSC,TSC,PTC,TTC, AXIC,RGC,GAMC, 

2 DAXIT ,SCALEA,TTE,CHOTST 

RfcAL MACH A ( I ) , MACHC 

LOGICAL AXIA.AXIC 

LOGICAL CHOTST 

LOGICAL AXI 

REAL MACHO l I ) 

EQUIVALENCE ( MAC HO, MACHC) , ( PSO , P SC » , ( TSO , TSC ) , (PTO, PTC) , 

x ( T TO , T TC ) , (AXI , AX I C ) , (RG»RGC) • ( GAM , GAMC ) 

COMMON /CSS / SSFML,SSEF, SSEANG, SSDF ,SSFEND,SSFNDI 
1 , S SOLE, AAF ACT, BRLX,CURRLX,T SIC, RHOC,RHOCSS 

INTEGER SSFML 

LOGICAL SSEF, SSDF, SSDLE 

SSFML = SUPERSONIC CURVATURE FORMULA NUMBER 
SSEF = SUPERSONIC ENTERING FLOW, T OR F 
SSEANG= ENTERING FLOW ANGLE (DEGREES) FOR SSEF=T 
SSDF = SUPERSONIC DISCHARGE FLOW, T OR F 

SSFEND= SUPERSONIC BEAM DOWNSTREAM END CONDITION, =0,1 FOR PARABOL 
SSFNO i= SUPERSONIC BEAM UPSTREAM END CONDITION, =0,1, FOR PARABOLA 
SSDLE = SS FLOW BELOW AND AFT OF LE PT, T OR F 
AAF AC T= CENTRAL POINT INFLUENCE COEFFICIENT FACTOR 
BRL X = B-REL AX AT I ON FACTOR 
CURRL X = CURVATURE RELAXATION FACTOR 

RHOC = ACCELERATION FACTOR ON CURVATURE ITERATION AT 
SUBSONIC POINTS 

RHUCS S= ACCELERATION FACTOR ON CURVATURE AT SUPERSONIC POINTS 
COMMON /IXQRIG/ LHO, LHE , LBDQ *LBDE , LTO*LTE , LWO,LWE, LFO,LFE , 

* LO.LESTA, LDUMA ( 2 ) »LDO»LDE , LDUMU), 

« MO,NM, N J ,NFCOLS, MAXNJ ,MAXOL , MAX NM, M AXLE , 

* LEO, LEE, LRO, LRE »LRD 

DIMENSION LIMITS124) 

EQUIVALENCE (LIMITS, LHO) 

COMMON /SLTAB / W ( 128 ) , X2 ( 128 ) , SLCHN ( 128 ) 

INTEGER SLCHN 

COMMON /BENDIN/ NBC IN ( 2 ) , ACF ( 2 ) 

COMMON /CB / B ( 300 ) 

COMMON /CB I TS / BITS, BLANK 

COMMON /CCRX / CR XSL , CR XOL , CRXSS ,CRXE ,CRXC ,CRMACH 
DIMENSION CRX ( 6 ) 

EQUIVALENCE (CRX,CRXSL> 

COMMON /CHDATA/ T ABLES ( 2046 ) 

COMMON /CEND / T BLEND ( 2 ) 

COMMON /CIADIN/ RHOBAS , RHOAMP , I ADM 
COMMON /CINNER/ I NRC TR ,RDUM, N INNER ( 1 6) ,CNVF ( 16 ) 

COMMON /CISBOT/ F ARFLD ( 2 ) , FREE ( 2 ) , PRES < 2 ) ♦ RFF , N2P, 

1 ZP(10),PPS(10>, A1 , A2» ADUM(6) 

INTEGER FARFLD, FREE, PRES 

COMMON /CLINES/ L I NE S ,0M I TFK , PT I TLE ( 6) 
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COMMON / C L W 0 S V / 
COMMON /CM / 

COMMON /CM Ax I T/ 
LOGICAL 

COMMON / CNORM / 
C OMMOt'l / C N I K L / 
COMMON /CPRINT/ 
COMMON /CPRPRN/ 
INTfc GER 

COMMON / C PT MOV / 
LOGICAL 

COMMON /CR / 

COMMON /CREFIN/ 


COMMON 

COMMUN 

COMMON 


/CS1 / 
/CS2 / 
/CT APOS/ 


LOGICAL 


COMMON 

COMMON 


/CTHICK/ 

/CTOLRL/ 


I 

* 


COMMON 
COMMON 
COMMON 
COMMON 
LOGICAL 
COMMON /TAPES 


/CVM / 
/CZ / 
/SPACER/ 
/TROUBL/ 


L WOS V 
J M S ( 300) 

MAXI T » MA JCTR #GREF I N»EDUM 
GREFIN 

R HL ♦ RM t AHL , ARM 
Kill ) ,STA(6) ♦ INSERT 

PRTES2 tPRTB, PRTA,PREFIN,PREFN2,SS0NIC,PDUM( 20) 

PRPRN 

PRPRN 

VELPQT, ICOB, NODENS*FBASTG 
VELPQT 
R F ( 300 ) 

SLS, SG2itVMGi»VMG2» NGRtNGZ » SGR( 10) ,GR<10> , 

SGZ ( 1 0 ) * G Z ( 1 J ) 

S i ( 300 ) 

S2 ( 300 ) 

RESTRT ,EN0BDT,STCFIL f K6SV 
RESTRT »ENDBDT »STCF IL 

NTHKX, NTHKY, THKX(20) ♦THKYI20) ,THIK2D(78) 

TOLRL»MAXSWP,CLEN,DS2MX,TOLES2,NSWP, 

DS1DMP,DS1MXA,DS1MXB,0S1RMS,ES2MX 

0S1RM0,SG1MIN,T0LINR 

VMF ( 300) 

Z F ( 300 ) 

MAXLH.MAXLT, MAXLF ,MAXLW 
ERR,ERRMAJ # INERR, PRERR 
ERR » ERRMAJ, INERR » PRERR 
NTAPO»NT APN 


COMMON / BL BOV / BLBI60) 
DIMENSION I BL B ( 60 ) 

EQUIVALENCE I I BLB » BLB ) 

COMMON /VISCOS/ TREF , MURE F , SCON 
REAL MUREF 

COMMON /REBL / RESTBL 
LOGICAL RESTBL 

COMMON /CGRAV/ CO 


LOGICAL 


F IRST 


DATA KA/1HA/, KBDY/3HBDY/, 
DATA FIRST/ .TRUE./ 


KCHN/3HCHN/, KSTA/3HS T A/ 


ENDCRD= T IF END OF CARD INPUT 
ENDBDT = T IF END OF BOUNDARY DATA 
STCF l L= T IF A STC-SUBFILE EXISTS 


TAPE 

TAPE=ORGF. 
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* PLOT, I PLOT, SAME XY, XSC ALE, YSC ALE, 

A R HO BA S , RHOA MP , I ADM , 

B NTHKX,NTHKY, THK X,THKY»THIK2D» 

C LBL,MAXL8L, T0LL8L 

* ,TREF,MUREF,SCON , CG 


C** INITIALIZE AND READ OVERALL (A) INPUT DATA 

I F t .NOT. FIRST .AND. (K5.NE.KA .OR. ENDCRD) ) GO TO 200 
I F ( FIRST .AND. K5.EQ.KA) GO TO 100 
WRITt (6,1000) 

ERR = .TRUE. 

PROGS V= 0. 

GO TO 200 


100 PROGS V = 0. 

END8D T= .FALSE. 

FIRST = .FALSE. 

LINES = 64 
MA J C I K - 0 
R E S TR T - .TRUE. 

STCFIL= .FALSE. 

CALL SETMI 1 , B I TS » MACHO, 8) 

C DETERMINE FIELD ARRAY SIZE 

MAXLE = L0C2( TABLES, TBLEND) 
MAXNM = L OC 2 ( RF , ZF ) 

GO TO 120 


C READ CARD INPUT 
130 READ (5, A) 

DO 135 1 = 1,8 

135 IF( MACHO( I ) .NE.BITS) MAC HA ( I I =MACHO( I ) 
DATA AOOOOO/6HAOOOOO/ 

K6S V = AOOOOO 


C DEFINE THE CHARACTERISTIC LENGTH, CLEN 
CLEN = SGR 
IF(NGR.LE.i) GO TO 146 
DO 144 1=2, NGR 
144 CLEN = CL EN + SGR ( I ) 

146 IF(NGZ.LE.O) GO TO 149 
DO 146 1=1, NGZ 
146 CLEN = CL EN+ SGZ < I ) 

149 CLEN = CLEN/FLOAT (NGR+NGZ) 


C 


>ET UP INDEX-ORIGIN TABLE 
ORDER OF TABLES IN BLOCK 

LH /CHDATA/ 

LB /BDYTAB/ 

LT /CONVTB/ 

LW /WAKETB/ 

LF /CADJWF/ 

L /STATAB/ 

IF(STCFIL) RETURN 
RES TRT= .FALSE. 

LBDO = LHO+MAXLH 
LBDE = LBDO 
(OTHER INDEX LIMITS 


IF THERE 
COMMON 


IS NO STC-TAPE INPUT 


ARE SET IN SUBROUTINE BLDTBS ) 
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RETURN 


C Rt AD INPUT FILE 

120 IF C .NUT .f IL IN ) GO TO 130 
REWIND NT API) 

READ (NTAPO) S TC F IL , ( L I M I TS ( I) , I = l , 24 ) 

LWOSV = L WO 
IF(STCFIL) GO TO 125 
E NDBD T = .TRUE. 

WRITE (6,1120) 

GO TO 130 

125 RE AO (NTAPO) ( ( I DENT < I ) , I =1 , 6 ) , AX I ,RG , GAM , MACHO, PSO , T SO , PT0, TTO , 

1 PRPKN,TTE,CHOTST,MAXIT ,MAJCTR, (NINNERII ) , I = 1, 16) , VELPOT , I COB, 

2 NUDENS,RN,NGR,NGZ, (SGR ( I ), 1=1,40) , VMG1 » VMG2 , INRCTR, SLS,SG21, 

3 NBC IN ( 1 ) , NBC I N ( 2 ) , ACF ( 1),ACFI2) , SSF ML , SSEF , SSE ANG , SSDF , S SF END, 

4 SSFND1 , SSDLE, A4F AC T , BR L X , C URR LX , TS I C , (FARFLDI I ), 1 = 1,8 ) , 

* RHOC , RHOCSS »RHL »RM, 

* TREF, MUREF, SCON, (BLB ( I ) , 1=1 ,60 ) , 

5 (2P( I ), 1=1,28), (TABLES! I ) , I=1#LESTA) , ( B ( I ) , I = 1 , NM ) , (JMS(I), 

6 I = . » NM ) , (SKI), 1 = 1, NM), ( S2(I ) , 1 = 1, NM) , C ZF (I ) , 1 = 1, NM) , ( RF (I) , 

7 1 = 1, NM), ( VMF ( I ) * I = 1 *NM ) , ( W< I ) , I = 1 , NJ ) , 1X2(1 ) « I = 1 ,NJ ) , 

8 ( SLCHN ( I ), I=1,NJ),T0LRL,MAXSWP,T0LES2,T0LINR,SG1MIN,DS1DMP, 

A DS1RM0, (CRX( I ) , 1 = 1,6), RHOBAS , RHOAMP , I ADM , NTHKX , NTHKY , 

B (THKX( I ), 1=1, 118) ) 

C CHECK TO SEE IF STC-A INPUT DATA EXCEEDED DIMENSIONS 

IF (NM.GT.L0C2(RF,ZF) .OR. LE STA .GT .L0C2 1 TABLE S ,T BLEND ) ) ERR=.TRUE . 
I F ( LDE.NE.O ) RESTBL=.TRUE. 

GO TO 130 

C READ BOUNDARY DATA 

200 CALL RBD 

IF(ENOCRD) GO TO 700 
IF ( K5.E0.KBDY ) RETURN 

C READ CHANNEL DATA 

300 I F ( K5 *NE . KCHN ) GO TO 400 

C IF RESTRT, UNPACK TABLES TO MAKE ROOM FOR NEW CHDATA AND CONVTB. 

IF ( .NOT.RESTRT .OR. LBDO.GT . (LHE+ll ) GO TO 350 
MOVE 1 = L0C2(TABLES,S1)-LESTA 
MOVE c = MOVE 1/2 
L W TO = LW0+M0VE1 
LBTO = L8D0+M0VE2 

CALL MOVE ( 2 , T ABLES ( L WO ), TABLE S ( LWTO) , LWO-LE STA-1 , 1 , 

1 TABLES (LBDO), TABLES (LBTO) ,LBD0-LTE-1 , 1 ) 

LBUO = LBDO+MOVt 2 
LTE = LTE+M0VE2 
LBDE = LBDE+M0VE2 
LTO = LT0+M0VE2 
L WO = L WO+MOVE 1 
350 CALL RCD 
RETURN 

400 WRITE (6,1690) 

E RRMA J = .TRUE. 

RETURN 

C CONSTRUCT LETEPT, OrTCHN, CONVTB, SLTAB, STATAB AND THE FIELD TABLE 
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7 JO I F ( EKKMAJ .OK. LBOE.EQ.LBOO) ERR*. TRUE. 

900 RETURN 

1000 FORMAT ( /1X73HERRUR- THE K5=A INPUT DATA DOES NOT IMMEDIATELY FOLLO 
*W THh PROGM=STC CARD) 

1120 FORMAT ( //1X43H*** NO STC DATA FOUND ON THE INPUT TAPE.//) 

1690 FORMAT ( // 1X44H** PLEASE CHECK THE INPUT VALUE OF K5 CK5*A6.18H). 
* IT MUST BE 0NE/6X37H0F THE FOLLOWING- A* BDY, CHN, STA.//) 

END 
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♦DECK RBD 

SUBROUTINE RBD 

♦RBD READ IN BOUNDARY DATA -RBD- 


iNPur- 

ENDBDT = END Of BDY/STC TAPE RECORDS, T OR F 

ENDCRD= END OF ALL S TC CARD INPUT, T OR F 

K6SV = VALUE OF K E Y ( 6 ) OF LAST RECORD READ FROM TAPE 

RLSTRI= RESTART (WITH EXISTING TABLES) IS TRUE ONLY 

IF CARD BUY- DAT A HAS NOT YET BEEN ENCOUNTERED 
STCFIL= T IF A STC-SUBFILE EXISTS ON TAPE=QRGF . 

OUTPUT - 
ENDBDT= 

K6SV = 

RE S TR T = 


INTEGtR 

COMMON / BCOMMN/ 
LOGICAL 


RE FS » BDY , CHN 

PROGM ( 8 ) tPROGSV.F ILIN»FILOT , REFS (5 ) 

FILIN, FILOT 


COMMON /ALLCOM/ 


REAL 

LOGICAL 

BOUNDARY TABLE 


MACH A »PSA,TSA,PTA,TTA, AX I A.RGA, GAMA, 
MACHC,PSC ,TSC,PTC,TTC, AX I C ,RGC, GAMC , 
DAXI T,SCALEA,TTE,CHQTST 
KACHA ( 1) , MACHC 
AXIA,AXIC,CH0TST 


INDEX- LB=LBDO, LBOE 

LBNEX T= INCREMENT TO NEXT BOUNDARY 

LBZi = INCREMENT TO THE FIRST BOUNDARY POINT (=0 BEFORE COALLATIO 
CHN AM E = CHANNEL WITH WHICH THE BOUNDARY DATA IS ASSOCIATED 
UP = T OR F FOR UPPER OR LOWER BOUNDARY 

LEDEX = RELATIVE INDEX OF L.E . POINT WHEN LOWER AND UPPER SURFACE 
CONTOURS ARE CONNECTED 

BDNAME, L6A, LBB=NAME AND INDEX LIMITS OF SPECIFIC BOUNDARY 

DATA WHEN BOUNDARIES ARE COALLATED 
COMMON /CHOATA/ BDT( i)»LBNEXT(l),LBZl(l) , 

1 CHN A ME (1),UP(1),LEDEX(1), 

2 ZBT(l) ,RBTU ) ,ANGBTU2) 


LOGICAL 
INTEGER BDT , 
DIMENSION 
EQUIVALENCE 
COMMON / I XOR 1 

♦ 

* 

* 

DIMENSION 

EQUIVALENCE 


UP 

CHNAME, BDNAME 

BDNAME(l) , LB A ( 1 ) , LBB ( I ) 

( BDNAME, ZRT) , ( LBA,RBT ) , ( LBB , ANGBT ) 

G/ LHO, LHE , LBDO »LBDE , LTO, LTE , LWO,LWE, LFO,LFE , 
LO,LESTA , L D UM ( 8 ) , 

MO ,NM , N J » NFCOL S » MAXN J , MAXOL , MA XNM , MAXLE , 

LEO, LEE, LRO ,LRE , LRD 
LIMITS(24) 

(LIMITS, LHO) 


COMMON / ADAM02/ 
LOGICAL 

COMMON /CBITS / 
COMMON /CLINES/ 
LOGICAL 

COMMON /CNTRL / 
EQUIVALENCE 
COMMON /CPI / 
COMMON /CREDIN/ 


endjob,numplt,ploted,endcrd 

ENDJOB, PLOTED,ENDCRD 

BITS, BLANK 

LINES,0MITFK,PTITLE(6) 

OM I TFK 

K 5, BDY (6 ) , INSLRT,CARRY,CHN 
( BDY, IBDY ) 

PI , T WOP I , PIQ2 ,PIQ4,T0DEG,T0RAD 

Z TRANS, R TRANS ,ROT ATE , ZPi VOT , RP I VOT , SC AL E , NB , T AB ( 9 ) 
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fcCU [ VALENCE 


( X TRANS, Z TRANS) , ( Y TRANS » RTRANS ) f ( XP I VOT , Z P I VOT) , 

* l YPI VOT.RPIVOT) 

COMMON /CTAPOS/ R E STR T , E NDBD T , S TCF I L , K6SV 
LOGICAL RESTRT,ENDBDT»STCFIL 

COMMON /ERASE / B(800) 

COMMON /SPACER/ M A XL H , MA XLT , MAX LF , MAXL W 
COMMON / TROUBL / ERR , ERRMA J f I NERR , PRERR 
LOGICAL ERR,ERRMAJ,INERR»PRERR 

C SMOOTH COMMONS 

COMMON /ADAMOl/ NAME ( 6 ) , ADDR ESI 6 ) , T I TIE ( 6 ) , I DENT! 6 ) 

COMMON /CALCPT / DX,XMOD 
COMMON /CELLPT / DZETA 

COMMON /CSEGME/ I A ( 10 ) , I B ( 10 ) , I MA ( 10 ) , I MB ( 10) , JT YPE ( 10 > , N » NSEG , 

1 N I I , N I M 

EQUIVALENCE (NI,NII) 

COMMON /CSMOOA/ DE VA ( 20 ) , FEND A( 20 ) , ANGA ( 20 ) ,CURV A ( 20 ) .NARB 
COMMON /CSMOOB/ XA ( 100 ) , YA( 100) .DEVI ( 100 ) 

DIMENSION ZA( 100) ,RA( 100) 

EQUIVALENCE ( ZA.XA), (RA,YA) 

COMMON /CDS2 / X < 100 ) , Y < 100 ) , ANG < 100 ) , ANGD < 1 00) , CURV < 100 ) . S ( 100 ) , 
I EQKl 100) »DEV( 100) »CURVB(100) 

DIMENSION Z ( 100) ,R ( 100) ,DUM( 100) 

EQUIVALENCE (Z,X),(R,Y), <DUM,CURVB) 

COMMON /BLBDY / BLB(60) 

DIMENSION I BL B ( 60 ) 

EQUIVALENCE ( I BLB ♦ BL B ) 

LOGICAL BL 
DATA LBLB/1/ 

LOGICAL DATAIN,ENDBDC,UPPER,ZRONLY 


DATA K BDY / 3HBDY / , KHIGH/6H / 


NAMELIST /A/ 

B, 

NB, 

TAB, 

DBLPTS, 

ZRONLY 

i 

BDY, 

CHN , 

UPPER, 

X , z , 

Y, R, 

ANGD 


ROTATE, 

ZP I VOT, 

RPI VOT, 

ZTRANS , 

RTRANS, 

SCALE 

3 

FLIP, 

XP I VOT, 

YPI VOT, 

XTRANS, 

YTRANS, 

DUM 


, IDENT, 

DX, 

XMOD , 

DEVA, 

FENDA, 

ANGA 

5 

CURVA, 

Z A , X A , 

RA ,YA, 

DEVI , 

Nil, 

DEV 

6 

ANG, 

CUR V, 

CURVB, 

FOK, 

S, 

NIM 


7 UPPER 
♦ , C APX 1 , BL 

C DEFINTE DOUBLE POINT TOLERANCE, DPTOL 

DPTOL = l.E-5 

C INITIALIZE 

C E N D B D C = END OF BDY CARD INPUT, T OR F 

ENDBDC= .FALSE. 

IF ( K5.NE.KBDY .OR. ENDCRD ) ENDBDC= . TRUE . 
IB DAT A I N= .FALSE. 

DBL PT S= .01 
JFOUND= 0 
CAPXi = 0. 

BL = .FALSE. 

C READ BDY INPUT CARDS 
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i5 I F ( t NOB DC ) GO TG AO 
FLIP = 1. 

RUT A I E = 0 . 

Z P I VG I = 0. 

RPIVU(= 0. 

SCALt = SCALEA 
ZTRANS= 0. 

RTRANS= 0. 

ZRONL Y= .FALSE. 

CALL SETM(l,.l, DEVI , 100 ) 

CALL SETK< 3 , 6 I T S , XA, 200, DE VA , BO , B , 300 ) 

CALL SETMI 1, BITS, X, 200) 

READ (5, A) 

IF { ZRGNLY ) CALL I SOR T ( XA , YA , DUM ,B , 200, 2 ) 

I F ( . NOT . Z RQNL Y ) CALL I SORT ( X, Y , ANGD , B , 300 , 1 ) 

IF( .NOT . ZRONL Y) CALL I SOR T l X , Y, ANGA , B,3D0) 

IF(INERR) ERRMA J = . TR UE • 

DA T A I N = .TRUE. 

RE S TR T = .FALSE. 

C COUNT THE LENGTH OF THE Z-LIST 
AO I F ( .NUT .DATA IN > GO TO 900 
I F { J FOUND .EQ.L ) GO TO A3 
N I =0 

00 Ai 1=1,100 

I F { X A ( I ). 60. BITS) GO TO A2 
Ai NI = I 

A 2 I F ( N I . EG . 0) GO TO A3 
LINES = 6 A 
CALL SMOTH 
J F OUN D= 1 
A3 NZ =0 

DO AS 1=1,100 
IF( Z< I ) .EQ.BITS ) GO TO 50 
A5 NZ = I 

50 IFINZ-2) 55,100,100 
55 WRITE (6,1055) BDY(l) 

ERRMA J= .TRUE. 

RETURN 

C DELETE DOUBLE POINTS FROM SMOOTH BOUNDARY RECORDS 
100 OMITFK= .TRUE. 

CALL FHE AD( NZ+ 10 ) 

WRITE (6,1090) IBDY,CHN, UPPER 

IF( JFOUND.NE.l .OR. DBLPTS.EQ.O. .OR. NZ.LE.2) GO TO 150 
WRITE (6,1100) DBLPTS,DBLPTS 

1 = 1 

110 I = 1+1 

I F ( I.GT.NZ ) GO TO 150 
IF ( ABS( Z( I )-Z( 1-1 ) KGE.DPTOL .OR. 
i ABS(R( I)-R(I-l) J.GE.DPTOL) GO TO 110 
ANGD I F= ABSIANGDl I )-ANGD( I-ll) 

IF (ANGOIF.GE.DBLPTS) GO TO 110 
NMOVE = NZ- I 

ANGSV = .5*( ANGO( I )+ANGD( 1-1 ) ) 

IF( ANGD( I )*ANGD( I-i) .EQ.O. .AND. ANGDIF .LE . .0005 ) ANGSV=0. 
ANGD ( 1 - 1 ) =ANGSV 

CALL MOVE ( 3 , Z ( I + 1 ) , l ( I ) , NMOVE , l ♦ 
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non oooooooooooooo 


1 R(I+1),R(I)« NMO VE , 1 1 

2 ANGD ( I + l ) ( ANGD ( I ) » NMOVE • 1 ) 

NZ = NZ-1 

GO TO 110 

C CALCULATE CURVATURES FOR PRINTOUT 
150 I =1 

CURVI 1 ) =0.0 
155 CURVB ( I ) = B I T S 

CUR V ( 1 + 1 ) =CUR V ( I ) 

DX = Z( 1 + 1 ) — Z ( I ) 

DY = R{ 1 + 1 )-R( I ) 

CHO = SORT ( DX*DX+DY*DY ) 

IFtCHD.LT.. 00001) GO TO 160 
ACHD = ATAN3 ( DY» DX, ANGD ( I ) *TQR AD ) 

YPA = ANG( I )*TORAD-ACHD 

YPB = ANGO( I +1 ) *TORAD- ACHD 

CURVB ( I )=(4.*YPA+2.*YPB)/(CHD*( l. + l. 5* YPA* YPA ) ) 

CUR V ( 1 + 1 )=(-2.*YPA-4.*YPB ) / ( CHD* ( 1 . + 1 . 5*YPB*YPB) ) 

GO TO 165 

160 I F ( I . EQ . 1 ) GO TO 165 

I F ( CURVB ( 1-1 ) .EQ.BITS) C UR VB ( I- 1 ) =CURVB ( I ) 

165 I = 1+1 

I F ( I.LT.NZ) GO TO 155 
CURVB ( I )=0.0 

♦REL013 RELOCATE FROM A ONE TO A THREE DIMENSIONED ARRAY -REL013- 

C SUBROUTINE REL013 

INPUT- 

Z , R = BOUNDARY COORDINATES 
ANGD = ANGLE OF THE BOUNDARY (DEGREES) 

NZ = NUMBER OF BOUNDARY COORDINATE POINTS 
FLIP = SCALER ON RID BEFORE ROTATION OR TRANSLATION 
ROT A T E = ANGULAR ROTATION IN DEGREES 

ZPI VOT,RPIVOT=PIVOT POINT FOR ROTATION BEFORE SCALING 
SCALE = MULTIPLICATIVE CONSTANT ON INPUT COORDINATES 
ZTRANS= Z-TRANSLAT ION AFTER SCALING 
RTRANS= R-TRANSLAT ION AFTER SCALING 
BDY = BOUNDARY NAME 

UPPER = T IF UPPER BOUNDARY, = F IF LOWER BOUNDARY 
CHN = CHANNEL NAME 

LBDE = NEXT AVAILABLE LOCATION IN THE BOUNDARY TABLE 
OUTPUT- 

BDT = TABLE OF Z »R , ANG IN 3-D ARRAY FORM 
LBDE = NEXT AVAILABLE LOCATION IN THE BOUNDARY TABLE 

IF ( EL IP.NE . 1. .OR. ROTATE. NE.O. .OR. SCALE. NE.l. .OR. ZTRANS.NE.O. 

1 .OR. RTRANS.NE.O. ) WRITE (6,1151) F L I P, ROTATE ,Z P I VOT ,RPI VOT , 

2 SCALE , ZTRANS,RTRANS 
WRITE (6,1152) 

LB 1 = LBDE 

LB2 = L Bl + 3* ( NZ - 1 ) 

LB = LB1 

BDT (LB )=BDY 
CHN AME ( LB )=CHN 
LBZi(LB)=0 
UP ( LB ) = UPPER 
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LEDEX(LB)=0 
I = i 

LBDEL = 3 
AOOPI = 0. 

IF ( .NOT .UPPER > GU TO 240 
LB = LB 2 
LBDEL = -3 
AODPI = PI 

240 ROTAT = ROT ATE* TUKAD 
SM = SIN(RGTAT) 

CS = COS(ROTAT) 

250 IFIROTATE.NE.O. ) GO TO 260 
ZBT(L6)=Z( UPSCALE + ZTRANS 
RBT (LB ) =R ( I )*FL IP*SCAL£ + RTRANS 
GO TO 270 

260 RFLP = R ( I ) *FL I P 

ZBT ( L B ) = ( ZPIVOT+CS*(Z( I ) - ZP I VOT )-SN* (RFLP-RPIVOT ) )*SCALE + ZTRANS 
RBT(LB)=<RPIVOT+CS*<RFLP-RPIVOT>+SN*(Z(I)-ZPIVOT> )*SCALE ♦ RTRANS 
270 ANGD( I ) =ANGD( I ) *FL IP + ROTATE 

ANGBT (LB)=ANGDm*TORAD + ADDPI 

WRITE (6,1280) l.ZBT(LB) , RBT ( LB > , ANGD ( I ) , CURV < I ) ,CURVB( I) 

IF ( I .GE .NZ ) GO TO 300 
I =1+1 

LB = LB+LBOEL 

GO TO 250 
300 LBDE = LB2 + 9 

LBNEXT ( LR1 ) =L BDE-L Bl 
BOT (LriDE )=BLANK 
C END SUBROUTINE REL013 


C SET UP BOUNDARY LAYER INPUT TABLE 

IBLB(LBLB)= I BDY 
IBLB( LBLB+1 )= 0 
I F ( BL) IBLB(LBLB+1)=1 
BL B ( L BL B + 2 ) = CAPX1 
LBLB = LBLB+3 
900 RETURN 

1055 FORMA T ( // IX43H** NO COORDINATE INPUT WAS FOUND FOR BDY = A6 » // ) 
1090 FORMAT ( // / IX,45HB OUNDARY COORDINATES, BDY = A6, 

* 5X4HCHN=A6 , 5X6HUPPER= 

*L2,6X,3HBL=L2, ) 

1100 FORMAT ( /6X46H DOUBLE POINTS WITH ANGLE DIFFERENCES LESS THANF6.3,lX 
*24HARE ELIMINATED ( DBLP T S=F 5 . 3 , 2H ) . ) 

1151 FORMA r ( /6X5HFL I P = F 7. 3, 3X7HR0 TATE =F 8. 3, 3X7HZPI VOT =F 1 0. 5 , 3X 7HRP I VOT = 
*F11 .5, 3X5HSCALEF7. 3, 3X7HZTRANS-F10.5, 3X7HRTR ANS=F 10*5, ) 

1152 FORMA T/9X48HI X,Z Y , K ANGD CURV- CURV+ ) 

12 80 FORMAT ( 1 1 0 , 2F 10 . 5 , F 10 . 3, 2F 10 . 4) 

END 
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♦deck stcb 

PROGRAM STCB 

COMMON /CHNFPT/ I CHN < 1 0 ) , WTF S 110 ) , WTF A QO ) , WPTO( 10 ) ♦ WTTO ( 10 ) * IC 
COMMON /SELECT/ LENTRY 
GO TO ! 10,20) , LENTRY 

C NORMAL ENTRY-- STATION LOOP, FLOW BALANCE 
10 CALL OVERLAY! 3HSTC, 2,1, 6HRECALL ) 

GO TO 30 

20 CALL OVERLAY! 3HSTC, 2, 2,6HRECALL) 

30 RETURN 
END 



oooooooooo 


-EDUMPX- 


♦ UECK FRRORX 

SUIWUGT INI tKKlJIU 

CLDUMPX I DUMP FOR STC EXECUTE SECTION 


C 

C 


C 


C 


COMMON /ALLCDM/ 

i 

R E AL 

LOGICAL 

LOGICAL 

CUM MO. 4 /CFG / 


* 


MACH A ,PSAtTSA»PTA,TTA, AX 1 A , RGA , GAMA , 

MACHC.PSC ,TSC»PTC,TTC, AX I C , RGC , G AMC » 
DAxITfSCALEAtTTE »CHOTST 
MACHA ( 1 ) ,MACHC 
A X I A t AX I C 
C HOT ST 

L» MA, MB ,PLB, PUB, WF ♦ CHOKE, SUBSON, NK » PLBC » PUBC « 
XCHOKE, TARfcA , VMBC , WRQST,WCALC, OV ( 8 ) , QVP ( 8 ) , 
J SUM » VML B SQ 


LOGICAL 

COMMON /ERASE / 
COMMON /ERASE2/ 


2 

2 


REAL 

DIMENSION 

EQUIVALENCE 

DIMENSION 

EQUIVALENCE 


CHOKE, SUBSON 

ERASEC (800) 

AREA (96) , ARE A0( 96 ) *DISP(96) , PT (96 ) , LAMBDA ( 96 ) , 
RHO( 96) , SQRTVV<96> ,TS(96) ,TT(96> ,VMSQ(96) , 

VVKQKPI96), 

WQA(96),WSTA(96), RG ( 96 ) ,C2CP (96 ) ,FGR(96) 
LAMBOA 

E S2( 96), SDNQRM(96) 

( ES2 , VVKQKP ) , ( SDNQRM ,RHO ) 

RCU( 96) 

(RCU, LAMBOA) 


FIELD TABLES 



INDEX- 

M=MO , NM 


COMMON 

/Cl 

/ 

Z( 300) 

COMMON 

/CR 

/ 

R ( 300 ) 

COMMON 

/CS2 

/ 

S2 ( 300 ) 

COMMON 

/CS1 

/ 

S I ( 300 ) 

COMMON 

/CPHI 1 

/ 

PHIK300) 

COMMON 

/CM 

/ 

JMS( 300) 

COMMON 

/CCURV 

/ 

CURVI300) 

COMMON 

/CB 

/ 

B( 300) 

COMMON 

/CIDEX 

/ 

M, J,MU,MD,ISTAG 

TABLE OF INDEX 

LI 

MITS 

COMMON 

/IXORIG/ 

LHO, LHE* LBDO.LBDE , LTO.LTE, LWO,LWE, LFO,LFE, 

* 



LO,LESTA »LSO*LSE,LDO,LDE,LDUM(A) , 

* 



MO , NM , N J ,NFCOLS, MAXN J , MAXOL , MA XNM , MAXLE , 

* 



LEO, LEE, LRO,LRE,LRD 

DIMENSION 


LIMITS(24) 

equivalence 


(LIMITS, LHO) 

COMMON 

/CVM 

/ 

VM ( 300 ) 

STREAMLINE TABLE 


COMMON 

/SLTAB 

/ 

W( 128) ,X2( 128), SLCHN (128) 


INTEGER SLCHN 
BOUNDARY TABLE 
INDEX- LB=LBDO,LBDE 
LBNEXT= INCREMENT TO NEXT BOUNDARY 
LBZi = INCREMENT TO THE FIRST BOUNDARY POINT ( 

CHNAME= CHANNEL WITH WHICH THE BOUNDARY DATA IS 
UP = T OR F FOR UPPER OR LOWER BOUNDARY 

LEDEX = RELATIVE INDEX OF L.E. POINT WHEN LOWER AND UPPER SURFACE 
CONTOURS ARE CONNECTED 

BDNAME,LBA,LBB=NAME AND INDEX LIMITS OF SPECIFIC BOUNDARY 

DATA WHEN BOUNDARIES ARE COALLATED 


0 BEFORE COALLATIO 
ASSOCIATED 
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1 

2 


DIMENSION BDT ( 1 ) ,LBNEXT (I) ,LBZ 1 ( 1 ) , 

CHNAME d ) ,UPd) tLEDEX(l) , 

Z B T ( 1) ,RBT( 1 ) » ANGST ( 42 ) 

LOGICAL UP 

INTEGER BDT, CHNAME, BDNAME 

DIMENSION BDNAME ( 1 ) »LBA (1 ) , LBB (1 ) 

EQUIVALENCE < BDNAME , ZBT ) , ( LBA , RBT ) t ( LBB t ANGBT ) 

FLOW ADJUSTMENT TABLE 
INDEX- L F = L FO , L FE 
NFCOL S= B 

X 1 F = ORTHOGONAL COORDINATE 

X 2 F = STREAMLINE COORDINATE OF SL EMINATING FROM T.E. 

XIBF = Xl-COORDINATE OF CHOKE STATION OF FLOW BELOW T.E. 

X I A F = Xl-COORDINATE OF CHOKE STATION OF FLOW ABOVE T.E. 

S1F = SI-COORDINATE OF T.E. (UPPER SURFACE). THIS ITEM 
IS USED WHEN INTERPOLATING FOR WAKE DELTA-STAR. 

LFB,L FA= INDICES OF STATIONS BELOW AND ABOVE T.E. 

NCHB,NCHA=NUMBER OF CHANNELS BELOW ANO ABOVE T.E. 

LRF = INDEX OF DUMMY ORTCHN LIST FOR THE T.E. 

LRXF = INDEX OF LAST CHANNEL BELOW THE T.E. 

JORDER= 0 IF TOTAL FLOW AT X1F IS GIVEN 
= 2 IF FLOW ABOVE T.E. IS GIVEN 
= 1 IF FLOW BELOW T.E. IS GIVEN 
JORDER= -1 IF FLOW AT XIF IS CHOKED AND SINGLE CHANNEL 
DIMENSION X1F( I) ,X2F(1) #XIBF(1) tXlAF(l) , 

1 S 1 F ( I ) . NCHB ( 1 ) , NCHA( I ) , JORDER ( 1 ) »VNR(I2) 

EQUIVALENCE ( LFB ,X1BF ),( LFA.X1AF ), (LRF ,NChB> , (LRXF, NCHA) 
DIMENSION LFBC I),LFA(I),LRF(1) ,LRXF(1) 

STATION TABLE 
INDEX- L =LO , L E S TA 

S CHUK E- STATION CHOKE INDICATOR ( ADJ WF , BRHS , WR IOUT ) 

MCL = SHARP CORNER INDICATOR ( BLDTBS ) 

MCL = FIELD INDEX OF CONTROL STREAMLINE ( PTMOVE , FLOBAL ) 

COMMON /CHDATA/ X 1 ( I ) , LNE XT ( i ) , MLB ( 1 ) , MUB ( I ) , PRI M ( 1 ) , 

1 TYPELBd ) ,NAMELB( I ) , ILB( I) ,FLB(1 ) ,S1LB( 1 ), 

i TYPEUBd ),NAMEUB( 1), IUBd) ,FUB(i ) ,SIUB( I) , 

3 VMB( 1 ) ,DWDV(1), X2CL ( I ) » VCL ( I ) » MCL (48 1 ) 

LOGICAL PRIM 

DIMENSION SC HOKE ( I ) 

EQUIVALENCE ( SCHOKE , DWDV ) 

EQUIVALENCE ( B DT , X IF , X L ) , ( LBNE XT , X2F , LNEXT ) , ( LBZI , XIBF , MLB ) 

EQUIVALENCE (CHNAME ,XIAF, MUB ) , ( UP, S1F , PRI M) 

EQUIVALENCE ( LEDEX,NCHB,TYPELB) , (ZBT, NCHA, NAMELB) 

EQUIVALENCE ( R BT , JORDER, I LB ) , ( ANGBT ,VNR,FLB) 

COMMON /CTABPR/ UTAB 
COMMON /BLBDY / IBLB(60) 

CALL TABPRT(6HALLCOM,MACHA,2J,8) 

CALL TABPRT ( 3HCFB,L» 33,4 ) 

CALL TABPRT(bHCIDEX,M,5,5) 

I 1 T Art = LBDO 

CALL TABPRT(6HBDYTAB,BDT ,LBDE ,3 ) 

I1TAB = LFO 

CALL TABPRT(6HCADJWF,X1F,LFE,8) 

UTAH = LO 

CALL TABPRT ( 6HSTATAB»X1»LESTA,5) 
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1 -Jj W«Ill it,, 1 1 *J» 0 I ( Jf X2U) fSLCHMI J) tW( J) tJ«l*NJ) 

CALL ) ABPRT { ‘>HLRASL,ERASLC,800, 5) 

CALL JMSt'KT 

CALL I AHPRr<2HSl,51,NM,Ij) 

CALI JAliPRT < 2HS2,S2,NM,10) 

CALL I ABPRT < 1HZ,Z ,NM,10) 

CALL I ABPRT ( 1 HR , K , NM , 10 ) 

CALL TABPRK4HPHI 1,PHI l.NM, 1 J) 

CALL TA8PRT(4HCURV,CURV,NM,10) 

CALL TA8PRT ( 2HVM, VM,NM,10) 

CALL IA6PRT ( 1HB»B»NM*10) 

CALL T ABPRT ( 6HER A SE2, ARE A ,1536,8) 

IK IBLBID.NE.O ) CALL TABPRT < 5HBLBDY, I BLB , 60 ,3 ) 
IK LUE.LG.O ) GO TO 1321 
I 1 T Ah = LDO 

CALL TABPRT (5HBLT AB, CHNA M,LDt ,3 ) 

1321 CONTINUE 

LSTGP = 5 

GO U) (999,999) , LSTOP 
999 RETURN 

1182 FORMA T ( ///1X17HSTREAMLINE TABLE-/ 17X32HJ X2 

* /./< IlB,F12.6,6X,A6,F12.6 f ), ) 

EMO 


SLCHN 
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♦DECK STCwl 

PROGRAM STCW1 


WRITE 

the over-all 

CALL 

WR I A 

CALL 

WR I OUT 

CALL 

WR IBDY 

CALL 

WR I ATP 

RETURN 

END 



STC OATA RECORD, 


KEY(5)=A. 
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nooooo oo onooooooooo 


• ni.CK WK1A 

SUBKUUT INE WRIA 

♦WKIA-- WRITE THE K E Y ( b ) *A STC DATA RECORD -WRIA- 


CilMKfJfM /ALL COM/ 


RLAL 

LOGICAL 

LOGICAL 

TEAL 


MACH A, PSA, TSA,PTA,TTA, AXI A ,RGA, GAMA, 

MACHC.PSC ,TSC,PTC,TTC, AX IC ,RGC, GAMC , 

DAXIT , SC A LEA, TTE,C HOIST 

MACH A ( 1 ) , MACHC 

AX IA, AXIC 

CHOTST 

MACHO 


tOUI 

COMMO 

1 

INTE 
LOGI 
SSFML 
SSLF 
S S E AN 
SSDF 
$$FE\ 
SSFNO 
SSOLt 
A4F AC 
BRL X 
CJRRL 
TS IC 
COMMO 

* 


dime 

lQUI 

COMMO 

INTE 


(MACHO,MACHA) , ( PSO.PSA ) , ( T SO, TSA ) , 

(AXI , AXIA ) , (RG» RGA ) , ( GAM, GAMA) 

S S FML, SS E F, S SEANG, SSDF ,SSP END, SSFNDi 
, S SDL t » A4F ACT , BRLX»CURRLX»TS IC »RHOC , RHOCSS 
SSFML 

SSEF , SSDF, SSOLE 

= SUPERSONIC CURVATURE FORMULA NUMBER 
= SUPERSONIC ENTERING FLOW* T OR F 
G« ENTERING FLOW ANGLE (DEGREES) FOR SSEF*T 
= SUPERSONIC DISCHARGE FLOW, T OR F 
D= SUPERSONIC BEAM DOWNSTREAM END CONDITION, =0,1 FOR PARABDL 
.= SUPERSONIC BEAM UPSTREAM END CONDITION, =0,1, FOR PARABOLA 
= SS FLOW BELOW ANO AF I OF LE PT, T OR F 
1= CENTRAL POINT INFLUENCE COEFFICIENT FACTOR 
= B-RtL AXAT ION FACTOR 
X= CURVATURE RELAXATION FACTOR 
= NUMBER OF POINTS 10 BE READ FOR TRANSONIC INTERPOLATION 
N /IXORIG/ LHO, LHE , L8D0,LBDE , LTO,LTE , LWO,LWE, LFO,LFE, 
LO,LESTA, LDUM ( 8 ) , 


VALENCE 

N /CSS / 

oER 

CAL 


NS ION 
VALENCE 
4 /SLTAB / 
GER SLCHN 


MU , NM, N J »NFCOLS, MAXNJ»MAXOL,MAXNM,MAXLE, 
LEO, LEE, LRO,LRE,LRO 
LIMITS (24) 

(LIMITS, LHO) 

W( 128) ,X 2 (128), SLCHN (128) 


COMMON /BCOMMN/ 
LOGICAL 

COMMON /ADAMQI/ 
COMMON /BENDIN/ 
COMMON /CB / 
COMMON /CBITS / 
COMMON /CCRX / 
DIMENSION 
EQUIVALENCE 
CHANNEL INPUT 


PR0GM19) , FILIN, FILOT 
FILIN, FILOT 

NAME ( 6 ) , ADDRL S (6), TITLE (6), IDENT (6) 
N BC I N ( 2 ) , ACF ( 2 ) 

B( 300) 

BITStlBLANK 

CRXSL , CR XOL , CRX SS ,CRXE ,CRXC ,CRMACH 
CR X ( 6 ) 

(CRX, CRXSL) 

DATA TABLE 


INDEX - LH=LHQ»LHE 

COMMON /CHDAT A/ CHNAM ( 1) , LHNt XT ( i ) , WTFLOW (10 ) ,NR ( 1 ) , NC ( I ) ,TAB( 6 ) , 
I BB ( 7 5 ) 

DIMENSION TABLES(998), TTO(I),PTQ(l) 

TARLE OF CONVECTED PROPERTIES 
INDEX- LT=LTO,LTE 


CH = CHANNELNAME 

LTNEXT= INDEX INCREMENT TO THE NEXT- CHANNEL 
LPSI = RELATIVE LOCATION OF PS I LIST 
NPT = NO. OF P SI, TT, PT AND RCU VALUES 
DIME NS ION CH ( I ) » LTNEXT ( i) ,NPT ( 1 ) ,LPS I (1 ) *LTT(495) 
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IN1LGER CH » C HN AM 
DIMENSION XCHl l ) 

fcOU I VALENCE ( C HN AM , TA BLE S » CH ♦ XCH ) t (LHNE XT , LTNE XT ) » 

* ( W T FLOW t NP T ) « 

* (WTFLOW(2 ) ♦LPSItTTO) , 

* ( W T FLOW ( 3 ) » LT T * P TO ) 

COMMON /CIADIN/ RHOBAS ,RHOAMP, I AOM 

COMMON /CINNER/ I NRC TR , RDUM , N INNER l 1 6 ) ,CNVF ( l 6 ) 

COMMON /CISBOT/ FARFLD ( 2 ) ,FREE ( 2 > ♦ PRES ( 2 ) , RFF ,NZP » 
i 2P ( 10) ,PPSI 10), A L * A2 1 ADUM( 6) 

INTEGER FARFLD,FREE, PRES 

COMMON /CLINES/ L I NE S ♦ OM T IFK , PT I TLE ( 6 ) 

LOGICAL OMTIFK 

COMMON /CM / JMSI300) 

COMMON /CM AX IT/ M AX I T , MA JCTR , GREF I N, EDUM 
LOGICAL GREFIN 

EQUIVALENCE ( MA JC TR , NREF I N ) 

COMMON / CNORM / R HL . RM, AHL , A RM 

COMMON /CPI / PI ,TWOPI ,PIQ2,PIQ4,T0DEG,T0RAD 

COMMON /CPRINT/ PDUM 1 ( 3 ) , PRE F IN , PREFN2 > PDUM ( L l ) 

COMMON /CPRPRN/ PRPRN 
INTEGER PRPRN 

COMMON /CPTMOV/ VELPOT » I COB»NODENS »CPTDUM 
LOGICAL VELPOT 

COMMON /CR / RF ( 300 ) 

COMMON /CREFIN/ SLS, SG21 , VMG1 ,VMG2 , NGR,NGZ,SGR( 10) ,GR( 10 ) , 
i SGZ( 1 0 I » G Z ( 10 ) 

COMMON /CS1 / S 1 ( 300 ) 

COMMON /CS2 / S2 ( 300 ) 

COMMON /CTOLRL/ TOLRL , MA X SHP , CL EN* DS2MX , T0LES2 ,N SWP , 

1 0S10MP,DSIMXA,DS1MXB,DS1RMS,ES2MX,DS1RM0* 

c SGIMIN.TOUNR 

COMMON /CTHICK/ N THK X , NTHKV , THKX ( 20 ) , THKY I 20 ) , TH I K2D< 7B ) 

COMMON / C VM / VMF ( 300 ) 

COMMON /CZ / ZF ( 300 ) 

COMMON /CHNEPT/ I CHN < 10 ) , WTF S <10 ) , WTFA ( 10 ) , WPTOl 1 0 > , WTTO ( 10 ) , IC 

COMMON /TAPES / NTAPO,NTAPN 
COMMON /BLBDY / BLB<60> 

DIMENSION I BL B ( 60 ) 

EQUIVALENCE ( I BLB * BLB ) 

COMMON /VISCOS/ T REF , MURE F , SCON 
REAL MUREF 

LOGICAL STCF1L 
DATA STCFIL/T/ 

DATA KA/1HA/ 


OMT I F K= .TRUE. 

I F ( F I LOT ) OMT IFK=. FALSE. 

CALL EHEAD( 6A ) 

T SC = TSA 

TTC = TSC*< l.+( GAMA-1. ) * . 5*MACH0*MACH0) 

PTC = PSC*(TTC/TSC)**( GAMA/ ( GAMA-1 . ) ) 

‘3 6 WRITE (6,1 000 ) A X I , MACHO , RG, T SC , GAM , P SC * TTE , PTC tCHOTST, TTC , 

1 NBC IN* ACF 

1000 FORMAT (/15H GENERAL INPUT-// 6X,7HAXI = , L8 * 26 X , 7HMACH0 *,F8.4/ 
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16X./HRG =,F8.2,26X,7HTS0 *,F8.2/ 6X,7HGAM *,F8. 4, 2bX, 

* 7MP SO = » F 8 . 3 / 6X * 7HTT6 = , F 8 . 3 » 26X , 7HPT 0 =,F8.2/ 6X,7HCHOTST= 

. ,Lb,.6X, 7HTTU =,F8.3// 27H STREAMLINE END CONDITIONS-/ 6X.7HNBC 

i I N =,218/ 6X.7H4CF =,2F8.3/ ) 

WRITt(6,l005) SSFML,SSFENO,SbFNDl,SSEANG,SSEF,SSDF,SSDLE 
1005 FZ)RMA1(43H CURVATURE CALCULATION FOR SUPERSONIC FLOW-/ 

16X , 7HSSFML =,I8,19H (FORMULA NUMBER I / 

c6X,7HSSFtND=,F8.3,43H (DOWNSTREAM END CONDITION, SSFML*2 ONLY)/ 
36X, 7HSSFND1=,F8.3,41H (UPSTREAM END CONDITION, SSFML*2 ONLY)/ 
46X, 7H S S E ANG= , F 8 . 3 » 43H (INLET FLOW ANGLE, DEGREES, SSEF-T ONLY)/ 

b / 3 8 h SUKSON IC/SUPtRSONIC BRANCH SELECTION-/ 

66X , 7HSSEF =,L8,37H (SUPERSONIC ENTERING FLOW, T OR F)/ 
76x,7hSSDF =,L8,56H (SUPERSONIC FLOW DOWNSTREAM OF CHOKE STATION 
*, T OR F)/ 6X, 7HSSDLE =,L8,58H (SUPERSONIC FLOW BELOW AND 

BAFT OF A L.E. POINT, T OR F) ) 

WRilr (6,1010) (GK ( I ) , 1=1, NGR ) 

KRITl (6,1011) ( SGR ( I ) » I = 1 ,NGR ) 

IF(NGZ.EC.O) GO TO 65 

WRITE (6,1012) ( GZ ( I ),I = 1,NGZ) 

WRITE ( 6 , 1C 1 3 ) ( SGZ ( I ) » I = 1 ,NGZ ) 

65 WRITE (6,1014) VMG1, VMG2 ,CRX 

1010 FORMAT ( /1X19HGRID SIZE CR ITER IA-/6X7HNGR/GR*10F8 « 2 ) 

1011 FORMAT (6X,7HSGR =,10F8.2) 

1012 FORMAT ( / 6X , 7HNGZ /GZ= , 10F 8. 2 ) 

1013 FORMAT ( 6X, 7HSGZ =,10F8.2) 

10.4 FORMAT ( /6X,7HVMG1 =,F8. 2 , 25X,7HVKG2 * ,F8 . 2/ /6X , 7HCRX *,6F8.3) 

WRITE (6,1030) NM,MAXNM, LESTA,MAXLE, NJ,MAXNj 
10 .30 FORMAT ( /1X19HMEMQRY UT IL I ZATI0N-/24X17HUSED A VAI LABLE/ 6X 1 1 HGRI D 

* P0|NTSI11,I10,/ 6X6HTABL ES 1 16 , 1 10, /6X1 1HSTREAMLI NESI11, 110, ) 

ATLDS2= CLEN* TOLE S2 

WRITt (6, 1040) MAXIT.NREFIN, I NRC TR , TOL I NR , TOLE S2 , C LEN , ATLDS2 , ES 2MX , 
1 DS 1DMP , NODENS 

104u FORMAT </l8H CONVERGENCE DATA-/ 

1 6X , 7HMAX I T = , 1 8 , 3 X , 20H ( M A X I MUM ITERATIONS)/ 

26X,7hNREFIN=, I8,34H - NUMBER OF REFINEMENT ITERATIONS/ 6X,7HINRCTR 
3= , 18, 56H - NUMBER OF ADDITIONAL ITERATIONS AFTER LAST REFINEMENT// 
46X, 7HT0LINR=, E8.1,47H (INNER ITERATION TOLERANCE ON S.L. MOVEM 
4E NT ) / 6X , 7HTULE S2 = » E8.1,37H (FINAL TOLERANCE ON S.L. MOVEMENT) 

o/6X,7HCLEN =, F8 . 3 , 52H - CHARACTERISTIC LENGTH BASED ON GRID SIZ 
6E CRITERIA/ E21.1,38H - ABSOLUTE TOLERANCE ON S.L. MOVEMENT/ 
76X,7HMAXES2=, E8.1,42H - LARGEST S.L. MOVEMENT ON LAST ITERATION/ 
8/6X, 7 HD Si DM P®, F8.3,54H (STREAMWISE PT MOVEMENT DAMPING, *0 FOR 

9 NU DAMPING)/ 6X , 7HN0DEN S = , I 8 ,58H (REFINEMENT LEVEL TO WHICH CON 

ASTANT DENSITY IS ASSUMED)) 

LINES = 64 
CALL F H F A D ( 13) 

WRITE (6,1090) FARFLD 

WRI1L (6,1092) IADM,RHOBAS,RHOAMP,TOLRL 
109j FORMAT ( / 2 6H SPECIAL BOUNDARY OPTIONS-/ 6X , 7HFARFLD* ,2 ( 2X , A6 ) ) 

1092 FORMAT ( / 28H MATRIX SOLUTION PARAMETERS-/6X , 7H I ADM * , I 8 , 3X , 70H ( ■- 
11,0,1, FOR STREAMLINE, ALTERNATING, AND ORTHOGONAL LINE RELAXATION 
2)/ t>X , 7 HR HO BAS=,F8.3,3X, 3 3H ( ACCELERATION FACTOR, BASE LEVEL)/ 

36 X , 7URH0AMP = ,F8.3,3X,45H(ACCtLERATI0N FACTOR, AMPLITUDE OF VAR I AT I 
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'VON ) / 6X , 7HTOLRL =, E8. 1 , 3X , 30H{ TOLERANCE RELATIVE TOMAXDS2) ) 

C PRINT HIGHLIGHT AND MAX. BODY RADII AND AREAS 
AHL = RHL 

IF(AXIA) AHL=PI*RHL*RHL 
ARM = RM 

IF(AXIA) ARM=P I *RM*RM 

WRITE 16,1091) RHL , AHL, RM, ARM 

1091 FORMAT (//6X, 1 7HH IGHL IGHT RADI US = t F 8. 3, AX , 15HHI GHL I GHT ARE A= , 

* F8.3/6X, 17HMAX. BODY RADIUS*, F8 . 3 , AX , 1 5HMAX . BODY AREA=,F8.3) 

C PRINT CHANNEL TABLE OF CONTENTS 
CALL FHE AD ( 2 ) 

WRITE (6,1060) 

LH = LHO 

80 IF ( LH.GE.LHE ) GO TO 96 
MOREL * A 

IF(NRILH) . N E • 0 ) MORE L=M0REL+2+NR ( LH ) 

CALL FHE AD{ MOREL ) 

LH2 = LH+9 

WRITE (6, 1070) CHNAM(LH) , IWTFLOW(LHX)*LHX=LH» LH2 ) 

NCX = NC(LH) 

IF(NR(LH) .LE.O) GO TO 95 
WRITE (6,1080) ( T A B ( I), 1=1, NCX) 

CALL TA8PRTI 2HB=,BB( LH) ,NCX*NRI LH) ,NCX) 

95 LH = LH+LHNEXT ( LH ) 

GO TO 80 

96 CONTINUE 

1060 FORMAT < / 1 X26HC0NTENT S OF CHANNEL TABLE-) 

1070 FORMAT ( //6X7HCHN =2X,A6,5X7HWTFL0W= El 2 . A, /6X7HTT0 = F8.2t5X 

♦7HPTG =F8 . 3 , 5X7HTS0 =F8. 2 ,5X7HPSO -F8. 3 , /6 X7HMACH0 =F8.A,5X7 
*H AO = E12.A,1X7HVARY =L8,/6X7HRG = F8.2.5X7HGAM »F8.A,) 

10 aO FORMAT ( /6X7HN B/ TA B=2 X , A6 , IH , 5X, A6 , 1H , 5X, A6 , 1H , 5X , A6 , 1H,5X,A6,1H, ) 

C LOOP THROUGH CHANNELS TO PRINT FLOW RATES, PRESSURES, AND TEMP 

RHO 1 N F = PSA/ ( RGA*TSA ) 

VINF = SQRT(GAMA*RGA*TSA )*MACHA 
WT NQRM= RHOINF*VlNF*PI 
J 2 =0 

IC =0 

100 J2 = J2+1 

JCHN = S L C HN ( J 2 ) 

1 j5 I F ( JCHN.NE.SLCHN( J2+1) .OR. J2.EQ.NJ) GO TO 110 
J 2 = J2+1 

GO TO 105 
110 IC = IC+1 

WTF A ( 1C >=W( J2 )/WTNORM 
I F ( RGA.NE.l. ) WTFA( IC) =W( J2) 

I CHN ( IC )=JCHN 
LT = LTO 

115 IF ( JCHN.EQ.CH(LT ) ) GO TO 120 
LT P = LT + LTNEXT ( LT ) 

IF ( LT P.GE .LTE ) GO TO 120 
LT = L TP 

GO TO 115 

120 LTP = LT + LPSI ( L T ) +NPT < L T) -i. 

WT F S ( IC ) = XCH( LTP )/WTNORM 
I F ( RGA.NE.l. ) WTFS( IC)=XCHILTP) 

WPTO( IC ) = PTC 
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WTTUl IC ) =TTC 
LH = L HO 

122 IF ( JCHN.I C.CHNAM(LH) ) GO TO 12A 
LHP = L H+LIINC XT ( LH ) 

IMLHP.GL.LHU GO TO 128 
LH = LHP 

gu ru 122 

12A IFCPffjCLHJ.Nfr.BlTS .AMO. PTO I LH ) . NE . 0 . ) WPTO( IC) *PTO( LH) 

IF ( TICK LH > .ML. BITS .AND. TTO ( LH ) . ME . 0. ) WTTOI IC) *TTO( LH ) 

128 IFIJ2.LT. NJ) GO TO 100 

110 WRITE (6,1130) (ICHN(I)»WTFS(l)»WTFA(I),WPTO(I)*WTTO(I),I=l»IC) 
1130 FORMAT </49H CHANNEL FLOW RATfcSt PRESSURES, ANO TEMPERATURES-// 

* 16X, ^SPECIFIED, 5X,8HADJUSTfcD,7X,6HPT/PS0,7X,6HTT/TS0 / 

* ( 6X, A6,AF13.4, ) , ) 


RETURN 


ENTRY WRIATP 
REWIND NTAPN 

WRITE (NTAPN) STCFIL, (LIMITS! I) ,1=1,24) 

WRI T t (NTAPN) ( ( I DENT ( I ) , I = 1 , 6 ) , AX I , RG ,GAM , MACHO , PSO, TSO, PTO, TTO , 
x PRPRN,TTE,CH0TST,MAXIT,MAJC1R, ( NINNER (I), 1*1,16) , VELPOT , ICOb, 

2 NOUtNS,RN,NGR,NGZ, CSGR( I ) ♦ I * 1 ,40 ) , VMG1 , VMG2 , INRCTR, SLS,SG21, 

3 NBC. IN ( 1 ) , NBC I N ( 2 ) , ACF ( 1),ACF(2>, SSFML , SSEF , SSE ANG, SSOF , SSFENO, 
A SSFNO 1 , SSDL E , A4F ACT , 8RL X, CURRLX* TSIC, (FARFLD(I ) ,1 = 1,8) , 

♦ RHUC,RHOCSS,RHL,RM, 

* T REF, MURE F, SCON, ( BLB ( I), I *1,60), 

0 (ZP( I ), 1 = 1,28), (TABLES! I), 1*1, LESTA), I B ( I ) , I *1 , NM ) , (JMS(I), 

6 1*1, NM), (Si( I),I = 1,NM), (S2(I),I*1,NM),(2F(I),I=1,NM), ( RF ( I ) , 

7 1*1, NM), ( VMF( 1 ),I=i,NM), ( W( I ) , I * 1 , N J ) , (X2(I ) • I * 1 ,N J ) • 

O ( SLCHNl I ),I=1,NJ),T0LRL,MAXSWP,T0LES2,T0LINR,SG1MIN,DS1DMP, 

A DS1RM0, (CRX( I ) , 1*1,6), RHOBAS ,RHQAMP, I ADM ,NTHKX , NTHKY , 

B ( THKX( I ), 1=1,118) ) 

NTSAV = NTAPO 
NTAPO = NTAPN 
NTAPN = NTSAV 
RETURN 
END 
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♦ DECK USECDW 

BLOCK DATA USECDW 


♦USECDW 

RE PL 

ACE 

STCW USE CARDS 

COMMON 

/ERASE 

3/ 

WDUM (318 ) 

COMMON 

/CP SM 

/ 

PSM{ 768) 


EMD 
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♦ DICK w < i un i 

S ( J B R 0 U T I N f wKIUUI 

♦wkiuut hkirt stc output data -wriout- 


CUMMUN /ALLCUM/ 


REAL 

LOGICAL 

LOGICAL 

COMMON /CFB / 


LOGICAL 

COMMON /CSS / 
INTlGLR 


MACHA,PSA,TSA,PT A, TTA, AX I A , RGA , G AMA , 

MACHC,PSC »TSC»PTC»TTC, AX I C ,RGC, GAMC , 

UAXIT , SC ALE A , TTEtCHOTST 
MACHAI 1 ) f MACHC 
AX IA ,AXIC 
CHOTST 

L,MA,MB, PLB,PUB,WF , CHOKE , SUBSQN, NK , PLBC , PUBC , 
XCHOKE, T ARE A , VMBC , WRQST,WCALC, QV ( 8 ) , QVP ( 8 > , 
J SUM » VML B SO 

CHOKEtSUBSON 

SSFML , SSEF, SSEANG, SSDF , SSF END , SS F ND1 

SSDLE ,A4F ACT, BR LX, CURRLX 

SSFML 


LOGICAL SSEF, SSDF, SSDLE 

SSFML = SUPERSONIC CURVATURE FORMULA NUMBER 
SSEF = SUPERSONIC ENTERING FLOW, T OR F 
SSEANv,= ENTERING FLOW ANGLE {DEGREES) FOR SSEF = T 
SSDF = SUPERSONIC DISCHARGE FLOW, T OR F 

SSFLNQ= SUPERSONIC BEAM OOWNS TRE AM END CONDITION, =J,I FOR PARABOL 

SSFND.= SUPERSONIC BEAM UPSTREAM END CONDITION, =0,1, FOR PARABOLA 

SSDLE = SS FLOW BELOW AND AFT OF LE PT, T OR F 

AA F AC T = CENTRAL POINT INFLUENCE COEFFICIENT FACTOR 

BRL X = B-REL AXAT ION FACTOR 

CURRLX= CURVATURE RELAXATION FACTOR 

COMMON /ERASE 2/ AREA (96) , ARE AO ( 96 ) ,D I SP ( 96 ) , PT (96 ), LAMBDA ( 96 ) , 

1 RHO( 96) , SQRTVV< 96 ) , T S ( 96 ) , TT ( 96 ) , VMSQ(96) , 

VVKQKP ( 96 ) , 

2 WCJA(96), WSTA(96), RG ( 96 ) ,C2CP ( 96 ) , FGR ( 96 ) 


C 

c 


N E k 
CAN 


REAL LAMBDA 

COMMON /ERASE 3/ J 1 ( 1 0 ) , K 1 ( 10 ) ,CHANLS( 10) ,PS(96) , MACH (96) , FLOW (96) 
DIM- NS ION X I2( 96 ) , Z (96) ,R (96) ,PHI ( 96) ,CURV(96) ,PSOPO(96) , 

» VM(96) ,FVX(96),FVY(96) ,FPX(96) ,FPY(96) ,SVX(96i , 

SVY( 46) , SPX( 96) ,SPY(96) »STX(96) , STY(96) 

(ARE AQ ,XI2,FVX,STX), (DISP,Z,FVY,STY), 

( SQRTVV,R ,FPX) , ( VMSG, PHI , F PY ) , ( VVKQKP ,CURV , SVX ) , 

( WQA,PSQPO,SVY) , (C2CP,VM,SPX ) , 

MACH 

X( 1) ,Y( 1 ) 

(X,Z), ( Y , R ) 

VARIABLES FOR NASA VERSION ONLY 
USE FGR IF NEEDED 


E UU l VALENCE 


REAL 

DIMENSION 
EGUI VALENCE 


(FLOW, SPY) 


DIMENSION 


EQUIVALENCE 


P F LO W ( 96 ) ,PSQPT (96) ,TSQTT(96) ,CP(96 ) , AQAREF ( 96 ) , 
PTQPTQ(96) , FLOWMX ( 10) 

( FLOW,PFLOW) , (LAMBDA, PSQPT) , ( T S , TSQTT ) , 

( R HO , C P ) , (FGR, AQAREF ) , ( RG, PTQPTO) 

COMMON /IXORIG/ L HO, LHE , LBOG , LBDE , LTO,LTE , LWO,LWE, LFO, LFE , 

LO,LESTA, L DUM ( 8 ) , 

MO ,NM , nj,nfcols, MA XN J , MA XOL , MA XNM , MAXLE , 
LlO,LEE, LRO , LRE , LRD 
L IMI TS( 24) 

( LIMITS, L HO) 

NBCBU),ANGE(c) , CURVE (2) ,FB(2) 

BITS, BLANK 


D I Mt NS ION 

equivalence 

COMMON /CBEND / 
COMMON /CHITS 


/ 


131 - 



o o o o o 


CUMMuj /CCUBE / NBC ( 2) *C 1(2) » C 2 C 2 ) » FEND (2 ) 

COMMON /CGRAV / CG 

COMMON /CPI / P I , TWO PI ,PIQ2»PIQ4,T ODEG » TOR AD 

COMMON /CREFIN/ SG I , SG2, VMG1 , VMG2 

1, NGR.NGZ, SGR( 10 ) ,GR( 10) , SGZ ( 101 , GZ ( 10 ) 

COMMON /SL TAB / W ( 12 8 I , X 2 ( 128 ) , SLCHN ( 128 ) 

INTEGER SLCHN 
STATION TABLE 
INDEX- L=LO»LESTA 

SCHOK E= STATION CHOKE INDICATOR ( ADJWF ,BRHS,WRIOUT ) 

MCL = SHARP CORNER INDICATOR (BLDTBS) 

MCL * FIELD INDEX OF CONTROL STREAMLINE ( PTMOVE ,FLOBAL ) 

COMMON /CHDATA/ X 1 ( l ) , LNE XT ( 1 ) , MLB I i ) , MUB I II • PRI M (1) » 

1 TYPELB ( 1 ) »NAMELB ( 1 ) « ILB(L) .FLB (I) ,S1LB( I > , 

i T YPEUB (1)»NAMEUB(1)» IUB (l)»FUB(l), S IUB (II* 

3 VMB( 1) ,DWDV(1)» X2CL (1) «VCL(1 ) »MCL(48I) 

LOGICAL PRIM 

INTEGER TYPELB, TYPEUB 
DIMENSION SCHOKE(I) 

EQUIVALENCE ( SCHOKE , DWDV ) 

COMMON / BCOMMN/ PROGM ( 9 ) , F IL I N, F I LOT 
LOGICAL FILIN, FILOT 

COMMON /ADAM01/ N AME ( 6 ) , AODRE SI 6 ) • T I TLE ( 6 1 * I DENT ( 6 ) 

COMMON /CCURV / CURVF ( 300 ) 

COMMON /CDS2 / MACHM ( 300 ) 

REAL MACHM 

COMMON /CPHli / PHIK300) 

COMMON /CPRINT/ PRTES2,PRTB,PRTA f PREFI N, PREFN2 ,SSONI C ,PDUM( ID I 
COMMON /CPRPRN/ PRPRN 
INTEGER PRPRN 

COMMON /CPSM / PSM(300) 

COMMON /CS2 / PTMI300) 

COMMON /CR / RFI300) 

COMMON /CRHS / T TM ( 300 ) 

COMMON /CVM / VMFI300) 

COMMON /CZ / Z F ( 300 ) 

COMMON /CIDEX / M , J , MU , MD , I S T AG 
COMMON /CLINES/ L I NE S , OM I TFK , PT I TLE ( 6 ) 

COMMON /CFRFLD/ FSAV(300), S TXU ( l 28 ) , STXD ( 1 20 ) , S T YU { 128 ) , STYD ( 1 28 ) 
COMMON /CHNFPT/ I CHN I 1 0 ) , WTF S 1 1 0 ) , WTF A ( 10 ) , WPTOU 0 ) , WTTO (10 ) , IC 

INTEGER DBSTARfSUB, SUPE R , BLANK , BRANCH, CHANLS • AS TE RP • TE 
LOGICAL UPSTRM,DNSTRM 

DATA TE/2HTL/ 

PIINV = I. /PI 
QO =0. 

IF ( MACHA.LE . . I ) GO TO 95 
IFIGAMA.NE.O. ) GO TO 92 
QO * I RGA*TSA) / ( PSA*MACHA*MACHA > 

GO TO 95 

92 QO = 2./ (GAMA*PSA*MACHA*MACHA) 

C BEGIN LOOP THROUGH STATIONS 
95 CHOKE = .FALSE. 


- 132 - 



IF I lL J= J 
JSUP = ■; 

L I M L > - 
L I N I B = 0 
L - L 0 

VJj ML B =0. 

pub = . 

WF = 0 . 

C SUPSGN lC/SUPhKSOM IC BRANCH SELtCTION 

M = KLB(L) 

CALL j t T IX 
JA = J 
P A A = P 
M = M U B ( L ) 

CALL UETlX 
Jl> = J 

M i?. B = M 

I F ( JSUV .1 0.0 ) subson = .truf: . 

IF(SSLF) SUbSON=. FALSE. 

IF ( SCHOKE (L ) .NE.XCHOKE > GO TO 510 
IF ( SS OF ) SUB SUN = .FALSE. 

JSUK = JA+256*JB 

C EXECUTE FLOW BALANCE 

5 i J CALL FLOBAL 

IF( TVPELB(L) .EO.TE .OR. T YPE UB ( L) . EQ. TE ) JSUM=0 

C BRANCH AND ASTERP ARE PRINTOUT INDICATORS 

DATA DBSTAR/2H**/, SUB/3HSUB/, SUPER/5HSUPER/ ♦ I CHOKE /5HCH0KE / 
5 J i A S T E < P = BLANK 

I F ( PR I M ( L ) ) ASTERP=DBSTAR 

B R A N C H = SUPER 

IF(SUBSON) BR ANCH= SUB 

IF ( SCHOKE ( L ) .EG.XCHOKE ) B RANC H* I CHOKE 

CALL SfcTM( 1, BLANK, CHANLS,10) 

CALL VOVE ( 2 ,ZF(MA ) ,Z,NK, I , RF ( M A ) » R , NK , 1 ) 

CALL MOVF ( 2 ♦ CURVE ( MA ) »CURV,NK» 1 * VMF < MA ) , VM , NK , l ) 

L 0 =0 

K = i 

P = R A 

5^0 F L 0 w ( K ) = W S T A { K ) *C G 

PHI ( K ) = PHI1(M)*T0DEG 
OGAM = FGR(K )/ ( i .+FGRIK ) ) 

MACH<K)=VP(K)*SURT(QGAM/(RG{K)*TS(K> )) 

A 0 A R L F ( K ) = R ( K ) 

IF ( AXIA ) AQ ARE F ( K > = PI*R(K)*R(K) 

PS(K> = RHO(K)*RG(K)*TS(K) 

PS0PO(K)=PS(K)/PSA 
PSOPT { K ) = PS ( K ) / P T IK) 

TSOH ( K > = r S ( K } / T T ( K ) 

C CP MUST FOLLOW USE OF RG 

C P ( K ) = ( P S ( K ) -P SA ) *00 
CALL GETIX 
X I c ( K ) = X 2 ( J ) 

IF ( SLCHN( J) .LQ.CHANLS(LQ) ) GO TO 530 
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LG = LQ +1 


J 1 ( LO ) = J 
K 1 ( LU ) = K 

CHANL S ( LQ )=SLCHN ( J ) 

IF(LO.GT.l) FLOWMX(LQ-l)=FLOW(K) 

I = 0 

525 I « 1+1 

IF ( SLCHNU ) .NE. ICHN( I ) . AND. I . LT . I C ) GO TO 525 
QPT 0 = 1 . / WP TO ( I ) 

5 3 D PTQPTO(K)=PT(K)*QPTO 


K = K+l 

M = M + l 

IF(K.LE.NK) GO TO 520 
J 1 ( LQ + 1 ) = J + 1 
K 1 ( LQ+1 >=K 

FL0WMX(LG)=FL0W(K-1) 

LQS = 0 

533 LQS = LQS+1 

KB = K 1 (LQS) 

KE = K 1 ( LQS + 1 )-l 

FL MX = 1 ./FLOWMX ( LQS) 
DO 5 3 5 K = K B » K E 
5 35 PFLOW(K ) =FLOW ( K ) *FLMX 

IF( LQS.LT.LQ) GO TO 533 


XII = X1(L) 

IF( PRPRN.EQ. ( -1 ) ) GO TO 610 
CALL FHEAD(LINEA+NK> 

LINEA = A 

I F ( .NOT.PRIM(L) ) L INEA=8 

WRITE ( 6t 1600 ) XI 1 , ASTERP .CHANLS, BRANCH, 

* ( X lc (K ) ,PFLOW(K) t Z (K ) |R ( K I t PH I (K ) .CURVI K > , PSQPO ( K ) .PSOPT(K) , 

2 TSQTT(K ) »CP< K ) , M ACH ( K ) f AQAREF ( K ) ? PTQPTO(K) ,K=l,NK) 

1600 FORMAT (/25H STATION COORDINATE. X 1 1 = » F7 . 3 . A2 , 13 H CHANNELS- . 
i 10 ( A6.2X ),A5// 5X» 1 3HX I 2 STRM FNCT .6X, 3HX , Z , 8X , 3HY , R , 8X , 3HPH I , 

16X,AHCURV»6X,21HPS/P0 PS/PT TS/TT.6X.2HCP.6X. A-HMACH.6X » 

3 6H AREA*3X»6HPT/PT0 / ( 2 X ,F 6. 3 , F 10 . 3 , F 1 2 . 5 ,F 1 1 . 5 , F9 .3 , F 1 1 . 5 , 
A F9.3»2F8.3»F10.3»F9.A.F 1 l . 3 » F9 . 3 . 7X , ) . ) 


610 IF( .NOT .PRIM(L) ) GO TO 800 
M = M A 

DO 620 K = 1 » NK 
COS PH I = COS ( PH I 1 ( M ) ) 

S I NPH I = SIN(PHIKM)) 

F VX ( K ) =VM( K ) ♦COSPH I 
F V Y ( K ) =VM ( K ) *S INPHI 
FPX ( K )=(PS(K)-PSA)*COSPHI 
FPY (K ) = <PS(K)-PSA )*S INPHI 
620 M = M+i 

SV X ( 1 ) = 0. 

S V Y ( 1 )= 0. 

spx m= o. 

SPY ( 1 )= 0. 

CALL LSPFIT( WSTA.FVX.NK, WSTA.SVX.NK, -1) 
CALL LSPF IT ( WSTA.FVY »NK» WSTA.SVY.NK, -1) 
CALL LSPFIT (AREA»FPX f NK . AREA »SPX» NK . -1) 
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FALL L SPF I T I ARE A, FPY ,NK, AREA, SPY, NK, -1) 

Of .' 6 i J K = 1 , N K 
S T X I K ) = SVX (K ) + ‘>PX (K ) 

6 iu ST Y ( K ) = SVY ( K ) +SP Y (K ) 

K A = i 

ntJ 64 j L L = 1 » L 0 

J = J 1 ( LL + 1 )- 1 

K = K1CLL+1I-I 

IM MH.NE.O) GO TO 63 5 
ST XU ( J )=STX(K)-STX(KA) 

S T YU ( J)=STY(K)-STY(KA) 

635 IF(MU.NE.O) GO TO 640 
STYDI J )=STY(K)-STY(KA) 

STXDl J )=STX(K)-STX(KA) 

640 KA = K 

IF ( PKPRN.EQ. ( -1 ) ) GO TO 800 

WKlTh (6,1700) SVX(NK) ,SVY(NK),SPX(NK) ,SPY(NK) ,STX(NK) ,STY(NK) 

LINES = LINES+4 

17 JO FORMAT ( /6X25HSUM~VM*COS< PHI ) *OFLOW =F10. 2 , 36X, 25hSUM-VM*S I NI PHI ) 
**DFLUW =F10.2,/6X25HSUM-(P-PS0)*C0S(PHI )*DA =F 1 0 . 2 , 36X , 25HSUM- ( P 
*-PSU) *S IN ( PHI )*0A =F10.2 , /6X45HT0T AXIAL MOMENTUM FLUX =F10.2,36X, 
♦25HTUTAL Y-MOMENTUM FLUX =F10.2t) 

C RELOCATE DATA INTO THE M-ARRA Y A 

8.iO CALL MOVE ( 2 , MACH , MACHM ( MA ) , NK, 1 , PS ,PSM{ MA ) , NK, 1 ) 

CALL MOVE! 2»PT»PTM{MA) ,NK ,1, TT , TTM ( MA ) , NK , 1 ) 

C FILL IN STAGNATION POINT VALUfcS 
IF(MLBIL).EQ.MA) GO TO 820 
M = M L B ( L ) 

CALL GETIX 
MACHM ( V ) =0 • 

PTM( v, ) - PTM ( MU ) 

PSM(k ) = P T M ( M ) 

TTM ( M ) = TTM( MU ) 

V M F ( V. ) = 0 . 

8 2‘J IF ( *Utt(L ) .EG.MB ) GO TO 830 
M = MUBIL) 

CALL GETIX 
MACHM ( y J =0. 

PTMIM ) = P T M ( M U ) 

PSM ( M )=PTM(M) 

TTM(V)=TTM(MU) 

VMF(M)= ). 

C INDEX 10 NEXT STATION 

830 L = L+LNEXT (L ) 

IF ( L.LT .LESTA) GO TO 500 

RE TURN 
END 
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-WRIBDY- 


*DECK WR I M) Y 

SUBROUTINE WR I BO Y 

*WR 1 BOY WRITE OUTPUT FOR EACH BOUNOARY 


COMMON /BCOMMN/ 
LOGICAL 

COMMON /ADAM01/ 
COMMON /CLINES/ 
COMMON / CNORM 
COMMON /ERASE2/ 

* 

* PSQPTI1 

NEW VARIABLES FOR 
COMMON /ERASE 3/ 
REAL 

DIMENSION 
EQUIVALENCE 
COMMON /CFRFLD/ 

COMMON /ALLCOM/ 


REAL 
LOGICAL 
LOGICAL 
BOUNDARY TABLE 
INDEX- LB=LBDO»LBDE 
LBNEX T= INCREMENT TO 
INCREMENT TO 
CHANNEL WITH 


PR0GMI9) , FILIN, FILOT 
FILIN, FILQT 

NAME < 6 ) , ADDKcSI 6) , TI TLE (6) , IDENT (6) 

LINES, 0M1TFK.PTITLEI6) 

/ RHL ,RM, AHL, ARM 

X 1 1 ( 100) ,SW(iOO)«ZW( 100 )»RW(100) , ANGW ( 100 ) , 

CURVWI 100 ) ,VE 1 100 J , MACH II 00) ,PSQPO( 100) ,CP( 100 ) , 
OJ) ,PTQPTO(10j), TT(IOO) ,AW( 100 ) , SPDA ( 100 ) 

NASA VERSION ONLY — PSQPT AND PTQPT0 
AQAN( 100) ,CDP I ( 100 ) , PSMPOI 100 ) 

MACH 

XWU) , YW ( 1) 

( XW, ZW) , ( YW, RW) 

FSAVI 300 ) , STXU< 128) ,STXDI 128) ,STYU( 128) ,STYD( 128) 

MACHA,PSA,TSA,PTA, TTA, AXI A , RGA, GAMA, 

MACHC.PSC ,TSC ,PTC, TTC, AX I C ,RGC, GAMC , 

DAXIT,SCALEA,TTE,CHOTST 

M ACHA ( 1 ) , MACHC 

AXIA.AXIC 

C HOT ST 


POINT (=0 BEFORE COALLATIO 
DATA IS ASSOCIATED 


UPPER SURFACE 


NEXT BOUNDARY 

LBZ1 = INCREMENT TO THE FIRST BOUNDARY 
CHN AM E= CHANNEL WITH WHICH THE BOUNDARY 
UP = T OR F FOR UPPER OR LOWER BOUNDARY 
LEDEX = RELATIVE INDEX OF L . E * POINT WHEN LOWER AND 
CONTOURS ARE CONNECTED 

BDNAMt ,LBA,LBB=NAME AND INDEX LIMITS OF SPECIFIC BOUNDARY 

DATA WHEN BOUNDARIES ARE CDALLATED 
DIMENSION BDT(l) , L BNE X T U),LBZ1(1), 

1 CHNAME (1 ) »UP( 1) , LEDEX ( 1 ) , 

2 ZBT( 1) »RBT< 1) , A NG B T ( 4 2 ) 

LOGICAL UP 

INTEGER BDT, CHNAME, BDNAME 

DIMENSION BDNAME ( 1 ) , LB A < 1 ) , LBB ( 1 ) 

EQUIVALENCE ( BDNAME , ZBT ) , (LBA.RBT), ( LBB , ANGBT ) 

STATION TABLE 
INDEX- L=LO, LESTA 

SCHOK.E= STATION CHOKE INDICATOR ( AD J WF , BRHS , WR IOUT ) 

MCL = SHARP CORNER INDICATOR ( BLDTBS ) 

MCL = FIELD INDEX OF CONTROL STREAMLINE ( PTMOVE , FLOBAL ) 

COMMON /CHDATA/ X 1 ( 1 ) , LN E XT ( 1 ) , ML B ( 1 ) , MUB C l ) * PRI M (1) , 

1 T YPELB ( i ) .NAMELB ( 1 ) , I LB ( 1 ) , FLB ( l ) , S 1LB( 1 ) , 

1 TYPEUB(l),NAMEUBU),IUB(l) , FUB ( 1 ) , S 1UB { l ) , 

3 VMB ( 1 ) »D WDV t 1 ) , X2CL(1),VCL(1),MCL(481) 

LOGICAL PRIM 

INTEGER TYPEL8 , TYPEUB 


DIMENSION 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


SCHOKEl 1 ) 

( SCHOKE , D WDV ) 

( XI, BDT) , (LNEXT,LBNEXT) , (MLB, LBZl) , (MUB, CHNAME) 
( PRIM, UP ) , (TYPE LB, LEDEX) , (NAMELB, ZBT) , (ILB,RBT) 
( FLB, ANGBT) 
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c 


COMMON / C F B 


* 


LUGI CAL 


L , MA , MB » PLB » PUB » WF , C HOKE , SUB SON, NK , PLBC , PUBC , 
XCHOKfc, T ARE A , VMBC , WRQST,WCALC, OV I H ) , QV P ( 8 ) , 
J SUM ,VMLBSQ 

CHOKE , SUBSON 


INULX- 

M- Ml.) , NM 


COMMfjN 

/U 

/ 

Common 

/CR 

/ 

COMMON 

/CS1 

/ 

COMMON 

/CPHI i 

/ 

COMMON 

/CM 

/ 

COMMON 

/CCURV 

/ 

COMMON 

/C I DEX 

/ 

C UMMuN 

/C VM 

/ 

COMMON 

/CDS2 

/ 

REAL 



COMMON 

/CPSM 

/ 

COMMON 

/CS2 

/ 

COMMON 

/CRHS 

/ 

COMMON 

/ I XUR IG/ 


* 

* 


D I MEN 

SION 


EQUIVALENCE 


common 

/CBEND 

/ 

COMMON 

/CBITS 

/ 

COMMO j 

/CCUBE 

/ 

COMMON 

/CGRAV 

/ 

common 

/CPI 

/ 

COMMO m 

/CREF IN/ 

, 

COMMON 

/SLTAB 

/ 

I NT L G 

er slch; 


common 

/CHNFPT/ 

COMMON 

/BLBDY 

/ 


L ( 30U) 

R ( 300 ) 

Si (300) 

PH II ( 300) 

J M S ( 300) 

C URV ( 300 ) 

M , J, MU,MD, I STAG 
VM ( 300 ) 

MACHM ( 30 j ) 

MACHM 
PSM( 300) 

P T M ( 300) 

T TM( 300) 

L HO, LHE , LBOU ,LBDE , LTO,LTE, LWO,LWE, LFO»LFC 
LO,LLSTA, LSG,LSE, LDO , LDE , LDUM(A), 

MO , NM , N J , NFCOL S ♦ MAXN J , MAXOL , MA XNM , MAXLE , 
LEO, LEE, LRO , LRE , IRD 
LIMITS! 24) 

( LIMITS, LHO) 

NBCB ( 2 ) , ANGE (2) ,CURVE( 2) , F B C 2 ) 

BITS , BLANK 

N6C( 2) ,C I ( 2 ) , C2 ( 2 ) , FEND (2) 

CG 

P I , TWOPI ,PIQ2»PIQ4 , TODEG , TOR AD 
SG1,SG2,VMGI ,VMG2 

NGR , NGZ , SGR( 10) ,GR( 10) , SGZ ( 10) , GZ ( 10 ) 

W( 128) ,X2( 128 ) , SLCHN ( 128 ) 

ICHN ( 10) , WTFS( 10) , WTFA( 10) , WPTO( 1 0 ) , WTTO ( 10 ) , 
bLB( 60) 


DIMENSION IBLB(60) 


IC 


EQUIVALENCE (IBLB.BLB) 

COMMON /BLOT A / BNAME , LOWER, I BTYPE ,N1 ,NI ,CAPX i 
INTEGER BNAME 


INTEGER HLL.HTE, ASL,BDY,TSL,CHNN,CHN,XK5SV,XKEYB,BLANK 

LOGICAL DOUBLE, LOWER, UPPER 

DIMENSION L 0 WUP ( 2 ) 

DATA L0WUP/5HL0WER ,5HUPPER/ 

DATA HLE,HTE/2HLE,2HTE/» ASL , BDY, TSL/3HASL ,3HBDY ,3HTSL/ 
NT RY = i 

C DEFINE REFERENCE DYNAMIC PRESSURE, ETC 
00 = 0 . 

IF ( MACHA.LE. . 1 ) GO TO 95 

IF(GAMA.NE.U. ) GO TO 92 

CO = ( RGA*TSA)/ (PSA*MACHA*MACHA ) 

GO TO 9 5 

Nc CO = 2./(GAMA*PSA*MACHA*MACHA) 
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C BEGIN LOOP THROUGH CHANNELS 

95 LINES = 64 
IUP =4 
NCHN = 1 
J 2 = 1 

105 CHNN = SLCHN ( J 2 ) 

LOWER = .TRUE. 

I = 0 

I j7 I = 1+1 

IF ( CHNN.NE. ICHNI I ) .AND. I.LT.ICJ GO TO 107 
QPTO = l./WPTO(I) 

QTTU = l./WTTOII) 

GO TO 122 
110 J2 = J 2+ 1 

IFIJ2.EQ.NJ .OR. SLCHNl J2 + 1) .NE.CHNN) GO TO 120 
GO TO 110 

120 LOWER = .FALSE. 

C BUILD I -SUBSCR I P TED ARRAYS 

122 M = MBEG I N ( J2 ) 

L =0 

SPDAS V= 0. 

XK5SV = BDY 

123 I =1 

SWORG = S1(M) 

PTO = PTM(M) 

TTO = TTM(M) 

TTQT T 0= TTM ( M ) *QTTO 

124 DOUBL E= .FALSE. 

125 SW( I) = S1(M) - SWORG 

7 W ( 1 ) = Z(M) 

RW ( I ) = R(M) 

ANGW ( l ) =P H I l ( M ) *TODEG 
CURVW ( I )= CURV(M) 

PS = PSM(M) 

PSQPT ( l ) = PS/PTM(M } 

PTQPTO ( I )=PTM(M)*QPTO 
MACH( 1 ) =M ACHM ( M ) 

VE ( I ) — V M ( M ) 

AW ( I ) = RW ( I ) 

I F ( AXIA ) AWI I >=PI*RW(I)*RW< I ) 

PSQPO ( I ) =PS/P SA 

PSMPQI I )=PS-PSA 

CP ( I ) = P SMPO ( I )*Q0 

I F ( LOWER ) PSMPO(I)=-PSMPOm 

CALL GET I X 

CALL STANO(M,L, UPPER) 

X I 1 ( I ) = X 1 ( L ) 

N I = I 

I = 1 + 1 

IF(NI.EO.l) GO TO 160 

C CHECK FOR LEADING EDGE POINT 

IF( ISTAG.NE.l) GO TO 140 

IF ( TYPELBI L ) .EO .HLE .OR. TYPtUB IL ) .EQ.HLE ) GO TO 170 
C I S TAG= 1 

IF(DOUBLE) GO TO 160 
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muji'i i = .Tkut . 

GO !U 1 2‘> 

C CHLf.r. F0° TRAILING EDGE POINT 

i A j IP ( I'j T AG. N L . 2 ) GO TO 160 
C I S T A , = 2 

IP ( I Y k L L ti ( L ) .LG.HTE .OR. T YPt UB ( L ) . EQ . HTE ) GO TU 190 

I S T A G = G i 3 OR DOUBLE = T 
M = v 0 

IHM.J.-I GO TO 124 

G 0 TO 1 « 0 

C APPROACH S T R E A M L I N L 

1 i G X K L Y i'j = A S L 

GO TO 200 
C BOG Y S UR P AC E 

in:- X K t Y G = X K A S V 

GO Tl' 200 

C 1 railing STREAMLINE 
190 X < L Y '> - X K 6 S V 

X< b i>V = TSL 

2 u j IP { X R t Y B .to. TSL) GO TO 2 2<J 

i P ( • N U T . L 0 W E R ) GO TO 2 20 
LB = LPF ( NAMELB (L ) ) 

IF tLtOEX(LB).EO.O) GO TO 220 
C LOOP TO FIND BOUNDARY NAME OF UPPER SIDE OF L.E. 

L B X = LB 

2iA IP ( LKAl LBX ) .GE.LEDEX(LB) ) GO TO 220 
L 0 X = LBX + 3 

IF (Li X.LI .(LB + LBZ1ILB) )) GO TO 214 

CALL lRRORI 
2 w u S P D A ( * ) = SP D A S V 

CALL L SUM ( AW.PSMPUfNIt SPDA) 

S P D A S Y = S P D A ( N I ) 

ARK = RM 

IF ( AXIA ) ARM = PI*RM*RM 
DO cc'.'i I = 1 » N I 
AW ( I ) = I ARM- AW ( 1 ) ) / AR M 
2_b C0Pi:i) = SPOAII) *QO/ARM 
ADDG = SPDASV*QO/ARM 
2 iJ LINbS = 64 

CALL PHEADINI+6) 

KUP = 2 
IFtLJrtERJ K UP = 1 
CHN = SLCHN ( J 2 ) 

X 1 2 = X 2 ( J 2 ) 

SWORu = 0. 

WRITE (6,12)0) LOWUP(KUP ) ,CHN,XI2, ( XI l (I) « SW (I) * Z W (I) • RW ( I ) • 

* ANGW ( I ) ,CURVW( I ) ,PSOPO( I ) ,CP ( I ) , PSQPT ( I ) * MACH { I ) ,CDP1 ( I ) , AW ( I ) , 

* P T 0 P T 0 ( I ) , 1 = 1, NI ) 

.2^0 FORMAT </2X,A6,17H BOUNDARY TU CHN=,A6,31H, STREAMLINE COORDINAT 
*L, XI2=»F7.3, 1H.// SX, 3HXI i ,6 X , 3HS 1 W , 7X , 5HX W , Z W ,6 X , 5H YW , R W , 5X , 

* 4 HAN GW, 5X, BHCURVW,5X, 5H P S/PO , 5 X , 2HC P t 4X , 5H P S / PT , 4X , 4HMACH , 5X , 

* MHC0PI,14H ( AMAX-A) /AMAX f 8H PT/PTO / I 2 X 1 2 F8. 3 , F 1 2 . 5 , F 1 1 . 5 , 

* PB. P,F11.6,2F9.3,F7.3»2F9.4fF14.3*F8.3*)» ) 
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WRITE (6,1210) TTQTTO 
1210 FORMAT ( / 6X , bHT T / T TO =,F9.3) 

IF ( XKEYB.LQ.ASL ) WRITE (6,1220) ADOG 
1220 FORMAT ( / 6X , 1 SHADE) I T I V£ DRAG = ,F 9.4) 

I F ( XK.fcYB.6Q.ASL .OR. XKfcYB.tQ.TSL ) GO TO 309 
£*#**$*$♦«*******$**:*#*:***#*** BOUNDARY LAYER *******♦******♦♦*♦#**♦♦$ 


NAMt = NAMfcLB(L) 

I F ( . NUT. LOWER ) N AME =NAMEUB ( L ) 

LBL = L BDYBL ( NAME, LOWER ) 

I F ( LBL. tQ. 0 ) GO TO 309 
CAPX i = BL B ( L BL + 2 ) 

BN AML = I BL R ( L 8L ) 

CALL SAB(NTRY) 

NTRY = 2 

309 I F ( MD.GT.O ) GO TO 123 


C INTEGRAL MOMENTUM BALANCE ON THE CHANNEL 
1 F ( .NOT. LOWER) GO TO 310 
PFLB = SPDASV 
GO TO 110 

310 PFUB = SPDASV 

F TO T = STXU( J2 ) +PFLB + PFUB 
FfcRR = F TO T- ST XD ( J2 ) 

WRITE (6,1300) CHN,STXU( J 2 ) , PFLB , PFUB , F TOT , ST XD( J2) ,FERR 
1300 FORMAT ( / IX 32H INTEGRAL MOMENTUM BALANCE, CHN = A6, 2 X 19H ( AX I AL FORCES 

* ONLY 1/6X31 H ENTERING MOMENTUM -F11.^,/6X31HL0WER BOUND 

*ARY PRESSURE FORCE = F 1 1 . 4 , /6X 31 HUPPE R BOUNDARY PRESSURE FORCE = F11 
* • A , / 1 2X 1 2HSUM OF A BO VE F 2 4 . A , / 6X3 1HLE AV I NG MOMENTUM = F 

# 1 1 . 4 , /12X26HERROR =F11.4,) 

J 2 = J 2 ♦ 1 

IFtJ2.LE.NJ) GO TU 106 

I F ( NTRY.EQ.2 ) CALL SAB ( 3 ) 

RETURN 

END 


-l40- 



OOOOOOOOO o o oo 


♦ DECK FRROkY 

SUBKOuT INE ERROR! 

ClOoMPX FDuMP FJK STC EXLCUTF SECTION -EDUMPX- 


C 1 


ALL COM 


C f 1 M M U N / A 1. 1 CO V/ 


kL AL 
L G G I C A L 
LOGICAL 

COMMON /CFB / 


LOGICAL 

COMMON /ERASE / 

common /erase2/ 


REAL 

01 Ml NS ION 

equivalence 

D l MENS ION 
EQUIVALENCE 
HELD TABLES 


MAGMA, PSA ,TSA ,PTA, Tl A, AX I A ,KGA, GAMA , 

MACHC i PSC , TSC ,PTC ,TTC , AX I C » RGC , GAMC • 

DA X I T , SC ALE A , TTE .CHQTST 
MACHAI I ) .MACHC 
AXIA.AXIC 
CH0TST 

L .MA.MB. PLB.PUB .WF . CHOKE . SUB SON , NK , PLBC , PUBC , 
XCHOKE, TAREA.VMBC, WRQST.WCALC. QV ( 8 ) . OVP ( 8 ) , 
JSUM, VML8SQ 

CHOKE .SUBSON 

ERASEC(BOO) 

area (96) , ARE AO ( 96) .DISP(96) , P T ( 96 ) , L AMBQA ( 9 6 ) . 
R HO (96)»SQRTVV(96) »TS(96) »TT(96) »VMSQ(96) . 

V VKQ KP ( 96 ) » 

WQA{ 96 ) , WSTA (96 ) , RG ( 96 ) ,C2CP ( 96 ) , FGR ( 96 ) 

L AMBDA 

L S2 ( 96 ) » SDN(JRM(96) 

(ES2.VVKQKP) . (SDNQRM.RHU) 
kCU( 96 ) 

( RCU.LAMBDA ) 


INDEX- M = MO , NM 
COMMON /CL t Z ( 300 ) 

COMMON /CR / R ( 300 ) 

COMMON /CS2 / S2 ( 300 ) 

COMMON /CSl / S i ( 300 ) 

COMMON / C P H I 1 / PH 1 1 ( 300 ) 

COMMON /CM / JMS( 300 ) 

COMMON /CCURV / CURV(300) 

COMMON /CR / rf ( 300 ) 

COMMON /CIDEX / M. J.MU.MD.ISTAG 
TABLE OF INDEX LIMITS 

COM Mi j N /IXORIG/ LHO.LHE, LBDO.LBOE, LTO.LTE, LWO.LWE, LFO.LFE. 

* LO,LtSTA,LSO,LSE,LDO,LDE,LDUM(A) , 

* MO.NM, NJ.NFCOLS, MAXNJ .MAXOL . MAXNM , MAXLE . 

* LEO. LEE, LRO.LRE.LRD 

DIMlNSICN LIMITS(24) 

EQUIVALENCE (LIMITS, LHO) 

COMMUN /CVM / VM ( 300 ) 


STREAMLINE TABLE 

COMMON /SLTAB / W ( 12 8 ) , X 2 ( 12 8 ) , SLCHN ( 1 28 ) 

INILGLR SLCHN 
BOUNDARY TABLE 
INDEX- LB=LRDO,LBDE 

L R N E X T = INCREMENT TO NEXT BOUNDARY 

LUZ1 = INCREMENT TO THE FIRST BOUNDARY POINT (=j BEFORE COALLATIO 
CHN AM fc = CHANNEL WITH WHICH THE BOUNDARY DATA IS ASSOCIATED 
UP = T OR F FOR UPPER OR LOWER BOUNDARY 

LEDEX = RELATIVE INDEX OF L.t. POINT WHEN LOWER AND UPPER SURFACE 
CONTOURS ARE CONNECTED 

BDNAMfc ,LBA,LBB=NAME AND INDEX LIMITS OF SPECIFIC BOUNDARY 
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DATA WHEN BOUNDARIES ARE COALLATED 
DIMLNSION BDT( I) ♦LBNEXTd ) • LB Z 1 ( 1 > • 

1 CHNAME ( 1 ) ,UP ( 1) fLEOEX(l) • 

2 ZBT(l) ,RBT(1) ,ANGBT<42) 

LOGICAL UP 

INTEGER BDT, CHNAME, BDNAME 
D I ME NS I ON BDNAME ( 1 ) ,LB A ( 1 ) , LBB ( I ) 

EQUIVALENCE (BDNAME, ZBT), (LBA.RBT), ( LBB , ANGBT ) 

FLOW ADJUSTMENT TABLE 
INDEX- LF=LFO»LFE 
NFCOL S= 8 

X i F = ORTHOGONAL COORDINATE 

X2F = STREAMLINE COORDINATE OF SL EMINATING FROM T.E. 

XI8F - X 1-COORDINATE OF CHOKE STATION OF FLOW BELOW T.E. 

X1AF = X 1 -COORD I NATE OF CHOKE STATION OF FLOW ABOVE T.E. 

S1F = SI-COORDINATE OF T.E. (UPPER SURFACE). THIS ITEM 
IS USED WHEN INTERPOLATING FOR WAKE DELTA-STAR. 

LFB.LF A= INOICES OF STATIONS BELOW AND ABOVE T.E. 

NCHB, NCHA=NUMBER OF CHANNELS BELOW AND ABOVE T.E. 

LRF = INDEX OF DUMMY ORTCHN LIST FOR THE T.E. 

LRXF = INDEX OF LAST CHANNEL BELOW THE T.E. 

JORDE R= 0 IF TOTAL FLOW AT X1F IS GIVEN 
= 2 IF FLOW ABOVE T.E. IS GIVEN 
= 1 IF FLOW BELOW T.E. IS GIVEN 
JORDE R = -I IF FLOW AT X1F IS CHOKED AND SINGLE CHANNEL 
DIMENSION X1F( 1 ) ,X2F( I ) ,X1BF ( 1 ) ,X1AF U ) , 

1 S I F ( I ) ,NCH8 ( 1 > , NCHA ( 1 ) , JORDER ( I ) ,VNR(I2) 

EQUIVALENCE ( LFB.XIBF ) , ( LFA , XI AF ) , (LRF .NCHB) , (LRXF. NCHA) 
DIMENSION LFB( 1) ,LFA(I) f LRFCl) »LRXF ( 1 ) 

STATION TABLE 
INDEX- L=LO,LESTA 

SCHOK L = STATION CHOKE INDICATOR I ADJWF .BRHS.WRIOUT ) 

MCL = SHARP CORNER INDICATOR (BLDTBS) 

MCL = FIELD INDEX OF CONTROL STREAMLINE ( PT MOV E , F LOBAL ) 

COMMON /CHDATA/ X 1 ( i ) , LNE XT ( 1 ) , ML B ( 1 ) , MUB ( 1 ) , PRI M ( 1 ) « 

1 T YPELB ( 1 ) ,NAMELB( 1 ) , ILB(l) , FLB ( 1 ) , S I LB ( 1 ) , 

1 T YPEUB ( 1 ) .NAMEUB ( l).IUB(l) ,FUB(I ) ,SIUB( l ) , 

3 VMB ( 1 ) » D WDV ( i ) » X2CL ( 1 ) , VCL ( 1 ) »MCL ( 481 ) 

LOGICAL PRIM 

DIMENSION SCHOKE(l) 

EQUIVALENCE ( SCHOKE , DWDV ) 

EQUIVALENCE ( BDT , X1F , XI ) , ( L BNE XT , X2F , LNE XT ) , ( LBZ 1 , XI BF , MLB ) 

EQUIVALENCE ( CHNAME , X 1 AF , MUB ) , ( UP , S1F , PR I M ) 

EQUIVALENCE ( L E DE X , NC HB , T YPE LB ) , ( ZBT, NCHA ,NAMELB ) 

EQUIVALENCE ( RBT , JORDER, I LB ) , ( ANGBT , VNR ,F LB ) 

COMMON /CTABPR/ 11TAB 
COMMON /BLBDY / IBLB160) 

CALL TABPRT ( 6H ALL COM , MAC HA , 2 0 , 8 ) 

CALL TABPRT ( 3HCFB»L,33,4) 

CALL TABPRT(8HCID£X,M,5,5) 

IITAB = LBDO 

CALL TABPRT (6HBDYTAB, BDT , LBDE , 3 ) 

IITAB = LFO 

CALL TABPRT(6HCADJWF,X1F,LFE,8) 

IITAB = LO 
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CALL I AHPR I ( 0HS1 AT AH ,X1 , Lt ST A ,b ) 

1-,'J WRIT, (A,11S0) I J , X2 ( J ) , SLCHNt J ) , W ( J ) , J= A ,NJl 
(.ALL l ABPRT ( SH! RASE, ERAS LC.HK, b) 


CALL JMSPRI 

CALL l ABPR I ( 2HS 1 , Si ,NM, 1 ■> ) 

CALL TABPRT ( 2 H S 2 » S 2 , NM , 1 9 ) 

CALL I ABPRT ( 1HZ , / t NM , 10 ) 

CALL TABPRT ( 1HR,R,NM, 10) 

CALL T ARPR T ( AHPHI L , PHI 1, N«» 1J) 

CALL TABPRT (AHCURV,CURV,NM, 10) 

CALL TABPRT( 2HVM,VM,NM,10) 

CALL TABPRT ( 1HB»B»NM,10) 

CALL TABPRT ( 6HER AS E2 , ARE A , 15 36, 8 ) 

1 F ( JDL0lll.ME.-3 ) CALL T ABPR T < 5HB LBDY , I BLB , 60 , 3 ) 
I F ( LOL.EO.j ) GO TO 132 1 
UTAH = LOO 

CALL TABPRT ( bHBL TAB, CHNAM,LDE, 3) 

1321 CONTINUE 


LSTOP = b 

GO TO (999,999) , LSTOP 

999 RETURN 

1 1 SO FOR RAT ( ///1X17H STREAMLINE TABLE -/ 1 7X 32HJ X2 

* */( I lb,F12.6,6X,A6,F12.6, ), ) 

EN U 


SLCHN 
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♦ DECK BDYPTM 

SUMRIIUT I N f BDYPTM (NAME , I NTVL , Z D, RD . FD , Si DD , DS L, DS1 GMA ) 

* HI) YP I M BOUNDARY POINT MOVE MEN T -BDYPTM- 

INPU I - 

BUI = BOUNDARY I ABLE 

NAMe; = BOUNDARY NAME 

INTVL = INDEX OF INTERVAL OF THE INPUT POINT IN THE BOUNDARY TABLE 

FO = FRACTION POSITION OF THE INPUT POINT IN THE INTERVAL 

S 1 DD = ARC DISTANCE FROM THE BEGINING OF THE INPUT INTERVAL 

DS I = REQ-D MOVEMENT IN THE CLOCKWISE DIRECTION FROM THE INPUT P 

OUTPUT- 

INTVL = INDEX OF INTERVAL OF THE OUTPUT POINT 
Z D» RD = COORDINATES OF THE CALCULATED OUTPUT POINT 
ANGD = ANGLE OF OUTPUT POINT 
CJRVD = CURVATURE OF OUTPUT POINT 
FD = FRACTION POSITION IN THE OUTPUT INTERVAL 
SIDD = ARC DISTANCE FROM THE BEGINING OF THE OUTPUT INTERVAL 
DS1GMA= -GET- MINUS -ASK- POINT MOVEMENT DISTANCE 

BOUNDARY TABLE 
INDEX- LB=LBOO,LBDE 

L B N E X T = INCREMENT TO NEXT BOUNDARY 

LBZ1 = INCREMENT TO THE FIRST BOUNDARY POINT (=0 BEFORE COALLATIO 
C HN AME = CHANNEL WITH WHICH THE BOUNDARY DATA IS ASSOCIATED 
UP * T OR F FOR UPPER OR LOWER BOUNDARY 

LEDEx = RELATIVE INDEX OF L.E. POINT WHEN LOWER AND UPPER SURFACE 
CONTOURS ARE CONNECTED 

BDNAME t LBA,LBB=NAME AND INDEX LIMITS OF SPECIFIC BOUNDARY 

DATA WHEN BOUNDARIES ARE COALLATED 
COMMON /CHDATA/ BDTU ) , L BNE X T <l) , LBZ 1 ( l > , 

: CHNAME (l),UP(l),LEDEX(l), 

2 ZBTIDtRBKII , ANGBT ( 42 ) 

LOGICAL UP 

INTEGER BDT , CHNAME , BDNAME 
DIMENSION BDNAME (l I »LBA ( 1 ) » LBB ( l ) 

EQUIVALENCE ( BDNAME , ZBT ) , ( LBA ,RBT ) * ( LBB » ANGBT I 
COMMON /CBEAM2/ DR , DZ , YP A , YPB ,F ,G • DX, YQDX, ZM,RM, ANGM,CURVM, SIM, 
i RZONLY, ANGCHD, SINTVL , YPASO, YPAB* YPBSQ 

LOGICAL RZONLY 

COMMON /IXORIG/ LHO, LHE , LBDO.LBDE, LTO.LTE, LWO,LWE, LFO,LFE, 

* LO,LESTA, LSO.LSE, LOO.LDE, L0UM<4), 

* MO » NM , N JtNFCOLS, MAXN J ,MAXOL ,MAXNM, MAXLE , 

* LEO, LEE, LRO»LRE»LRD 

DIMENSION LIMITSI24) 

EQUIVALENCE (LIMITS, LHO) 

COMMON /CBDYP T / ANGD » CUR VD 
COMMON /CBITS / BITS, BLANK 

COMMON /CFB / L,DUMCFB(33) 

COMMON /CPRINT/ PPDUM( 6) ,PDUM(6 ) 

COMMON /BLBDY / BLB(60> 

DIMENSION I BL B ( 60 ) 

EQUIVALENCE ( I BLB , BL B ) 

COMMON /REBL / RESTBL 
LOGICAL RESTBL 

COMMON /CPI / PI, DUMPI ( 5 ) 
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0 1 ME NS I UN 

1 4 T L i k 

i OU IVALLNCL 


LOGICAL LOWER 


" , OIJMX ( 3 ) , l STAG 

dNAME ( 1 ) ,LBLNXT ( 1 ). *NSEP( 2 ) ,SWKEF ( 1 ) , 

SIGNU) , SW( 1 ) ,DSTAR< 1) ,ODSTAR (1) 

BNAME 

( BNAME ,BDT) , { LBLNXT t LBNE XT )»(NSEP,LBZi)» 

( SWRLF ,UP ) , ( SIGN* LEDEX) , ( SW.ZBT) , ( DST AR , RBT ) , 
( ODS TAR, ANGd f ) 


0 I M EM S I ON NAVfcUfUi) 

EQUIVALENCE ( NAMEUB, ANGB T ( 4 ) ) 

DIMENSION S w T ( I JO) ,DSTAR T ( 10 j ) , DDSTRT ( 100 ) 


F = F D 

SiD = SIDD 

I F ( F . l Q . .OR. F.EQ.l.) F=BITS 

D S 1 G M A = 0. 

C SEARCH FUR MATCHING BOUNDARY NAME 

L >< = LBF(MAME) 

iriLB.EQ.u) CALL ERROR 1 

I ^ INDEX OF POINT WHICH BEGINS THE INTERVAL 

Sf I = DISTANCE FROM POINT (I) 

S F I P i = DISTANCE FROM POINT (1+1) 

MINI = EB + LBZHLB) 

I = MIN 1 + 3* < INTVL-1) 

MAXI = LB+LBNEXT ( LB ) — 12 
75 CALL B A R C ( I ) 

C IF -I- IS THE FIRST OF A DOUBLE POINTt BACK UP TO PREV INTERVAL 

IF ( SINTVL.NE.U. ) GO TO 80 
I = 1-3 

FO = 1 . 

IF l I .L T .MINI ) CALL ERROR 1 

GO TO 75 

dJ IFIFD.EG.l. .OR. S1D.GT.SINTVL) S1D=SINTVL 

SFI = DS1+S1D 
SFIPI = SFI-SINTVL 

C IS the new POINT WITHIN THIS INTERVAL 
1 )u IF(STI) 120, 114* 114 
114 IF(SFIPI) 160,160,140 

C (MOVE COUNTERCLOCKWISE) 

120 IF ( 1 .GT .MINI ) GO TO 125 
ns 1 GM A=-SF I 
SF I = 0. 

GO TO 23 
125 1 =1-3 

E = H I T S 

SFIPI = SFI 
CALL BARC ( I ) 

SFI = SFIP1+SINTVL 
GO TO 100 

C (MOVE CLOCKWISE) 

140 I F ( I.LT.MAXI ) GO TO 145 
OS 1 GM A= -SFIPI 
SFI = SINTVL 


-145 



go tu / 3 ■> 

1A5 I = I *3 

F = h I T S 

SE I = SFIP1 
CALL UARC(l) 

SFIP1 = SFI-SIMTVL 
GO TU 100 

C CALCULATE COORDINATES OF THE NEW POINT (PROPER INTERVAL FOUND) 
160 IF( F.tQ.BITS) GO TO 230 
IF(DSi) 210,220,220 
210 F = F*SFI/S1D 

GO TO 2 50 

220 F = ( (SFI-S1D)+(SINTVL-SFI)*F)/(SINTVL-S1D) 

GO TO 250 

C (NEW INTERVAL) 

230 F = SFI/SINTVL 


250 G 

RZ ONL Y = 
CALL BF 
ZD 
RD 

ANGD = 
CUR VD = 
S1DD = 


l.-F 

.FALSE. 

I 

ZBT( I ) +ZM 
R8T ( I ) + RM 
ANGCHD+ANGM 
CUR VM 
SIM 


FD = F 

INTVL * (I - (LB + LBZ1UB ) ) )/3 ♦ 1 


C ************** ********** BOUNDARY LAYER ADJUSTMENT ******************** 

I F ( LDE.NE.O .AND. PDUM( 15) .NE.O. ) WRITE (6,288) NAME, ZD, RD, 

* ANGD, CURVD, S IDD 

I F ( LDE.EO.O ) GO TO 300 
CALL GET I X 

I F ( I STAG • EQ . 1 ) GO TO 300 
LOWER = .TRUE. 

I F ( NAMEUB(L) .EO.NAME ) LOWER=. FALSE . 

LBL = LBDYBUNAME, LOWER) 

I F ( LBL.EQ.O ) GO TO 300 
NAMBL = IBLB(LBL) 

C SEARCH FOR NAMBL IN BL TABLE 

LO = LDO 

270 I F < LD.GT.LDE ) CALL ERRQR1 

I F ( BNAME(LD) .EQ. NAMBL ) GO TO 280 
LD * LBLNXT(LD) 

GO TO 270 

2 HO NVAL = (LBLNXT (LDJ-LD-6 ) /3 

LD1 = LD 

DO 281 1=1, NVAL 
SWT ( I )= SW(LDl) 

DST ART ( I )= DSTAR(LDl) 

DOS TK T ( I )= DDSTAR ( LD1 ) 

2 H 1 LD1 = LD1 + 3 
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C ; VALIJAII Va l F I IK IOTFKPOL All C3T>I 

S«I - S IGM( Lf; ) * ( BAKCS< NAME . i » INTVL ) ♦ Si DO- S WREF ( LO) ) 

IK J3FK LO) .1 Q.'J ) GU ID 28t> 

S a S i. K = SW( DGIP(ID) ) 
to R I f i ( 6, iOOi ) MAM.BL , SWSEP 
iOJl FORMAT ( /6X, an* * W A R N I rt G * *,6X, 

* 26H SEPARATED BL , BOUNDAR Y= , IX , A6 » 3X , 3HSW=, F 14.6// > 

2dd CALL LF I T 1 ( SWT,OSTART,NVAL,SWI ,DSTRCtl) 

CALL LFITK SWT,DDSTRT,NVAL,SWl,ANGC f 1) 

ANiGD = ANGD + SIGN(LD)*ANGC 
C A N l> — 0 • 

I F { .OUT. LOWER ) CANG=PI 

ZO = ZD-SIGNj{LD)*OSTRC*SIN(ANGD-CANG) 

R D = KD + SIG\HLD)*DSTRC*COKANGD-CANG) 

I F ( (' GUM ( 1 5 ) . E 0 . 0 . ) GO TO 3 JO 

WkI II ( 6 t 2 B 9 ) NAME.MAMBL ,ZD,RD, ANGD , CUR VD ♦ S 1DD , S W I ,DSTRC, ANGC 

2 db FORMAT ( / /‘3XfA6.2Xf5E16.8) 

289 FORMA T( //5X, A6,2X, A6.2X, 5F 16 . 8/ 21 X , 3E 16. 8 ) 

3 j j RCTUKxl 

F \ D 


-INSTA- 


INST A( LNEW,LBASE.L3,D0WNB,MA,MB) 

INSERT A STATION 

DOWNB 

c input- 

C LNEW = LOCATION IN STAT I ON- TABLE OF NEW STATION 

C LBASt = LOCATION OF BASE STATION 

C L3 LOCATION OF DOWNSTREAM (OR UPSTREAM) STATION 

C DOWNB = T IF L3 IS AN UPSTREAM STA, OTHERWISE =F 

C MA , MB = NEW STATION FILED POINT INDEX LIMITS 

C Z.R.PHll FIELD VALUES 

C OUTPUT- 

C LNEW = STATION FOLLOWING NEW STATION 

COMMON /ALLCOM/ M ACHA , PS A , TS A , P TA , TT A , AXI A,RGA, GAMA, 

1 MACHC,PSC,TSC,PTC,TTC, AXI C , RGC, GAMC , 

2 DAXIT,SCALEA,TTE,CHOTST 

REAL MACHA ( I ) , MACHC 

LOGICAL AXIA.AXIC 

LOGICAL CHOTST 

COMMON /CBEAM2/ DR , DZ , YPA , YPB ,F ,G, DX, YQDX , ZM , RM , ANGM ,CURVM, S 1 M, 
1 RZONLY, ANGCHD* SI NTVL , YPASQ, YPAB , YPBSQ 

LOGICAL RZONLY 

C INDEX- M=MQ , NM 

COMMON / CZ / Z ( 300 ) 

COMMON /CR / R ( 300 ) 

COMMON /CS2 / S2I300) 

COMMON /CSI / SI (300) 

COMMON /CPHI1 / PH 1 1 ( 300 ) 

COMMON /CM / JMS( 300) 

COMMON /CCURV / CURV(300) 

COMMON /CB / B ( 300 ) 

COMMON /CIDEX / M , J , MU ,MD , I S TAG 

COMMON /IXORIG/ LHO,LHE, LBDO ,LBDE , LTO » LTE , LWO,LWE, LFO,LFE, 

* LO,LESTA, LDUM ( 8 ) , 

* MOtNM, N J , NFCOL S, MAXNJ,MAXOL ,MAXNM , MAXLE , 

* LEO, LEE, LRO,LRE ,LRD 

DIMENSION LIMITS<24) 

EQUI VALENCE (LIMITS, LHO) 

COMMON /SLTAB / W ( 128 ) , X 2 ( 12 8 ) , SLCHN ( L28 ) 

INTEGER SLCHN 
C STATION TABLE 

C INDEX- L=LO,LESTA 

C S CHUK E = STATION CHOKE INDICATOR ( ADJWF ,BRHS , WRIOUT ) 

C MCL = SHARP CORNER INDICATOR ( BLDTBS ) 

C MCL = FIELD INDEX OF CONTROL STREAMLINE ( PT MOVE , FLOBAL ) 

COMMON /CHDATA/ X 1 ( 1 ) , LNE XT ( 1 ) , ML B ( 1 ) , MUB < 1 ) • PRI M ( 1 ) , 

1 TYPELBd ) , NAMELB ( 1 ) , I LB ( 1 ) ,FLB(1 ) ,S1LB( 1) , 

1 T YPEUB ( 1 ) »NAMEUB ( 1 ) , IUB(l) ,FUB(1 ) ,S1UB(1) , 

3 VMB( 1) ,DWDV< i ), X2CLU) ,VCL < 1 ) ,MCL ( 481 ) 

LOGICAL PRIM 

INTEGER TYPELB, TYPEUB 

DIMENSION SCHOKE(l) 

EQUIVALENCE ( SCHOKE , DWDV ) 

COMMON /CATAN3/ DANG 


♦DECK INSTA 

SJHRDUT INE 

* IN', T A- 

LOGI CAL 
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COMM ) „ 
CuMMON 
COMMON 
COMMON 
C n M Ml , N 
COMMON 
COMMON 
COM Ml i j 


/CBDYPT/ 
/CHITS / 
/ C V A X I T / 
/CPI / 
/CPR INT/ 
/CVM / 

/ r p a s t / 

/CI H / 


ANGD,CURVD 
HITS, [BLANK 

MAXI T.MAJCTR, GREF IN, EDUM 
PI ,TWQP1 ,PIQ2,PIQ4,T0DEG,T0RAD 
PIJUMK 3) , PRL F IN 
VM (300 ) 

ASLI BOO) 

L N » DUMCF D ( 33 ) 


I N T L G t w 
LOGICAL 


BUYNAM,F ARFLD,FREE ,F I E LD , PRE S , SOL 1 0 
UPU.UPD 


DATA F ARFLD/6HFARFLD/, F I ELD/ SHF I E LD/ , F REE /AHFRE E / * PRE S/4HPRE S/ , 
* SOL IC/5HS0L ID/ 


*** KtLOCAlE TO MAKE ROOM FOR THE NEW STATION 

INITIALIZE NEW-STATION VALUE TO THE BASE-STATION VALUES 
CORRECT THE STA-TABLE INOICIES- L-END, L-BASE, L-THREE, L-UPSTREAM 

L N = L N E W 

N^OV. = LN-i - LESTA 

LH = L BASF 

CALL MOVE ( 2 » X i ( L N ) , Xl ( LN+20 ) ,NMOVE , D, X 1 < L B ) , XI ( LN ) , 20 , 1 ) 

LESTA = LESTA+20 

LT = L3+20 

LU - LB 

IF ( .NUT .DOWNB) GO TO 60 
LH = LB + 20 

L T = L 3 

LJ = L3 

C UPDATE THE POINTERS TO THE FIELD-TABLE 
60 NPTS = MR-MA+1 
LNEXI ( LN ) =20 
CALL STTOF I (LN.NPTS) 

C*** OEFINh STAT ION-TABLE VALUES FUR THE NEW STATION 
X I ( L N ) = .5*(X1(LB) + X1(LT)) 

MLB ( LN ) = MA 
MUB(LN)=MB 
PRIM( LN)=. FALSE. 

X2CLI LN ) = 8 I T S 

C ** LOWER BOUNDARY STATION-TABLE VALUES 
M = MA 

CALL GET IX 
MX = MU 
IF(UOWNB) M X = M l) 

LX = LU 

CALL STANOIMX, LX, UPPER) 

IF IMX-ML8ILX) ) 2 10,220,250 
2'u CALL ERROR i 

C LOWER BOUNDARIES OF NEW AND BASE STATIONS ARE ON THE SAME SL 
2 0 IF ( lYPELB(LB) .EQ. FIELD) GO TJ 250 
If ( TYPELB(LH) .EQ.FARFLD) GO TO 260 

C FRLl " UUNUArt Y 

IF (1 YPELB(LH) .NE.FREE .AND. I YP ELB ( L T ) . NE . F RE E ) GO TO 224 
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TYPEL B ( LN ) =FREE 
GO TO 260 


C PRESSURE BOUNDARY 

2c<* IF< TYPELB(LB) .NE.PRES .ANO. TYPELB(LT) .NE. PRES) GO TO 230 
T Y P E L B ( LN)=PRES 
GO TO 260 

C SOLID BOUNDARY 

230 TYPELB(LN)=SOLIO 
BDYNAM= NAMELB(LX) 

NAMELBl LN )=BDYNAM 
ILB(LN)=ILB(LX) 

FLB(LN)=FLB(LX) 

S L L B ( LN)=S1LB(LX) 

LD = LU 

CALL STANO(MU»LU»UPU) 

CALL STANO(MD,LD,UPD) 

DS1 = . 5* ( BARCS ( BDYNAM , ILB ( LU ) * ILB(LO) ) • S1LB ( LD ) -S1LB ( LU!) 

I F ( UPU.OR.UPD) CALL ERROR l 
IF ( 1)0 wNB ) DS 1 =-DS I 

CALL BDYPTM(BDYNAM,ILB(LN),Z(M) ,R(M) ,FLB(LN),S1LB(LN) ,DS1 »GMA ) 

IF < GMA.NE .0. ) CALL ERROR 1 
PHI i( M ) = ANGD 
B ( M I = . 5* ( ft ( MU ) +B ( MD ) ) 

VM ( M ) = .5*(VM(MU)+VM(MD) ) 

IF ( VM( M ) .EQ.O . ) VM(M)=VMIMU+1> 

GO TO 300 

C INFIELD BOUNDARY 
ZbO TYPELB(LN)=FIELD 
ISTAG =3 
CALL SAVIX 
NAMELBl LN )= IBLANK 
260 ILB(LNI)=0 

FLB(LN)=BITS 
S1LB< LN) = BITS 

C ** UPPER BOUNDARY STATION-TABLE VALUES 
300 M = MB 

CALL GETIX 
MX = MU 

IF(DUWNB) MX = MD 
CALL STANOIMX, LX, UPPER) 

IF(MUBILX)-MX) 310,320,350 
31j CALL ERROR l 

C UPPER BOUNDARIES OF NEW ANO BASE STATIONS ARE ON THE SAME SL 
320 IF ( TYPEUB(LB) .EQ. FIELD) GO TJ 350 
IF( TYPEUB(LB) .EQ.FARFLD) GO TO 360 

C FREE BOUNDARY 
LD = LU 

CALL STANO(MU,LU,UPU) 

CALL STANOI MD,LD,UPD) 

IF (TYPEUB(LB).NE.FREE .AND. TYPEUBi LD) .N6 .FREE) GO TO 324 
TYPEUBI LN )=FREE 
GO TO 360 



C PRESSURE BOUNDARY 

324 IF (TYPEuBILR ) .NE .PRES .AND. T YPE UB I LD ) . NE . PRE S ) GO TO 330 
TYPEUBl LN)=PRES 
“ GO TO 360 

C SOLID BOUNDARY 

- 3 3j TYPtUB(LN)=SOLlU 

BDYNAM= NAMEUB(LX) 

NAMtUB ( LN ) = BDYN AM 
I U B ( L N ) = I U B ( L X ) 

F U B ( L N ) = FU0 ( L X ) 

S 1 U B ( LN) = S1UB(LX) 

LP = LU 

CALL STANO(MU,LU,UPU ) 

CALL STAND ( MD,LO,UPD ) 

IF ( .NUT.UPU .OR. .NOT.UPD) CALL ERROR I 

CSI = . 3*( BARCSI BDYNAM , IUB (LD ) , IUB(LU) ) + SI UB ( LU ) -SIUB ( LD ) ) 

IF (.NOT. DOWN B) DSi=-DSl 

CALL BDYPT M ( BDYNAM ,IUB(LN)»Z(M),R(M),FUB(LN) , SIUB ( LN ) » OS 1 » GMA ) 

IF ( GYA.NE .0. ) CALL ERR0R1 
• PH 1 1 ( M ) = ANGD-PI 

HIM) = . 5* ( B ( MU ) +B( MD ) ) 

VM(M) = .5*(VM(MU)+VM(MD) ) 

T IF ( VM( M) .EQ.O. ) VM(M)=VMtMU-U 

GO TO 400 

- C INFIELD BOUNDARY 

3 60 TYPEIJB(LN)=FIELD 

I STAG = 3 

_ CALL SAVIX 

NAKLUBI LN) = I BLANK 
360 I U B ( L N ) = 0 

FUB ( L N ) = B I TS 

- SiUBI LN ) = B I TS 

C DEFINE THE FIELD POINTS BY CUBIC POLYNOMIAL INTERPOLATION ON SL-S 

4 00 M = M A 

RZ ONL Y= .TRUE. 

I F ( TYPELB(LN) .EQ. SOLID) GO TJ 420 
4 1 j CALL GET I X 

DZ = ZIMD)-Z(MU) 

DR = R(MD)-RIMU) 

F = .6 

- G = .5 

ANGCH D= ATAN3(DR,0Z,PHI1(MU> > 

YPA = PHI i ( MU J-ANGCHD 
YPB = PHI 1 ( MD I-ANGCHD 
MSV = M 

MU S V = MU 

. MDSV = MD 

~ M = MD 

CALL GETIX 
1 S T AGD= ISTAG 

MD = M 

M = MS V 

MU = MUSV 

IF ( I STAGD.EQ. i ) YPB=-YPA 


RZONL Y = .FALSE. 

CALL BF I 

Z < M ) = Z (MU) +ZM 

R ( M ) = R ( MU ) + RM 

PH I 1 ( M )=ANGCHD + ANGM 
VM(M) = F*VM(MD)+G*VM( MU ) 

B ( M ) = F*B(MD)+G*B(MU) 

C CHECK FOR POINTS ON A SLIP LINE 

IF(M.IQ.MA .OR. WUJ.NE.O.) GO TO 420 
Z ( M ) = . 5* ( Z ( M- 1 ) +Z ( M ) ) 

M = M-I 

CALL GET I X 
M = MSV 

DZ = .25*<Z(MUSV)-Z(MU)+Z(MDSV)-Z(MD) ) 

DR = ,25*(R(MUSV)-R(MUM-R(MDSV)-R(MD) ) 

Z ( M - L )= Z ( M ) - DZ 
R( M-l ) = R ( M ) -DR 
Z(M ) = Z ( M ) +DZ 

R ( M ) = R ( M ) +DR 

A20 M = M+l 

IF(M-MB) 410,425,500 
425 IF( TYPEUB(LN) .NE. SOLID) GO TO 410 

C CHECK FOR OUT-OF-ORDER POINTS 
5 JO NORD£R= 0 
502 NORDER= NORDER+1 

I F ( NORDER . GE . 20 ) CALL ERR0R1 

MX 1 =0 

MAPI = MA+1 

MSV = MA 

S2(MA)= 0. 

DO 520 M=MAP 1 »MB 
DR = R(M)-RIM-l) 

DZ = Z (M)-Z(M-l ) 

S2(M) = S2(M-1 ) +SORT (DR*DR+DZ*DZ ) 

CALL GET I X 

IF ( W( J ) .EO.O. ) GO TO 518 
ANG2 = AT AN 3 < DR ,DZ»PHI1 ( M— 1 ) ) 

ADANG = ABS(DANG-PIQ2) 

IF(MXl.NE.O) GO TO 515 
IF ( ADANG. GE.PIQ2) MX1=MSV 
MSV = M-l 

515 IF( ADANG. GE.PIQ2) MX2=M 
GO TO 520 

516 IF( (M-1J.EQ.MX2I MX2=M 
520 CONTINUE 

C DEFINE THE FIELD PT LOCATIONS BY UPSTREAM AREA DISTRIBUTIONS 
IF(MXl.EQ.O) GO TO 999 
MX 1 = MAX0(MX1-N0RDER,MA) 

MX 2 = MIN0( MX2+N0RDER, MB ) 

WRITt (6, 1550) MX 1 ,MX2 
1550 FORMAT ( 14H IN ST A-MX 1 , MX2 , 2 1 6 ) 

MX 1 = MAXO ( MX1- 1 ,MA ) 

MX 2 = M INO ( MX2+ 1 ,MB ) 

C ADD UP UPSTREAM AREAS 

M = MX1 

CALL GET I X 

K = 1 
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ASi ( ; ) = u. 

5 62 MJMI = MU 
M = M + l 

K = K+ 1 

CALL UfcT I X 

AREA = SQRT( { R ( MU )-R ( MUM L ) ) * ( R ( MU ) -R ( MUM1 ) ) ♦ 
i (Z(MU>-Z(MUM1))*(Z(MU)-Z1MUM1) )) 

I F ( AX I A ) AREA= ( R ( MU) +R ( MUM1 )) *AREA 
AS L ( K ) = ASL(K-1 )•* AREA 
IFIM.LT.MX2) GO TU 562 
ASLMK = A SL ( K ) 

C INTERPOLATE FOR COORDINATES 

OZBA = Z ( MX 2 )-Z( M X 1 ) 

ORB A = R ( MX2 )-R I MX1 ) 

DRS<*HA= DRBA*(R(MX2)+R(MX1) ) 

RMASQ = RIMX1 )*R( MXl ) 

DVMBA = VM(MX2)-VM(MX1) 

M = MXl+1 

K =2 

564 F = ASL ( K ) / ASLNK 

Z(M) = Z(MXI)+F^DZBA 
R ( M ) = R ( MX I ) + F*DRB A 

IF ( AX I A) R(M)=SQRT(RMASU+F*DRSQBA) 

VM(M) = VM ( MX 1 ) +F ♦DVMBA 
M = M + l 

K = Ml 

IFtM.LT.MX2) GO TO 564 
GO TO 502 

94V LNEw = LN + 25 
RETURN 
END 


• — 


o n o o o 


♦DECK ERRORM 

SUBROUTINE ERROR 1 
CEDUMPM EDUMP FOR 


STCM LINK 


TABLE OF INDEX LIMITS 
COMMON /IXORIG/ LHO.LHE. 

♦ 

* 

♦ 


LBDOtLBDEt LTO,LTE * LWO.LWE, LFO,LFE, 

LO.LESTA , LS0,LSE,LD0,LDE,LDUM(4) , 

MO ,NM, NJtNFCOLS* MAXNJ » MAXOL * MAXNM , MAXLE , 

LEO, LEE* LRO,LRE»LRD 
LIMITS(24) 

(LIMITS, LHO) 


DIMENSION 
EQUIVALENCE 
STREAMLINE TABLE 

COMMON /SLTAB / W ( 128 ) ,X 2 t 128 ), SLCHN ( 1281 
INTEGER SLCHN 
STATION TABLE 
INDEX- L=LO,LESTA 

SCHOKE= STATION CHOKE INDICATOR ( ADJ WF , BRHS ♦ WR IOUT ) 

MCL = SHARP CORNER INDICATOR (BLDTBS ) 

MCL = FIELD INDEX OF CONTROL STREAMLINE (PTMOVE,FLOBAL) 

COMMON /CHDATA/ X 1 ( 1 ) , LNE XT ( 1 ) , MLB ( 1 J , MUB ( 1 ) , PRI M ( 1 ) , 

1 TYPELB ( 1 ) tNAMELB ( 1 ) , ILB(l) ,FLB(1 ) ,S1LB( 1 ) , 

1 TYPEUBI 1 ) tNAMEUB ( 1 ) , IUB(1) ,FUB(1 ) , S1UB( 1 ) , 

3 VMB( 1) ,DWDV( 1 ), X2CLU) ,VCL (1) ,MCL< 481 ) 

LOGICAL PRIM 

INTEGER TYPELB, TYPEUB 


DIMENSION 


SCHOKEd ) 

EQUIVALENCE 


( SCHOKE, DWDV ) 

COMMON 

/CA2 

/ 

A2 ( 300 ) 

COMMON 

/C A3 

/ 

A3 ( 300 ) 

COMMON 

/CA4 

/ 

A4 ( 300 ) 

COMMON 

/CA5 

/ 

A 5 ( 300 ) 

COMMON 

/CA6 

/ 

A6 ( 300 ) 

COMMON 

/CA7 

/ 

A7 ( 300 ) 

COMMON 

/CA8 

/ 

A 8 ( 300 ) 

COMMON 

/CB 

/ 

B( 300) 

COMMON 

/CCURV 

/ 

CURV( 300 ) 

COMMON 

/CDS2 

/ 

DS2( 300) 

COMMON 

/CDDS2 

/ 

DDS2 

COMMON 

/CFB 

/ 

L»MA»MB»PLB»PUB,WF , CHOKE , SUBSON 




XCHOKE , TAREA »VMBC , WRQST.WCALC 




J SUM , VML BSQ 

LOGICAL 


CHOKE, SUBSON 

COMMON 

/C I DEX / C 1 ( 5 ) 

COMMON 

/CIDEXR/ 

C2 ( 25 ) 

COMMON 

/CPHI1 

/ 

PH 1 1 ( 300 ) 

COMMON 

/CR 

/ 

R ( 300 ) 

COMMON 

/CRHS 

/ 

RHS( 300) 

COMMON 

/CS1 

/ 

S 1 ( 300 ) 

COMMON 

/CS2 

/ 

S2 ( 300 ) 

COMMON 

/CTABPR/ 

1 1TAB 

COMMON 

/CTOLRL/ 

C 3 ( 12 ) 

COMMON 

/CVM 

/ 

VM ( 300 ) 

COMMON 

/Cl 

/ 

Z( 300) 

COMMON 

/BLBDY 

/ 

I BLB ( 60) 


NK, PLBC, PUBC , 
QV ( 8 ) * QVP ( 8 ) » 


CALL TABPRT(3HCFB,L,33,4) 
CALL TABPRT ( 5HC I DEX ,C 1 ,5,5) 
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CALL TABPRT l 6HC I UtXR , C2 , 25 , D ) 

CALL IABPRT (6HCTGLRL,C3,6,6) 

I i I A 1 = LO 

CALL IABPRT ( 6EISTATAB, Xl ,LESTA, 5) 
f- CALL JM jPR I 

CALL IABPRT I <rMS 1 , S 1 , NM, 1 ' ) 

CALL r ABPR T ( 2HS2 » S2 , NM , 1 0 ) 

*_ CALL IABPRT ( 1HZ , Z ,NM, 10 ) 

CALL IABPRT ( 1HR , R ,NM , 10 ) 

CALL TABPRT(4HCURV,CURV,NM,1G) 

CALL TABPRT ( 2HVM , VM , NM, 10 ) 

WRITE (6,1000) 

DO 100 1=1, NM 

WRITL (6,1001) I,B( I) ,A2 ( I ) ,A3( I ) , AM I) ,A5(I ) , A6 ( I) ,A7( I ) ,A8( I ) , 
1 DS2 ( I ) ,RHS ( I ) 

loO CUNTlNUt 

WRITE (6, 1002) DO S2 

10OU FORMAT (4H1 M , 1 1 X , 1HB , 1 0 X , 2HA2 , 10X , 2HA3 , 1 OX , 2HA4 , 10X , 2HA5 , 1 OX , 

1 2HA6, 10X,2HA7,10X,2HA8,9X,3HDS2,9X,3HRHS) 

1001 FORMAT ( 1H , I 3 , 8F 1 2 . 3 , 2F 1 2 . 6 ) 

1002 FORMAT ( / / / 8 H DS2 M X= , F 12 . 6 ) 

* I F ( IBLBID.NE.O ) CALL TABPRTt 5HBLBDY, IBLB,60,3 > 

I F ( LDE.EQ.O ) GO TO 1321 
*_ I1TAB = L DO 

CALL TABPRT (6HBLT AB, CHNAM ,LDE ,3 ) 

1321 CONTINUE 
LSTOP = 5 

GO TO (999,999) , LSTOP 

9 99 RETURN 
END 
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